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Abstract—The transient rates of absorption of carbon dioxide into water, non-reactive solutions, 
sodium hydroxide and alkaline buffer solutions have been measured at 23°C. The liquid was 
carried through the gas in the form of a thin film ona rotating drum. the time of exposure varying 
from 0-01 to 0-25 see. The liquid was free from convection currents, so that the rates of absorption 
were controlled by diffusion and chemical reaction in the liquid. The experimental results are 
analysed in this light and values obtained for solubilities, diffusivities and reaction velocity 
constants in the various solutions. 

The relationship between the rate of absorption into reacting solutions and the homogeneous 
reaction kinetics is in general not a simple one, because of the depletion of reactant near the 
surface and other factors. Caution should be exercised against oversimplification when pre- 
dicting or explaining the results of gas absorption under industrial conditions. 


Résumé — Les vitesses d’absorption, en régime transitoire, du bioxyde de carbone dans l'eau, 
des solutions non réactives, des solutions caustiques et des solutions tampon alcalines ont été 
mesurées & 25°C, Le liquide était transporté 4 travers le gaz sous forme d’un film mince sur 
un tambour tournant. La durée du contact variait de 0,01 & 0,25 sec. Les courants de con- 
vection ctaient évités dans la phase liquide, de telle maniére que les vitesses d’absorption étaient 
controlées par diffusion et réaction chimique dans la phase liquide. Les résultats expérimentaux 
sont analyses a partir de ces considérations et les constantes de solubilité, de diffusion et de 
vitesse de réaction sont obtenues pour différentes solutions. La relation entre la vitesse dabsorp- 
tion dans les solutions réactionnelles et la cinétique des réactions homogénes n'est pas simple 
en général, & cause des diminutions de concentration en surface, ainsi que de influence d'autres 
facteurs. Il faudra prendre des précautions quant aux simplifications utilisées quand on prédit 
ou explique les résultats de absorption des gaz dans des conditions industrielles. 


Zusammenfassung Die nichtstationiren Geschwindigkeiten der Absorption von Kohlendioxyd 
in Wasser, nichtreaktiven Losungen, Natriumhydroxyd und alkalischen Pufferlsungen wurden 
bei 25 C gemessen. Die Fliissigkeit wurde durch das Gas in Form eines diinnen Films auf einer 
rotierenden Trommel geleitet, die Kontaktzeit lag zwischen 0,01 bis 0,25 sec. Die Flissigkeit 
war frei von Konvektionsbewegungen, so dass die Geschwindigkeiten der Absorption durch 
die Diffusion und die chemische Reaktion in der Flissigkeit bestimmt wurden. Die Versuchser- 
gebnisse wurden in dieser Hinsicht analysiert und es wurden Werte fiir die Léslichkeiten, die 
Diffusionskoeffizienten und die Geschwindigkeitskonstanten der Reaktion in den verschiedenen 
Ldsungen erhalten. 

Die Bezichung zwischen der Geschwindigkeit der Absorption in reagierende Lésungen und 
der homogenen Reaktionskinetik ist im allgemeinen nicht einfach wegen der Erschépfung des 
Reaktanten in der Nihe der Obertliiche und wegen anderer Faktoren. Die Vorginge diirfen 
daher nicht allzu stark vereinfacht werden bei der Vorausberechnung oder der Deutung von 
Ergebnissen der Casabsorption unter industriellen Bedingungen. 
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Tur experiments described here were carried out 


in connection with an investigation of the 


mechanism of gas-absorption in packed columns 


which has already been the subject of a per 
liminary report [1], and will be further reported 
elsewhere. An interpretation of the results in 
1 vsico chemical terms ilthough not necessary 


for the purpose of the origmal investigation, may 


be of some interest ind the technique used may 


find a wider application in the study of fast 


liquid phase reactions 
the rate of ab 


sorption of carbon dioxide into various aqueous 


It was necessary to measure 
solutions in the first trae tion of a second of contact. 
\ number ol pre VIrous worke rs have had the sate 


reviewed by Cun 


objective methods, which were recently 


and Davipso 2), have 


probably all been subject to error, Since the work 


deseribed mm thus paper Ww is done, the latter authors 


2 3) have carned out meticulous measurements 


of the absorption 1 gases into jets and into water 


running over the surface of a sphere. The 
shortest exposure time attained was about 0-0] 
se" is with the rotating drum described here. 


The internal consistency of their results suggests 


that the methods used were highly reliable, and 


the value for the diffusivity of carbon dioxide in 


water at 25 C deduced from their work (namely 
Bl 10 ci”. se has been adopt i here as 
ben a probably the most trustworthy, in the 


literature, The rm iting drum inically 


mere compli ated than the falling film or jet 
used by Cutten and Davipson, but it vields 
results which are more easily interpreted when 
the absorbed gas reacts with the liquid, as the 


complication introduced by the “ stretching ” of 
uidl surface is avoided. 

The compositions of the aqueous solutions used 
were dictated by the requirements Of the parallel 
series of experiments on absorpti m m packed 
columns. It was necessary for the viscosities and, 
as tar as possibl . the densities of all solutions to 
be the same. In the case of sodium hydroxide 
this was achieved by adding varying concentra- 
tions of sodium sulphat: to the solutions. The 
sodium hydroxide solutions also contained some 
carbonate, for reasons connected with the packed 


column experiments. 
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EXPERIMENTAI 


A diagrammatic view of the apparatus used 
for the principal series of experiments is shown 
in Fig. 1. The drum, of stainless steel, was 11-4em 


dia. and 12-4em long. The drum surface was 


ground smooth to give a cylinder truly con- 


The drum itself was 


centrie with the spindle. 
inside a Perspex box, through which the spindle 
stuffed 


asbestos. The spindle was supported externally 


passed via glands with graphitised 


on ball bearings and driven by an electric motor 


through a system of gears. 


Gas 
A te 
Vent djustable picte 
= 
| i 
Mercury ¥ 
seo! 
Liquid 
outiet 
Fic. 1. Rotating drum apparatus. 


In order to separate the inflowing and outflow- 
ing liquid, the lower part of the drum dipped into 
a trough of mercury, and a rubber roller bore 
against the drum so as to squeeze off any film of 
liquid which might otherwise have been carried 
through the mercury. The gap between the end 
faces of the drum and the sides of the Perspex 
box was so narrow that the rate of liquid leakage 
was negligible compared to the rate of transfer of 
Two horizontal 

the 
These plates could be set at 


liquid by the intended route. 


plates separated the gas from horizontal 
liquid surface. 
different heights, and each carried an adjustable 
plate which could be slid forward to cover the 
gap between the edge of the fixed plate and the 
drum. 

The liquid level on either side of the drum was 
adjusted to conform with the level at which the 
horizontal plates were set. It was thus possible 
to expose varying proportions of the drum’s 


surface to the gas. The adjustable horizontal 


The kinetics of absorption of carbon dioxide into neutral and alkaline solutions 


plates were knife-edged, and by means of feeler 
gauges they could be aligned close to and parallel 
to the drum surface. 

The drum speed could be varied by eight steps 
from 1-9 to 5-l rev sec. By varying also the 
fraction of the surface exposed, the same exposure 
time could be obtained with different drum speeds. 

The thickness of the liquid film on the drum 
must be suflicient for it to be in effect “infinitely 
deep”. It is easily shown from the diffusion 
equations [4] that the finite depth, 1, of a liquid 
layer will not affect the amount absorbed in time 
t by more than } per cent providing Dt/P < 0-25, 
where D is the diffusivity of the dissolved gas. 
For carbon dioxide in water, this corr sponds to 
l>4 10° em for ¢ 0-25 sec. The film thick- 
ness was kept well above this value ; it could be 
controlled by varying the clearance of the 
adjustable knife edge on the inlet side, and 
determined by measuring the rate of carryover 
of liquid by the drum. The gap between the 
knife edge and the drum was in practice set at 
5-7-5 10%em, and variation within these 
limits did not affect the absorption rate. (The 
thickness of the film was greater than the width 
of the gap.) 

With the drum running, the knife edge on the 
outlet side was advanced until it touched the 
liquid surface, a fine ripple being formed parallel 
to the knife edge; entrainment of gas on the 
outlet side then stopped. The actual width of 
the gap was not critical. 

The box was closed by a lid which rested in a 
mercury seal, The box was provided with a gas 
inlet and a vent for purging the gas space, and 
with pipes for the entry and exit of the liquid. 
Although the water used for the solutions was 
de-aerated by boiling under vacuum, it was 
found that the rate of absorption gradually fell 


off after purging, due to the accumulation of 


insoluble gases in the box. 

The carbon dioxide was held at atmospheric 
pressure in a large balloon, and allowed to flow 
into the box through a soap-bubble flow-meter. 
The carbon dioxide was taken from a cylinder, 
and contained 0.3 per cent of relatively insoluble 
gas. However, in some comparative runs using 
carbon dioxide from dry ice, and containing less 


than half as much impurity, the rates of absorp- 
tion were not significantly different. The carbon 
dioxide was humidified by passing it through 
water at 25 C before storage in the balloon. 

If the partial pressure of water at 25°C is p and 
the carbon dioxide is saturated with water- 
vapour, the measured rate of passage of carbon 
dioxide and water-vapour through the flow-meter 
will be [P (P 


rate of absorption of carbon dioxide and P the 


p)| q moles sec, where q is the 


total pressure. Now in general, subject to a 
qualification mentioned below, theoretical con 
siderations indicate that the rate of absorption 
into a particular solution is proportional to its 
physical solubility. Since the solubility obeys 
Henrv’s law under the conditions of these 
P — p)/P|C*, 
where C* is the solubility at a partial pressure P. 

Thus the effect of the dilution of carbon 


experiments, it will b equal to 


dioxide with water vapour was to increase its 
apparent bulk and decrease its rate of solution 
in the same ratio, so that the observed volumetric 
rate of passage of moist carbon dioxide was 
equal to the rate at which dry carbon dioxide 
would have been absorbed at the same _ total 
pressure P, All experimental results were 
corrected to a standard total pressure of 760 mm, 
and the whole apparatus was housed in a constant 
temperature room at 25°C. 

In interpreting the experimental results it 
is assumed that the liquid surface moves at the 
speed of the drum, and that end effects can be 
ignored. This assumption was justified by the 
fact that for a given exposure time there was 
generally no systematic difference between 
results obtained with different drum speeds, 
exposed areas and film thicknesses. However, 
at the lowest drum speeds it was noticed that the 
capillary ripples formed by the knife edge at the 
outlet side sometimes had their origin several 
millimetres upstream of the knife edge. Particles 
of tale scattered on the surface collected in the 
zone between the knife edge and the ripples, 
where the the surface seemed to be stagnant. 
This effect is undoubtedly due to the accumula 
tion of surface active contaminants, and has been 
noticed elsewhere [5]. This is a potential source 


of error, but precautions were taken to keep the 
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stagnant zone, in the worst cases, down to about 
2 per cent of the exposed surface. 


Technical grade sodium hydroxide, sodium 


carbonate and sodium bicarbonate were used 


because of the large quantities required for the 
The 


sodium 


packed column experiments. 
of the 
buffer solutions was determined by 


‘ omposition 


sodium carbonate bicarbonate 
titrating a 
sample with N/10 hydrochloric acid to an end- 
point with brom-—cresol green, thus determining 


({NaHCO,] + 2 [Na, CO,)). 
a known excess of sodium hydroxide was added, 


To another sample 


the carbonate in the solution precipitated with 
barium chloride, and the remaining excess sodium 
titrated with N/10 


acid, using equal amounts of phenolphthalein 


hydroxide hvdrochloric 


this 


indicator : 
found 


and thymol blue as 


In 


rave 


practice it was best to 


neutralize almost all the remaining excess of 
sodium hydroxide with acid before adding 


barium chloride. 


RESULTS AND Discussion 


The results given here supersede some reported 


earlier [1], which obtained with a 


a" a? 


were more 


primitive apparatus, and before the possibility 


of error from the “ stagnant region” had been 
appreciated. 

Suppose the drum has a width w and its surface 
a linear velocity v, and the length of the pe riphe ry 
Then the 


it has an 


vas 18 l. 


of the drum exposed to the 
strip of liquid surface formed in time 
area wv 
length of time l/v. 


. ot, and it is exposed to the gas for a 
If the amount of gas absorbed 
by unit area of surface in time ¢ is Q{t}, the 
M.Qil vt. 
the 
amount of gas absorbed in unit time by the total 


amount absorbed by the strip is w: 


Since 1 /8¢ such strips are formed in unit time, 


exposed area of the drum is wv . Q {1% 

Absorption into water. If the surface is saturated 
for t 
C*, the amount Q absorbed in time ¢ by unit area 
is 2C*4/ Dt), where D is the diffusivity of the 
dissolved gas. A plot of @ against y! should then 


- 0 with carbon dioxide at a concentration 


be a straight line of slope 20 "4 Din. The 
experimental results are plotted in Fig. 2. The 
full line is the calculated relation, taking 
C* = 3-37 x 10 * g mole ml [6] and D= 
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5 8 *— Full line calculated for 
L 

°g mole/m| 


Fie. 2. Absorption of carbon dioxide in water. 


The the 


found by Davipson and [2, 3), 


1-92 10° em? sec, latter is value 
which is 
considered more reliable than other values to be 
It will be seen that the 
observed systematically lower than 
the line. It 


suggested that this might be caused by a very low 


found in the literature. 
results are 
calculated has previously been 
dioxide into 


DAVIDSON [2] 


“invasion coefficient” of carbon 


water. However, CULLEN and 
absorbed carbon dioxide into jets of water, using 
exposure times down to about 0-01 sec, and found 
The rise in surface 
the 


carbon dioxide and its associated water-vapour 


no trace of such an effect. 


temperature caused by condensation of 
can be calculated [7] and is much too small to 
account for the discrepancy; in any case, this 
effect would also have been operative in CULLEN 
and Davipson’s experiments. The effect is no 
doubt a systematic experimental error; this is 
probably caused by end effects, such as the 
initial adjustment of the liquid velocity to that of 
the drum or “ stagnant” surface near the exit. 
The magnitude of the error gives a measure of 
the reliability of the other results, which cannot 
be so accurately checked, 

Absorption into inert salt solutions. The results 
for absorption into 1-05 M sodium sulphate and 
0-72 M magnesium sulphate solutions are shown 
in Fig. 3. The solubilities of carbon dioxide in 
these solutions are known [8] so the diffusivities 
can be calculated. For sodium sulphate, with 
c* 1-7 10° ¢mole/ml. D is 
1-1 x 10° em?/sec. 


found to be 
For magnesium sulphate, 


q 
-+ 
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vt 


Absorption of carbon dioxide in 1-05 M sodium 
sulphate and 0-72M magnesium sulphate. 


with C* 2.2 x 10° g mole ml, D is 1-2 Ww? 
em* sec. When estimating diffusivity it is often 
assumed [9] that Dy is constant for a series of 
aqueous solutions at a particular temperature. 
Since the viscosities of both solutions were 1-57 
times that of water, this rule would predict the 
diffusivity in both solutions to be (1-92 10° 
1-57), or 1-22 


subsequent calculations that D had the value 


10° em*/sec. It is assumed in 


10° em? /sec in all other solutions used, 
since these all had the same viscosity. 
Absorption into sodium hydrovide solutions. 
Typical plots of Q against t and 1/t, for 0.085 M 
and 1-5 M sodium hydroxide solutions, are shown 
in Figs. 4 and 5. In general Q is initially propor- 
tional to ¢ and later becomes linear in y/t. The 
intercept on the \/t axis increases with increasing 
concentration. 


J 


O 00501 OD O2 

t sec 

Fia. 4. Absorption of carbon dioxide in 1-50 M and 0-085 M 
sodium hydroxide (@ v. t). 


In the earlier stages of absorption carbon 
dioxide is diffusing into a solution containing an 
almost uniform concentration of OH’ ions, with 


which it undergoes the reaction : 


(b)O°O85 M Sodum 
,. Hydroxide 


(a) M Sodium-# 
Hydroxide | 7 


g mole/c 


10 Q 


> © 
| 


Absorption of « arbon dioxide in 1°50 M and 
0-085 M sodium hydroxide (Q v. 4/2). 


CO, + OH’ HCO,’ , (1) 


which is first-order with respect to CO, and to 
OW. This reaction 1s followed by 


HCO,’ + OH’ + CO,” + HO, (2) 


which has a much higher velocity constant than 
reaction (1). Thus, provided the concentration 
of OH’ is sutfliciently high (> 0-1 M) for the 
equilibrium concentration of HCO,’ to be neg- 


ligible, the overall reaction is 
CO, + 20H’ + CO,” + H,0, (3) 


being first-order with respect to CO, and to OH’, 
and having the velocity constant of reaction (1). 
So long as the OH’ ions are not substantially 
depleted in the neighbourhood of the surface, 
the reaction will be pseudo first-order with a 
k’ [OH’). The quantity 
absorbed in time ¢ by diffusion accompanied by 


\ elocity constant k 


chemical reaction is given by [10] 


D kt 
( ce t+ thje kt + 4 
[(kt + 4) erf 4 }. (4) 


When kt 


- 2-5, @ can be closely approximated by : 
D 
a! k (kt + 4) (5) 


As the order of magnitude of k’ is 107 ml/g mole 
sec [11,12], the value of kt is greater than 2-5 for 
the most dilute solution used, 0-045 M, and for 
exposure times greater than 0-01 sec. One would 
therefore expect the Qv.t curves to be linear 
with a slope of C* V/kD, until depletion of OH’ 
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ions near the surface disturbs the first-order 
character of the reaction. 

In the case of 1-5M sodium hydroxide the 
expected order of magnitude of k is 1-3 lo” 
10". or 1-5 10* sec”! the inte rcept of equation 
(3) on the Q axis would therefore be very small, 
and in effect the equation would represent a 
straight line through the origin. As may be seen 
from Fig. 4 it is possible to draw a str ught line 
through the origin and the experimental points 
for times less than 0-1 sec. - The slope of the line 
is 7-5 1? ge mole em* sec. If the solubility of 
carbon dioxide in this solution is assumed to be 
1-7 10°° » mole ‘ml (estimated by the method 
of vAN Kreveten and Horriuzer [13] for 


sodium hydroxide solution with an ionic strength 


of 2-25) and its diffusiy ity 1-2 10° ‘see. 
k is found to be about 1-6 10° sec, whe nee k' is 
about 10° ml/g mole see. h is the 


expected order of magnitude. The results with 
other solutions, down to 0-15 M sodium hvdroxide 
vield apparent values of k’ within 20 per cent 
of this. 

Although at first sight this analysis seems 
plausible enough, it is shown in Appendix I that 
if the value of k were in fact 1.6 10* sec! in 
the 1-5 M sodium hydroxide solution, the con- 
centration of OH at the surface would be 
depleted to the extent of 10 per cent within 
about 1-5 10°“ sec; in the more dilute solu 
tions the fractional rats of de ple tion is somewhat 
more rapid, Since the « xperimental observations 
refer to exposure times of 0-01 sec and more. it 
must be supposed that the depletion of OH’ ions 
renders the above analysis invalid, and that the 
true value of k’ is substantially larger than 107 
ml g mole sec. Only if the effective diffusivity 
of OH’ ions was very much greater than that of 
dissolved carbon dioxide would the depletion be 
sufficiently small to be ignored. This is not the 
case, as is shown below. 

If the reaction between the dissolved gas and 
the reactant in the solution were instantaneous, 
the rate of solution would be controlled by the 
diffusion of CO, and OH’ to a reaction site 
situated at a depth below the surface which 


increases with time. The concentration distribu- 


tion in this case is shown in Fig. 6. The quantity 
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Fic. 6. Instantaneous reaction. Concentration distribu- 
tion at time ¢, 


of carbon dioxide absorbed in time ¢ would then 


be given by fl 


Dt 
@=2aC* / (6) 
N 
where 
x= L/erf(B \/D) (7) 
and 8 is defined by 
por [8° D) erf [B/y D| 
D exp D ] erfe \ D (8) 


where C" is the concentration of the reactant, and 
n the number of moles which react with one mole 
of dissolved gas. When D DD, 


x= (1 + C’/nC*), (9) 


A rigorous treatment of the kinetics of absorp- 
tion of a gas which undergoes a second order 
reaction of finite speed with a reactant solute 
involves the solution of non-linear differential 
equations, Perry and Picrorp [15] have used 
“ computer to obtain solutions covering a certain 
range of parameters in the special case where 
D = D’. Their results are not expressed in the 
form used here, but show implicitly that after 
a suflicient time of exposure Q approaches the 
form 


Q = 2C* (1 =<) J — A, (10) 


Here A is independent of t. It may be expected 
that, by analogy, when D ¢ D’, Q will approach 
the form : 


B=S yt—B, (11) 


Carbon Dioxide 
zone OW’ zone 
VOL 
8 
195 
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with B independent of t, « having the value given 
by equations (7) and (8). This view is supported 
by dimensional analysis, which shows that the 
relation between Q and ¢t may be put in the form : 


As k’ is increased indefinitely the behaviour 
must become that of the “ instantaneous ” 
reaction, and the left hand side becomes equal to 
afor any value of t, however small. It follows 
that for any value of k’ the left-hand side becomes 
equal to « for sufliciently large values of t, so that 
u plot of Q v. \ t will approach asymptotically 
a straight line of slope S. 

It can be shown that when « > 1 equations 
(7) and (8) become : 


(13) 
~ nC* D 


whence from equation (11) : 


gmt /P 14) 
n ~ 7 Db’ 


where S is the limiting slope of the Q v. yt curve. 
The error in « is about 1 per cent when « = 10. 
Experimental values of S are plotted against 
the sodium hydroxide concentration C’ in Fig. 7. 
Assuming C* ,D and D’ to have the same 
value in all the solutions, the curve should tend 
to a straight line of slope y/ D’ ‘x (since n = 2) 


D 
and intercept 2C* J 
7 D 


The limiting 


77 | | 
< / | 
$30 | 
4 
| 


g mole iri 


Fic. 7, Absorption of carbon dioxide in sodium hydroxide 
(Sv. [OH’)). 
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slope of Fig. 7 is 2-8 
dD 2-5 
Using this value of D’, S can be calculated from 


10° em sec?, whence 


10° em? /sec. 


equations (7), (8) and (11) for the whole range of 
C’, and compared with the experimental results. 
Since the intercept of equation (14) is so small, 
the calculated values for the more alkaline 
solutions are virtually independent of C*, and it 
is therefore appropriate to use a value of C* 
estimated for the three or four least alkaline 
solutions. The compositions and ionic strengths 


of the solutions are shown in Table 1. The major 


Table Sodium hydroxide 


Initial conen. (g mole B 
mole /em* 
Run ‘6 
OH | CO, SO,’ 10°S 10°B 
115 0-045 0-03 0-96 0-17 0-004 
116 0-085 0-07 0-91 0-26 0.006 
117 O15 0-13 0-78 0-48 0-02 
118 O21 0-19 0-70 0-68 0-04 
119 0-62 0-30 0-42 1-7 0-16 
120 1-16 0-38 0 0-40 
122 1-50 0-25 0 4-2 0-61 


contribution to the ionic strength of the less 
alkaline solutions is made by sodium sulphate, 
added to increase the viscosity : treating them 
as 1M solutions of sodium sulphate, of ionic 
strength 3M, one finds C* 1-7 x 10° g mole ‘ml. 
Using the method of van KreveLen and 
Horriszer [13] to predict the solubility in 1M 
sodium carbonate one finds the same value. 
whereas if one uses the figures of Nusina and 
KRAMERS [12] one finds C* = 1-8 x 10°, 

These figures suggest that the substitution of 
CO,” for SO,” will have little effect on the 
solubility. We shall adopt here the value 1-7 x 
10° gmole/ml. (It may be noted that during 
absorption of carbon dioxide the OH’ at the 
surface is rapidly converted to CO,”; the 
resultant increase in the ionic strength is, how- 
ever, small in the less alkaline solutions). 

Taking D = 1-2 x 10° em? /see equations (7), 
(8) and (11) are used to calculate S as a function 
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of C’. In Fig. 8 the calculated curve showing 

10° (s 2c* [P D as a function of 10°C’ 
D’ 

is compared with the experimental results. The 

agreement is good, but it happens that over the 

range of experimental conditions the closeness of 

the agreement is not very sensitive to the valu 


assigned to C* even in the less alkaline solutions. 


Sena 


~0°046),g mole/cr? sec’? 


Fic. 8. Comparison of experimental and caleulated 
values of S. 


The agreement therefore gives little informa- 
tion about the accuracy of the estimate of C*. 
However, it may be noted that th experimental 
value of S for 1-05 M sodium sulphate is 0.063 
10° compared to the caleulated value of 
0-066 x 10° g mole cm* sec! for infinitely dilute 
sodium hydroxide. Since C* in the former solution 
is very nearly equal to that estimated for the 
weaker sodium hydroxide solutions, this indicates 
that the latter estimate was correct. 

The rate of absorption of carbon dioxide in- 
to sodium hydroxide solutions is not proportional 
to C* as it is in the case of buff r solutions and 
non-reacting liquids. Strictly speaking, therefore, 
both the solubility of carbon dioxide and the 
observed rates of gas absorption should be 
corrected for the diluting effect of the water- 
vapour present in the carbon dioxide. Since the 
vapour-pressure of water at 25 C is 23-5 mm. the 
true absorption rate of carbon dioxide is 3 per 
cent less than the observed rate of passage of 
moist carbon dioxide, and the solubility is 
3 per cent less than that of dry carbon dioxide 
at atmospheric pressure. In view of th experi- 
mental errors and other uncertainties. these 


corrections are small enough to be ignored. 
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particularly since they tend to counteract one 
another. 

The rate of absorption of carbon dioxide into 
the stronger sodium hydroxide solutions is so 
fast that the temperature at the surface may be 
suspected of rising appreciably and invalidating 
arguments based on an isothermal system. This 
temperature rise may be calculated approximately 
by assuming that Q — 8 \/t and that the heat of 
reaction, as well as the heat of solution, is 
liberated at the geometrical surface of the liquid, 
Both of these assumptions exaggerate the rise. 
If H is the sum of the heats, then the rise in 


temperature in time ¢ will be [7].: 


aT SH do (15) 
TAI VOL 0) 2pa x 
where p= density, o specific heat, and 
x thermal diffusivity of the liquid. Putting 
eal Cem’, x 1-46 « 10% see, 
Hi 6,260 cal/g mole [7], we find for 15M 
sodium hydroxide where S$ +2 10° mole 


em* sec*?, A T = 0-61°C. The rise in temperature 
is therefore not large enough to exert any sub- 
stantial effects, 

From the foregoing it appears that after a 
certain exposure time the limiting rate of absorp- 
tion into sodium hydroxide solutions is con- 
trolled by the diffusion of CO, and OH’ to the 
reaction site, and that the effective diffusivity of 
the latter is about twice that of CO,. However, 
the actual situation is more complicated than the 
model on which the “ instantaneous reaction ” 
equations are based. In the first place, OH’ ions 
cannot diffuse independently of the other ions 
present, as the requirements of electrical 


neutrality must be satisfied. 


Fic. 9. Concentration distribution during absorption of 
carbon dioxide in sodium hydroxide solution. 
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Fig. 9 shows the way in which the various 
solute species might be momentarily distributed 
in the neighbourhood of the surface. In every 
element of solution the concentration of Na* must 
be equal to the combined concentration (in 
equivalents) of CO,", OH’ and HCO,’. The fact 
that the mobility of OH’ is greater than that of 
CO,"’ means that OH’ can diffuse to the reaction 
zone more rapidly (equivalent for equivalent) 
than CO,” can diffuse away from it. 

As a result, the concentration of negative and 
hence of positive ions is greater in the reaction 
zone than in the bulk of the solution. The 
resulting concentration gradient of Na* ions is 
opposed by an electrical potential gradient, the 
diffusion of Na* ions into the bulk of the liquid 
being opposed by their bodily transport towards 
the surface. This potential gradient also opposes 
the transport of OH’ to the reaction zone (and 
expedites the transport of CO,"’ away from it). 

The net effect will be to reduce the effective 
diffusivity of OH’ to less than its “ self- 
diffusivity "’ in the sodium hydroxide solution. 
The value of the latter quantity may be estimated 
roughly as follows. The experimental evidence 
available concerning the self-diffusion of in- 
dividual ions [16] shows that the value, D,*, of 
the self-diffusivity at infinite dilution is given by : 

D,* = Uy, (16) 
where u, is the ionic mobility at infinite dilution, 
The value of u, for OH’ ions at 25°C is 2-06 « 10% 
em*/V sec [17], whence D,* = 5-3 


The effect of finite concentration can be partly 


10 em* /sec, 


accounted for by changes in viscosity ; thus the 
quantity (D* «/p,) varies less than D* itself. 
Other effects, due to interionic attraction, cannot 
be accurately allowed for in concentrated solu- 
tions. The measurements of the self-diffusivities 
of halide ions in alkali halide solutions made by 
Mitts and Kennepy [16] show that values of 
(D* u/o) are some 10-22 per cent less in 2-3 M 
solutions than at infinite dilution. For the 
purposes of rough estimation it is reasonable to 
assume that it is the ionic strength of the solution 
which is of primary importance in determining 
the extent of this reduction. 
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The ionic strength in the 1-5M sodium 
hydroxide solution was about 2-3 M (owing to the 


presence of sodium carbonate), so the estimated 
value of (D* is 48 — 41 10° em? /see. 
Since p/p, = 1-57, the estimated value of D* 
is 3-1-2-6 

value 2-5 
* effective ” 


The method of estimating D* is so rough that this 


10° em?/sec, compared with the 
10° cm?/sec found for the 
diffusivity D’ in gas absorption. 


results lends only slender support to the state- 
ment that D’ < D*, However it seems likely 
that D* can be predicted in this way (when the 
requisite data are available) and used in place of 
D’ with a precision suflicient for many purposes, 
particularly as the square root of D’ is the factor 
determining the rate of absorption. 

A further complication is the reaction of 
dissolved CO, with CO,” to form HCO,’ in the 
region between the reaction zone and the surface. 
The reaction which actually occurs is reaction (1) 
between CO, and the OH’ ions which are present 
in a solution of CO,” by virtue of the equilibrium : 


CO,” + H,O =>. HCO,’ + OH’. (17) 


If the non-ideal behaviour of the solutes is 


ignored it is easily shown that : 


K,,. 


OH’ 18 
[OH'] = [HCO,’ (18) 
where kK, is the second dissociation constant of 


H.CO,. Thus the local rate of reaction is given by : 


[HCO,'] 2k’ [OH’) [CO,] 


_ [CO,] [CO,”] 


K,  [HCO,'} 


The method of integrating this equation to give 
the total amount of HCO,’ formed in a given 
time is detailed in Appendix 2. The fraction of 
the CO, dissolving in 0-20 see which reacts to 
form HCO,’ is greatest in the most dilute solution, 
namely 0-045 M, in which it is less than 7-5 
per cent. 

To sum up the discussion of the experiments 
with sodium hydroxide, although the results are 
readily explicable, quantitative prediction would 


be difficult in many circumstances because of the 


complexity of the processes involved. The 
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initial homogeneous reaction to produce CO," 
is followed by depletion of OH’ at the surface, 
leading in time to a diffusion-controlled reaction : 
superimposed on this is the secondary reaction 
leading to the formation of HCO,’. These observa- 
tions should be taken as a caution against over- 
simplification when predicting or explaining the 
results of gas absorption under industrial condi- 
tions. 

Absorption into buffer solutions. Sodium 
carbonate- bicarbonate buffer solutions were used; 
The pH 


was in all cases greater than 10, and in these 


their compositions are given in Table 2. 


circumstances [18] dissolved carbon dioxide is 
Prov ided that the 
reaction does not deplete the buffer sulliciently 


removed by reaction (1). 
to alter the concentration of OH’ appreciably, 
the rate of reaction, R, of carbon dioxide per unit 
volume of solution will be given by : 


R = k’ (CO,} [OH’} = k [CO,}, (20) 


with & constant for a given buffer solution. The 
value of k’ is likely to vary somewhat from one 
solution to another, even though the ionic 
strengths of the solutions were all about the same 
(2-3-2-5 M), since specific salt effects can be 
expected to occur at these concentrations. In 
addition the concentration of OH’ is not lik ly 
to be strictly proportional to the buffer ratio 
[CO,”)},/[HCO,"} [as implied by equation (18) ]. We 
should not therefore expect the experin ntally 

observed value of k to be necessarily proportional 
to the buffer recently 


published work of Nusinc and Krawers [12] 


ratio: however, the 


suggests that departure from proportionality is in 


Table 2. 
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fact small, perhaps owing to a cancellation of 
effects. 

Neglecting for the moment the depletion of the 
buffer by the reaction, the relation between 
@ and ¢ is given by equation (4). If erf ykt and 


e™ are expanded in powers of kt this becomes : 


l 


x + kt (Kt)? 


30 210 (82) 

It transpires that the maximum value of kt 
encountered in these experiments was less than 
1.5 so that the fourth term contributes less than 
2 per cent of the total and may be ignored in view 
of the experimental errors. The values of C* VD 
and of k may be found in principle by fitting the 
first three terms of equation (21) to the experi- 
mental results by the method of least squares, 
The first term corresponds to absorption without 
reaction, and the greater the value of At — that is, 
the more alkaline the solution and the longer the 
time of exposure - the greater is the influence of 
the reaction and the accuracy with which k can 
be determined. 

However, the range of kt within which equation 
(21) can be expected to apply is limited by 
depletion of the buffer solution. It is shown in 
Appendix 1 that for the most and least alkaline 
solutions used (Runs 112 and 107) the values of 
kt leading to a change of 10 per cent in the buffer 
ratio at the surface are about 0-85 and 3-2 res- 
pectively. If a preliminary estimate can be made 
of the value of & it is possible to estimate the time 


of exposure above which depletion exceeds, say, 


Buffer solutions 


Conen. (g mole /1.) lonic 


Run — strength 


CO,” HCO, I (g ion 1.) 


Baffe r ratio 


C* /D 


0-49 2-29 
108 0-62 0-535 2-40 
109 0-68 Ow 244 
110 O-735 295 250 
111 0-765 0-225 2-32 
112 O-175 2-52 


[COs (g mole k (see) k 
[HCO’s| em? sec?) (see *) 

0-60 67 O-70 1-16 

116 6-8 1-51 | 1-30 

1-70 6-8 1-45 O-B8S 

2-47 67 2-25 ool 

340 6-8 2°55 O75 

445 746 374 Ost 


= 195 
| all 
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10 per cent and which may be arbitrarily taken 
as the limit above which equation (21) ceases to 
be applicable. 

There are various ways in which k can be 
estimated, starting from the observation of 
Nusinc and Kramers [12] that its value is 
closely proportional to the buffer ratio and in- 
dependent of the ionic strength, over a range of 
conditions not very different from those of the 
present work. Thus they report one value of k 
at 25°C, from which its values for the present 
Runs 112 and 107 may be estimated as 3-5 and 
0-5 sec’ respectively. Taking their value for 
k [ HCO,’ | at 20°C, namely 0.55 sec 
and employing the temperature coefficient for 
the reaction determined by Pinsent, Pearson 
and RouGuron [11], we find the values for Runs 
112 and 107 at 25°C again to be 3-5 and 0.5. 

On the other hand, if we assume that in the 
present work the degree of depletion in the runs 
with the lowest alkalinity (107 and 108) was small 
enough to be ignored (as will subsequently be 
shown) we can use the value of k in these solutions 
to predict the value of & in the more alkaline 
solution. The results for the least alkaline 
solutions cannot be expected to give the value of 
k with great accuracy, because of the small 
contribution made to @ by the reaction. 

It will be seen that the value of k [ HCO,’ | 
[COs } for the less alkaline solutions is con- 
siderably greater than the value derived from the 
work of Nigsinc and Kramers, which was about 
0-8; however, in case there was a systematic 
discrepancy between their results and ours it was 
considered appropriate here to extrapolate from 
Run 108, which gave the largest value of k 
[HCO,'| namely 1-30 sec’, giving the 
value of 5-7 sec"' for k in the most alkaline 
solution (Run 112), and 0-8 sec"' for the least 
alkaline (Run 107). For Run 107 the maximum 
time of exposure was 0-263 sec, giving a maximum 
value of 0-21 for At; this is so much less than the 
value, 3-2, at which the surface depletion reaches 
10 per cent that we can take it that equation (21) 
applies to the entire course of the run. 

For Run 112, with a maximum exposure time 
of 0-215 sec, on the other hand, the estimated 
maximum value of kt is 1-2, which is considerably 
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greater than the value of 0-85 for 10 per cent 
depletion ; only the observations below 0-15 sec 
would be expected to conform to equation (21) 
with reasonable accuracy if the estimate of the 
value of k were correct. Similar calculations show 
that the change in the buffer ratio was less than 
10 per cent in all runs except 112. 

When the first three terms of equation (21) 
were fitted to the experimental values of Run 112 
for times less than 0-15 sec, k was found in fact to 
be about 3-7 sec™'. It was concluded that, since 
the value of kt for the maximum exposure time 
of 0-215 sec was therefore in fact below 0-85, 
the whole range of experimental values should 
conform to equation (21). The values of C* , D 
and k shown in Table 2 were obtained from the 


whole range of Run 112, and of the other runs. 


% 
+ + + 
° 
1 + + + - 4 
1 2 3 5 
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Fic. 10. Absorption into buffer solutions (k v. buffer ratio). 


In Fig. 10 the values of k& are plotted against 
the buffer ratio. The results are compatible with 
the relationship : 

(COs) 
[| HCO,’ 


k= 08 (22 


] sec 


This is in agreement with the single result of 
NiJstInc and Kramers [12] at 25°C, and also 
with their results at 20°C if corrected by the 
temperature coeflicient of PiInsenr, Pearson and 
RoOUGHTON {11}. 

No direct comparison is possible between the 
values of C* yD obtained by Nuysinc and 
Kramers and by ourselves, since the ionic 
strengths and viscosities of the solutions were 
different in the two cases. Taking D in our case 
10° em? /sec, we find 


to have been 1-2 


VOL. 


C* = 2-0 
runs 107-111, in which the ionic strengths 
varied from 2-29—2-52 M. 

Using the method of van Kreve.en and 
Horrizer [13], in conjunction with the coefti- 
cients found by Nigsinc and Kramers [12], the 
predicted solubility in these solutions varies from 
2-13-2-02 


good, and provides some support for the validity 


10” gmole ml. This agreement is 


of the method of estimating solubilities and for 
the conclusion of Nusinc and Kramens that the 
solubility in  carbonate—bicarbonate buffers 
depends on the ionic strength and nature of the 
cation, but not on the buffer ratio. However the 
present findings represent a relatively small 
extrapolation from their experimental results. 
The much higher apparent value of C* in Run 112 
suggests that the relationship may not hold at 


high buffer ratios. 


CONCLUSIONS 


The results obtained for the absorption of 
carbon dioxide in water are consistently some- 
what lower than the values calculated from the 
most reliable values of C* and D. In view of the 
work of CuLLEN and (2, 3}, which 
showed no such discrepancy, this is thought to be 
due to systematic experimental error rather than 
to a low accommodation coeflicient, as previously 
suggested. 

The results for absorption into inert salt 
solutions yield values of D consistent with the 
assumption that D varies invers« ly with the 
viscosity of the solution. 

The early stages of absorption into sodium 
hydroxide solution are presumably accompanied 
by a pseudo first-order reaction between CO, 
and a virtually uniform concentration of OH’, 
However, calculation shows that depletion of 
OH’ near the surface will become appreciable in 
a time much less than the shortest exposure time 
observed in these experiments; this cd stroyvs 
the first order character of the reaction and makes 
it impossible to deduce the reaction velocity 
constant from the observations. 

At higher exposure times the rate of absorption 
is governed by the diffusion of OH’ and CO, to 


the reaction site ; the results agree quantitatively 
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10° gmole/ml as the average of 


with predictions, and the effective diffusivity of 
OH’ is found to be about twice that of C¢ »,. The 
effective diffusivity is expected to be less than 
the self-diffusivity of OH’ in the same solution ; 
a rough estimate of the latter quantity, based on 
the ionic mobility, supports this conclusion. 

At still higher exposure times, reaction of the 
CO, with the CO," 


give HCO,’ will become important and will 


formed near the surface to 


complicate the calculation of absorption rates, 
In these experiments it is calculated that less 
than 7) per cent of the C¢ ), absorbed was con- 
verted to HCO’, in the worst case. 

The absorption of CO, into CO,” HCO,’ 
buffer solutions is accompanied by a first order 
reaction so long as the buffer is not appreciably 
depleted near the surface. Calculation shows 
that the depletion did not become substantial 
during these experiments. The calculated values 
of the velocity constant were compatible with 
those found by other workers. 
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APPENDIX I 
Tue Derterion or Sopium Hyproxipe axp 


Burrer Souerions spy Reaction 


Sodium hydroxide solutions, Let the initial uniform 
concentration of OL be C(O), and its conce ntration at 
a distance a from the surface at time t be C (x. t). The 
concentration of CO, at the surface is C* at all times : 
it diffuses into the solution, in which it undergoes a second 
order reaction with OH’. During the initial stages of the 
process, when the depletion of OH’ is small, the concentra- 
tion of OH’ can be regarded as uniform and the reaction 
is first order in the local concentration of ( Os, the velocity 
constant & k’ C(O) having the same value everywhere. 
The method to be described makes use of this approxima 
tion, and thus only allows small degrees of depletion to be 
calculated. 

Putting 

V =) 0) 23) 
Co) 
the simultaneous diffusion and reaction of OH’ are ex 
pressed by the equation : 
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where D’ is the diffusivity of OH’ and ¥% is the concentra- 
tion of unreacted CO, at (x, t) and is determined by the 
equation for diffusion accompanied by a first order 
reaction, the solution of which has been given elsewhere 
{10}. If the expression for X is inserted into equation (24) 
it becomes mathematically analogous to that for the 
changes of temperature produced by a first order reaction 
accompanying gas absorption; this problem has been 
considered elsewhere {[7| and a comparison with the 
present problem shows that if V (O, ¢) is the value of 
Vat the surface at time ¢, and L (0, its Laplace 
transform, then : 


|» 
C(O) ND 


Liv 
2kc* D pip(b D) + kD’! 


Here p is the parameter of the Laplace transformation 
and D the diffusivity of dissolved CO,. In order to simplify 
the expression, it will be assumed that the diffusivities of 
CO, and OH’ are equal, so that D D. (This will 
overestimate the extent of depletion, since in fact D D 
as shown in the body of the paper ; the magnitude of the 
error so introduced is discussed later), Mquation (25) 


then becomes : 


Liv (0,0)} Vl (26) 


The inverse transform of the first term on the right has 


been found elsewhere [7]: using this result : 


C(O) — C (0,8) (;) (; 
. (27) 
2c¢ n! n! 
= (kt). 


The function f, (At) can be read from Fig. 1 of | 
by subtracting 1 from the ordinate ; however, it can be 
approximated with suflicient accuracy for the present 
purpose by ; 


fy (ht) kt/2, kt <1; (28) 


In order to calculate the time required for 10 per cent 
depletion at the surface, we put : 
C(O) — 2C 


| (At). (0) 
Coy) Coy) 


Taking C* = 1-7 « 10°° g mole/ml, we find for 1-5 sodium 
hydroxide, f, (At) t4. Using the approximation given 
in (29) we find At 21. If we take the value of & to be 
about 107C (O) sec”', t is found to be about 1-5 x 10° see. 
A similar calculation shows that for 10 per cent depletion 
of 0-045 M sodium hydroxide, /f, (At) (22, and using 


approximation (28) ¢ is found to be about 6-7 x 10 sec. 
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If D’ > D the inverse transform of equation (25) takes 
the form : 


Cc y 
(0) : V(O,) (ety, (31) 
2c* D 


= fo (kt). 


The coefficients »,, and a graph of fy (At) are given in 
[7]; fo (kt) tends to 2 \/kt/m for large values of At, 
and to At for small values. If the ratio D’/D were equal 
to 10, calculation shows that the time required for 10 per 
cent depletion at the surface would be about 9 x 107% sec. 
for 1-5M sodium hydroxide and 1-7 x 1078 see for 0-045 M 
sodium hydroxide. Actually since the effective value of 
D’/D is about 2, as shown in the body of the paper, the 
times required will be considerably less than this. 

Buffer solutions. As before, it is assumed that the change 
in composition of the solution is small enough for the OH 
concentration to be treated as uniform, so that the CO, 
undergoes a pseudo first-order reaction with a velocity 
constant k k{OW). If the initial concentration of 
CO, is C(O) and its concentration at a, tis Cy(v, 0), and : 

C, (4,0 


(32) 
' (0) 


then 


where D, is the diffusivity of CO,"’.. Assuming the diffusivi- 
ties of CO, and of CO,” to be equal, comparison with [7] 


shows the solution to be : 


C, (0,8) c* 
C,(0) ©, (0) 


Sy (ht), (34) 


where f, (Kt) is given by equation (25). If the initial 
concentration of HCO,’ is Cy (OQ), and its value at wv, f is 


(a, the equation expressing the change in concentra- 
tion of HCO,’ is: 


where 


(4, 


A a (36) 
(O) 


and Dy, is the diffusivity of HCO,’. If this is equal to the 
diffusivity of CO, : 
Cy (0, t 2c* 
C2 (0) (0) 
The fractional surface depletion, A, of the initial buffer 
ratio C, (O)/C, (O) is thus : 
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da 
0,z>0,t 0 : 
= dD. 54 (35) 
a C2 ( = O,a 
| 


(O) [Cg (O) + 20% f, (ht) 

Considering the most and least alkaline buffer solutions 

used, and putting C* = 20 « 10°5 g mole /mil, it is found 

that the values of kt required for A to reach a value of 

O-1 are 0-85 and 3-2 respectively. 


APPENDIX 2 
Tue Formation or Bicarponate Sopi™ 
Hyproxipe Souutrions 


In order to calculate the amount of ¢ Og which undergoes 
the secondary reaction with CO,” to form He O,', the 
actual situation is simplified by assuming “ instantaneous ™ 
give CO, 


primary reaction between CO, and OH’ to 
If the secondary reaction is left out of account, there will 
be a layer of depth 2 8 \ ¢ [8 being defined by equation 
(8)], in which all the OH’ has been converted to COs, 
and in which the concentration ¥ of unreacted COg is 
given by: 


erfe Dt) — erfe (B/y D) 
ce erf (8 D) 


x being the depth below the surface [14]. Thus at a 
depth the concentration of Og is zero at t = («7/2 
and thereafter rises progressiy ely to C*, The concentration 
of CO,” rises from C” for t < BP to (¢ 4c") 
for t> (7/2 BY, C’ being the initial concentration of OW 
and C” the initial concentration of CO,” present in the 
solution. In the present calculation the diffusion of both 
CO," and HCO,’ ions is neglected, the changes in their 
concentrations at a given depth being regarded as being 
caused sole y by the local secondary reaction 

The errors introduced by these two approximations are 
opposite in sign; that due to neglecting the diffusion of 
CO,” is probably the greater, and will lead to an over- 
estimate of the amount of HCO,’ formed. The depletion 
of the CO, concentration by the secondary reaction is 
ignored, and its concentration is assumed to be given 
by equation (39); this also gives an overestimate of the 
HCO, formed, as does the assumption of an “ instan- 
taneous’ primary reaction. The calculated extent of 
reaction will, however, be of the right order of magnitude. 

If in equation (19) we put a, (CO, r. 
§ + C’— a), the equation can be integrated 
with respect to ¢ for constant x : 


a 
ada 
ac’ a) 
t 
= 4 AK, erfe (x /2 Dt erfe (B/y D) dt. 
Ky erf (8 
(2/28) 


The concentration dy is the equilibrium concentration of 
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[HCO ’) in a solution of sodium carbonate alone in concen- 
tration (C 4 C’); ag was negligibly small in all cases, 
and may be ignored. Equation (40) then becomes on 
integration : 


k’ 
(2C’ c | = 
2C ( a Ky erf (B/ 


| 


t 
a e Di. (41) 
2, \ 


If we take ¢, the exposure time, to be 0-20 sec, a can be 
evaluated as a function of «. and a graphical integration 
from a Otor — 28 \ tygives the total amount of HCO,’ 
formed, 


For the purposes of the calculation we take 


Ce = 1-7 g mole dD 12 105 em* /see, 
A 107 ml/g mole sec. The total amount of HC Os 


formed is about 7-5 per cent of the total amount of ¢ Os 
absorbed in the case of 0-045 M sodium hydroxide, and less 
for more concentrated solutions. 


NOTATION 
a = concentration of HCO,’ 
A defined by equation (10) 
B defined by equation (11) 
Cc (with various subscripts and superscripts) 
concentration of substances other than ¢ O, 
ce saturated concentration of ¢ O, in water or 
solution at 25°C, partial pressure 760 mm. 
D = diffusivity of ¢ Oy in water 
D’, D,, Dy = diffusivity of solutes other than CO,g in 
water 
D* = self-diffusivity of OH’ ions in sodium 
hydroxide solutions 
Faraday’s number 
il sum of heat of solution and heat of reaction 
of CO in sodium hydroxides solution 


4 pseudo first-order reaction velocity constant 
4 second order reaction veloc ity constant 
kK, dissociation constant of water 
kK. second dissociation constant of H,COs 


1 = depth of liquid film; length of exposed 
periphery of drum 
nN no. of moles of solute reacting with one 
mole of CO, 
p — vapour-pressure of water; parameter of 
Laplace transformation 
P — total pressure 
q rate of absorption of ¢ O, 
Qit} = amount of CO, absorbed per unit surface 
area in time ¢ 
R = gas constant; rate of reaction per unit 
volume 
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defined by equation (11) «x = defined by equation (6); orthermal diffusivity 
time of exposure of liquid 
temperature i defined by equation (8) 


AT = temperature rise at surface = fractional change in buffer ratio 
te = ionic mobility of OL ions at infinite dilution = viscosity of liquid 


J velocity of drum surface. viscosity of water 
0, 1), U2 = quantities defined by equations (23), (32), = coefficient in equation (31) 
(36) density of liquid 
width of drum specific heat of liquid 
depth below liquid surface concentration of CO, in liquid 
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Equilibrium theory of water treatment 


K. Kiamer and D. W. van Kreveten 


Central Laboratory, Staatsmijnen, Geleen, Holland. 
( Received 12 October 1957) 


Abstract— It is possible to determine in a simple way — assuming the occurrence of an instan 
taneous establishment of equilibrium — such data as the break-through capacity and the relation 
between capacity and leakage of ion exchangers. A theoretical derivation is given of a relation 
determined experimentally by Wuearon and BauMAN, 

Application of the measured equilibria to the purification of water leads to a first approxi- 
mation of the capacity to be expected of, e.g., cation exchangers, if several cations are to be 
removed and a buffering ion is present. 


Résumé—Supposé qu'il y ait un établissement instantaneé dequilibre, il est possible de déterminer 
de fagon simple, des données telles que la capacité de traversée et la relation entre la capaciteé et 
la fuite des échangeurs d’ions. On donne la dérivation théorique d'une relation déterminée 
experimentalement par Wuearon et BAUMAN. 

L’application des équilibres mesurés & la purification de l'eau méne a une premiere approx 
imation de la capacité a attendre, par exemple, d'un échangeur de cations. si plusicurs cations 
doivent étre enlevés et qu'il y ait un ion-tampon. 


Zusammenfassung —Es ist méglich, auf einem einfachen Wege — unter der Annahme, dass sich 
das Gleichgewicht sofort ecinstellt — Daten wie die Durehbruchkapazitét und die Bezic hung 
zwischen Kapazitiét und Schlupf von lonenaustauschern zu bestimmen. Kine von Wheaton und 
Bauman experimentell gefundene Beziehung wird theoretisch abyeleitet. 

Die Anwendung der gemessenen Gleichgewichte auf die Wasserreiniguny fiihrt zu einer ersten 
Anniherung der zu erwartenden Kapazitét von z.B. Kationenaustauschern, wenn mehrere 


Kationen zu beseitigen und ein Pufferion anwesend ist. 


Part II [2] described the general method to may be used. The important point here is the 
be followed for determining equilibria on ion shape of the isotherm. 


exchangers; some of the data obtained were 
1.1. Break-through capacity without leakage ; 


[D, (0)] 


used for computing capacities of columns in We shall treat in suecession the various shapes 


presented and interpreted. The present part deals 


with the way in which these equilibria can be 


which instantaneous establishment of equilibrium — the equilibrium curve may assume, 
oceurs, and with the application of these equilibria 

in the ification of water. 

n purification of wa f, (cx) 


Break-Turoven Capaciry 


An important practical quantity of any ion 


exchanger is the capacity as a function of the 


permissible leakage. By way of example we shall O° 
lake the Na* H* equilibria, with Cl- as the atta, 

negative ion, for cation exchangers. For a buffering Fic. 1(a). f,” (a) <0. Graphic determination of the 
anion or anion exchangers similar derivations break-through capacity. 
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(a) The equilibrium curve is linear (f,,'° («) 0) 
or concave with respect to the x axis, f,'(a)< 0 
(Fig. la). The front always remains sharp and 
the break-through capacity is equal to the total 
capacity. 


= f, (1) 


(b) The equilibrium curve is convex with respect 
0 (Fig. 1b). 
that if a solution of hydrochloric acid and common 


to the « axis f,"’ (x) It is now clear 
salt is sent through a column containing a H 
ion-loaded exchanger, a diffuse front will be 


formed, The result is that after a short time Na 


| 


Fic. 1(b). f, (x) > 0. the 


break-through capacity. 


Graphic determination of 


ions will be contained in the effluent. The number 
of Na 


moment is called the break-through capacity. 


ions taken up by the exchanger until that 


The height of this break-through capacity can be 
determined as follows. Our starting point is the 
equation for the diffuse front : 


ai 


Sf, (x) 


ay 


from which follows, after integration : 
avt 


fy 3 
ay - Sf, (x) 


The constant d(x), which in- 


dependent of t, gives the distribution of Na 


integration 
ions 
over the column before the experiment was 
(a) for t 0). 
column to be loaded only with H 


started (a As we assumed the 


ions, (x) 


equals zero in the case in hand. If the length of 


the column be / and the break-through time tp, 
then after tp, the point of the diffuse front « = 0 
(0). 


has arrived at 1, where f, (x) From 


equation (3) it follows that : 


avlp 
ay + Sf, (0) 


/ 


(4) 


If we call the equivalent fraction of the influent 
%y, then xg avty represents the number of Na 
that the the 
break-through the Na 
ions contained in the free space is smaller than 
that the 
greater than a, lSf, (0). Ina first approximation 


ions have entered column until 


moment ; number of 


contained in exchanger is 


xq la y is negligible in comparison with x,» LSf, (0). 


Strictly speaking, x US/,’ (0) represents the load of the 


exchanger if the flow of liquid is stopped after 


hy 
and the column is then drained or the 


='tp 
solution is replaced by distilled water; in this case the 
last of the solution has just left the column after tp (this 
liquid was still completely free of Na* ions) and, conse- 


quently, all the Na* ions that have entered the column. 


tq tg Wty — ag lay = 2g ISS,’ (0) 


are present on the exchanger only. 


Consequently, the break-through capacity of 
the column is x, lSf, (0), and, as the volume of 
the exchanger in the column is lS, the break- 
through capacity per volume of exchanger is : 


(O) 


(c) The equilibrium curve contains an inflexion 
point (Fig. 1c). If we 
in a preceding paper [1] about this type of 


tof, (0) (5) 


combine what was said 
isotherm with the treatment given in 1.1 (b) in 
the present section, we see that for | eg (a)< O 
the break-through capacity is equal to that for 
f, >09, tof, (0). their 
origins both isotherms have the same shape. 


VIZ.. since near 


f, (cx) 
i 


| 


(x) > 0. Graphic determination of the 
break-through capacity. 


Fic. l(c). fy 


‘ 
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that for 


capacity is 


If, however, f,"’ («) > 0, 
0<a,< a, the break-through 


JS, (%o), and for 4, < ag< 1: 


| D, (0) toSe (x,) 


1.2. Dependence of the capacity on the permissible 
leakage ; (D, («) ) 


Here also, the cases treated under 1.1 may be 


we sce 


(6) 


distinguished. 
(x) 0) 
0). 


We are here concerned with a sharp front; the 


(a) The equilibrium curve is linear (f, 
or concave with respect to the « axis (f, (a) < 
break-through capacity is equal to the total 
capacity [equation (1))}. 

(b) The equilibrium curve is convex with respect 
to the x axis; (f.' (a) 0). At the start the 
column contains H”™ ions only. The loading with 
Na 


ion concentration in the 


ions is continued to the point where the Na 

effluent has reached a 
this point ts reached at the 
At this moment the 
number of Na* ions that have entered the column 


certain value «a; 


moment ¢, after the start. 
is : 
(a) ) 


Of the right term of this equation it is only the 
second part that is of interest (see Section 1.1). 


% avl, = agl(ay 


Consequently : 

(7) 
The number of Na* ions issuing from the column 
in the same period is : 


ty avt, = a, ISf, (a) 


avdt, = « avt, 
o 
which after elimination of t,’ from equations (7) 
and (8) results in : 


IS (x) —f, (9) 


so that : 


IS {ao fe (x) —af, (a) + Sa (a)} 


is retained on the exchanger [cf. equations (7) 
and (9)]. Consequently, the number of Na” ions 
taken up per unit volume of exchanger at a 
given leakage a may be expressed as (cf. Fig. 2): 


W. vAN KREVELEN 


D, = (xq — a) fa (a) + fy (x) (10) 


a 


Graphic determination of the capacity as a 
function of leakage. 


It is further possible to calculate the fraction of the 


exchanger loaded with Nat ions at the time t,; this is 


represented by (ef. equation (10)): 
Number of Na* ions present on exchanger 


Max. number of Na* ions the exchanger is : 
capable of taking up 


(x) 4 Sq 
(*o) 


(11) 
For the loading efficiency [¢f. (10) and (7)] we find : 


Number of Na* ions present on exchanger 


Total number of Na* ions that have entered the column 


(a9 (a) + Sq (a) 


(12) 


From equations (11) and (12) it follows that the maximum 
efficiency of the column is : 


loading efficiency 
fraction loaded 
Sa (xq) 


Sa 


Maximum efficiency 


This relation has already been found experimentally by 
Wueaton and Bauman [3] in the regeneration of anion 
exchangers, which, of course, is a completely similar 
process ; however, the theoretical derivation has not been 
given before. 


(c) The equilibrium curve contains an inflexion 
point. With these equilibrium curves there is a 
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series of possibilities, which will be discussed 
below. ‘It is easily proved that one of the equations 
(1), (6) or (10) always applies. If f, | (x) > 0, 
there are two possibilities. 

In the first place: 0 < < the front 
then has the same sharp form as for f, (x) 0 
and the break-through capacity is equal to the 
fa (%o) [ Fig. la, equation (1)]. 
Le 1. With 


equivalent fractions « in the effluent 0 x x, 


total capacity 
In the second place: 4, 
(so in the sharp part of the front) the break- 
through capacity is a, f, (x,) [ Fig. le, equation 
(6)| and for a, < « z, (in the diffuse part of 
the front) equation (10) (Fig. 2) applies. 
If f, («)< 0 we may likewise distinguish 


between two possible cases. In the first place we 


a 


may have 0 < a, x. When the break-through 
capacity in the diffuse front is given by equation 
(10) (Fig. 2). 

In the second place: a, < a < 1. Here also 


equation (10) applies if 0 < « to (%» has the 


same meaning as ¢,. in Fig. 7, [1]), or, in other 
words, if « lies in the diffuse part of the front. 
If, however, < then « lies in the 
sharp part of the break-through front and 
(x) = f,(%) [equation (1), Fig. 1a}. 

The practical conclusion from this section is, 
therefore, that if the equilibrium curve and a, are 
known, the capacity with and without permissible 
leakage can immediately be graphically derived 
(ef. Figs. 1 and 2 respectively). The formulae 
derived in this section are summarized in the 
Fig. 3. 


2. APPLICATION OF THE MEASURED 
EQUILIBRIA TO THI PURIFICATION OF 
WATER 


In this section we will examine how, (under 
equilibrium conditions), the break-through capa- 
city and the load condition for a given composition 
of the water can be determined from the measured 
isotherms. The rules to be observed will be set 
forth in brief with reference to an example taken 
from practice ; finally conclusions will be given 
with respect to the mixing of exchangers. 

Natural waters contain chiefly the following 
ions : the cations of calcium, magnesium, sodium 


and potassium, and as anions, bicarbonate, 
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chloride and sulphate; of the non-dissociated 
compounds we mention carbon dioxide and 
silicon dioxide. Before proceeding to discuss the 
break-through capacity as mentioned above we 
shall pay attention to the correction which the 
isotherms (of the cation exchanger) have to 
undergo if the water contains a buffering anion 
(e.g. bicarbonate). For the equilibria described 
in a preceding paper [2] are valid only if the 
counter ion is obtained from a strong electrolyte. 


2.1. Equilibria with buffering anions 

The HCO, -equilibria needed for the calcula- 
tions may be derived from Figs. 2-7 of the 
preceding paper (Part II) [2] by means of the 


following equations : 


(14) 
and 


(15) 


2 
f,( 
re 
Co: conc. O-OIN 
“in ——~ Ca"; conc. 0-003 N 
Bicarbonate (2) 
1500) 
/ 
/ 
/ 
/ 


~“e—Chioride (1) 


| 


Fic, 4. Equilibrium curves of Dusarit VKG with HCO, 
and CI” as counter ions. 


VOL 
8 
A Cy 
co 
2 il 
oe which are identical with equations (9) and (10) 
1000/4 
50C 
iz 


in Part Il [2]. In this case AK is the first dis- 
sociation constant of carbonic acid, For any pH 
(and hence c+) the corresponding x can be 
calculated by means of equation (14). Since the 
relation f, (x) vs. a, is known, it is then possible 
by using equation (15) to calculate also the 
relation f,, («) VS. a. 

These new isotherms are shown for Dusarit 
VKG and Dusarit S in Figs. 4 and 5. The dotted 
parts relate to estimated values; this is the 
region of low «, relatively unimportant since it is 
data such as break-through capacity, etc., which 
are most important. As the caleium equilibria 
are dependent on the total concentration in 
solution, one would be inclined to expect the 
same of the Ca curve corrected to apply to 
HCO,” ions. It was found, however, that this 


influence due to concentration is negligible. 


resp. | 
1500} OF 

in —»Ca**;conc. O‘OO3N 


est ail Bicarbonate (2) 
7 
/ 
/ 
t 


1-0 


04 


—- resp. a, 


Fic. 5. Equilibrium curves of Dusarit S with HCO,” and 


Cl as counter ions. 


Therefore, since the Na’ H* equilibrium is 
in general independent of the total ion concentra- 


tion in solution, in the following the index a to 


f(«) will be omitted for convenience. 


If the water contains salts of strong and weak 
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acids and only one kind of cation the second 
equilibrium curve will depend on the concentra- 
tion of these two salts. Instead of : 


Na CH Cx (16) 
a new equation for the electric neutrality in 
solution is obtained : 

Cx Coy (17) 


and the following equations for the concentration: 


NaHCO 


NaCl 


z+y (18) 


From the equations : 


and (17) and (18) it follows that : 


kK Cw 
K Cu K 


By means of equations (22) and (15) the equili- 
brium curve for, e.g., the Ca** ion can be calculated 
if both Cl’ and HCO, 
then happen — this was found with polyfunctional 


ions are present. It may 


exchangers — that an isotherm of the shape shown 
In this curve there is a point 


characteristic (a) 0, 


in Fig. 6 results. 


/ Ce - 


representing the 


a, 


Isotherm with two inflexion points ; 


sere 


(a2) <0. 


voL. 
"Na 
Cert € 
(20) 
CHX 
Na Cyt “Hx (21) 
m equiv! 
| 
| 
/ | 
fo (ad) 
4 | 
O-2 0-8 | 
o 1 
6 
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Application of the considerations given in a 
preceding paper (Part I) [1] leads, also in this 
case, to the change in the shape of the chromato- 
gram of this remarkable isotherm (cf. Fig. 10a). 
The correction for the bisulphate equilibrium : 
HSO, * SO,” + H* 
proves to be negligible for the concentration 
range we are here interested in. 

From the measured equilibrium curves shown 
in Figs. 2-10 in Part II [2] the break-through 
capacity under various conditions may be read. 
For Dusarit VKG this is : 

700 m.eq. perl. for NaCl (a,? = 1) 

1030 » » 0801 N CaCl, = 1) 
1100 » » 0-008 N CaCl, («,? = 1) 
1650 » » Ca (HCO,), 


a fw 0-75, see 


(pH 7, 
Fig. 4.) 


The chromatograms which are obtained when 


a H* exchanger as mentioned above is loaded with 


calcium bicarbonate, calcium chloride or sodium 
chloride, appear as shown in Figs. (7a), (b) and (c). 


a 


Fig. 7 (a). Ca (HCO Fig. 7 (b). Cath, 


Fic, 7 (c). Natl. 


representation of the chromatograms 
* and Na 


Diagrammatical 


for the Ca ions. 


2.2. Several substances (cations) are in solution 


If two substances A and B are in solution, 
while for both of them f’’ (c) < 0 and for each 
concentration ¢ the aflinity of A to the absorbent 
is greater than that of B, the chromatogram for 
the loading of an empty column will appear as 
shown in Fig. 8. Substance A is absorbed more 
strongly and will expel part of substance B. 


Fic. 8. Chromatogram of a mixture of arbitrary substances, 


For example if a solution containing sodium 
chloride and calcium chloride is led over a H* 
exchanger, a chromatogram like Fig. 8 will be 
obtained (A = Ca** and B = Na*). Now the 
question is how many Ca** and Na* ions are 
present side by side on the exchanger in the 
In principle this must be 
for separate 
case. It has been found, however, that with Ca‘ 


range 0 <@< 


determined experimentally each 
and Na* ions as normally occuring in water the 


Na* content of an exchanger such as Dusarit 


x 
x 
and Na* 


Fic. 9 Chromatogram of a mixture of Ca** 


ions, with Cl” as the counter ion. 
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VKG is negligibly small in the range O < 2 < @,. 
The experiment is carried out as described 
before [2] except that several cations are used. 
The chromatogram appears similar to that 
shown in Fig. 9. 

If bicarbonate is also contained in the water 
a much heavier loading of the exchanger with 
Ca** ions is possible. Two cases are illustrated 
in Figs. 10(a) and 10(b). 


HCO, 
f(a) 
: cv 
a, l 


Fic. 10 (a), 


f (cx) 


y 
Fig. 10 (b). 


Fie. 10 (a) and (b). Examples indicating the influence of 
the HCO, ion on the chromatogram. 


A third case may be imagined, viz. that due to 
the presence of HCO, ions also, more Na* ions 
are taken up; however since the calcium ion 
concentration is generally greater than the 
bicarbonate ion concentration this case will be 
infrequent. Whether the first or the second case 
is encountered depends on the ratio of the 
chloride and bicarbonate ion concentrations 


[ef. equation (22) }. 


Studies on ion-exchange—III 


223 


With the chloride and bicarbonate ion concen- 
trations normally found in water and with an 
exchanger such as Dusarit VKG the second case 
will be encountered. Then the break-through 
capacity can be calculated as follows : 


Assume the following concentrations in the 

, me (at? .. ‘ 
water : = ¢y,++, Na* = cy,+, HCO,” =tyco,- 
and Cl 4+ SO,” = c,,; further assume that 
Coat+ > > > Which is generally 
the case —then the load of calcium bicarbonate 
pertaining to the « (or pi/) of the raw water holds 


for the part 0 < 2< a@,, and the load of Na* 


pertaining to « = -“*- for the part 2, < < lL. 


This latter point will be clear if it is borne in 
mind that the liquid passing 2, contains ¢,, 
m. eq./l. of electrolyte, of which cy,+ m. eq./l. 
are Na* ions and (c,, — cy,+) m. eq./l. are H* ions, 
The number of unit volumes of water that can be 
treated per unit volume of exchanger then is : 


(l — a,) f(a = 


/ (293) 


From this equation 2, can be solved ; the break- 
through capacity is : 


(24) 

An example of such capacity calculations is 
included in Table 1. 

It will be appreciated that capacity calculations 
as described above only give an impression of the 
order of magnitude that can be reached under 
given conditions. This is partly due to the fact 
that equilibrium conditions do not always exist 
everywhere in the column. As a result of the 
kinetic effects the sharp fronts are spread and 
give lower capacity values than shown above. 
Conclusions of a comparative nature can always 
be drawn, however ; in the following an example 


will be given. 


2.3. Mixing of ion exchangers of the same kind 


By mixing two or more ion exchangers of the 


same kind, either cation or anion exchangers, a 
mixture may be obtained which behaves as a 


OL . 
> 
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Table 1 An example of the results that can be obtained 
with such mixtures is shown in Table 2 


Exchanger DusaritS  Dusarit VKG 

The mixing of exchangers of the same kind should be 

Coa 6-9 8-5 6-9 8-5 distinguished from the mixing of cation and anion 
exchangers in a so-called mixed bed The advantage of 

Na 0-6 O73 | O«6 OTS this modern technique is that the acid formed by a 
cation exchanger from a salt solution is immediately 

Cnat/Con 1 0-4 , 0-4 caught by the anion exchanger. As a result, the equili- 
brium will always lie in an extremely favourable position. 

ti a (HCO), (PH 7) 1000 1000 «16501630 Regeneration of such a bed is possible only after separation 


of the cation from the anion ex: hanger. 


Table 2 


ail 085 O48 O83 


(x) for the exchangers 


Break-through capacity : 


Calculated 20 5801400 970 


30°. Imac C 12 Imac C 12 
Dusarit Stamion Stamion VO 
As measured at / 0-88 m VAG Calculated Found 5 


and J 44 m/hr 670 1000 S80 


18S 200 
new ion exchanger. When the equilibria of the 0-20 743 720 10 


constituent exchangers (which can. in principle, 


be determined by means of the foregoing) are 


known, the equilibria of the new mixture can be NoTATION 
~alculated from these equilibria and the mixing 


ratio. This is achieved by summation of the ¢ = concentration in the liquid ML? 
products of the volume fraction of one exchanger cO — for definition see equation (21) 


in the mixture and the load of that exchanger for D,() = break-through capacity at equilibrium te 
a given composition of the water. J (x) concentration in adsorbent at equilibrium ML 3 


lissociation constant 
Mixing of the exchanger may be an important X = dissocia on 
! = length of column L 
aid for obtaining a bed with a certain load (e.g. S = cross-sectional area of column 12 
those of Fig. 10(b), but not those of 10(a)), or for ¢ = time T 
replacing a good exchanger that is no longer ! linear velocity LT! 
ace ve 
commercially available by a mixture with similar = space velocity of liquid “"t 
wonertis J distance to the point where the liquid enters 
& 
I the column L 
‘The advantage of the load condition shown in x equivalent fraction of a Liven ton im solu- 
Fig. 10(b) over that of 10(a) is that the regenera- tion = exy,+/a 
: tion of the first column, due to the homogeneous y ~ free space per unit length in column L? 
loading, proceeds more efliciently than the 
NDICE 
regeneration of the second. With the condition ious 
illustrated in Fig. 10(a) there is even the possibility b = inflexion point in the isotherm 
break-through 
that an ion regenerated from the upper layer may 
r tangent point 
settle again in some deeper layer (cf. Figs. 6(a) 


sz = strong acid 


and 6(b), [1)). 


0 condition before the liquid enters the column 
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Abstract—The heat input in a distillation column which is needed to obtain a given separation 


of the components in the feed for an infinite number of trays in both sections, the stripping section 


and the rectifying section of the column, was calculated in a previous article. The reflux varied 


from tray to tray. In the present communication the actual and the generalized liquid and 


vapour rates as well as the temperatures have been calculated for that separation, using the 


methed of generalized component rates. 


As a check the total generalized vapour or liquid rate has also been calculated and found to 
remain constant within 0-1 per cent whereas the actual rates vary up to 40 per cent. 


\ calculation for a lower heat input than that corresponding to minimum reflux has also been 


carried out. One then finds that the top section leads to a temperature of the feed tray which 


is below the one which is calculated from the bottom section. Thus no common feed tray for 


the two sections exists in this case. 


Résumé—-La quantité minimum de chaleur nécessaire pour effectuer une distillation a reflux 


variable dun systéme de six composants a été calculée dans une communication précedente. 


Ici les auteurs calculent les courants de vapeur et de liquide pour tous ces composants et montrent 


que quoique leur somme varie considérablement dans la colonne (jusqu’é 40 pour-cent), la somme 


des courants généralisés reste constante dans 0-1 pour-cent. 


Les auteurs ont aussi fait le calcul en utilisant une quantité de chaleur inférieure & la quantité 


minimum, Dans ce cas en effectuant le calcul dans la section de téte, on aboutit & une température 


pour le plateau de charge qui est inférieure 4 celle qu'on obtient en faisant les calculs pour la 


section de queu. Dans ce cas évidemment ces deux sections ne peuvent pas avoir un plateau 


commun, 


Zusammenfassung—In einem friiheren Aufsatz wurde der Wirmeaufwand einer Destilla- 
tionskolonne berechnet, der fiir cine bestimmte Trennung bei unendlicher Bodenzahl in beiden 


Teilen, dem Abtriebs — und dem Rektifizierteil, erforderlich ist. Der Riickflass war von Boden zu 


Boden verinderlich. In dieser Mitteilung werden die wirklichen und die generalisierten Fliissig- 


keits - und Dampfstréme und die Temperaturen fiir diese Trennung berechnet unter Benutzung 


der Methode der generalisierten Komponentenstréme. 


Zur Priifung wird auch der gesamte generalisierte Dampf — oder Fliissigkeitsstrom berechnet. 


Man findet ihn innerhalb 0,1° konstant, wihrend die wirklichen Stréme bis zu 40° variieren. 


Auch fiir einen geringeren Wiirmeaufwand als dem, der dem Mindestriickfluss entspricht, 
wurde eine Rechnung durchgefiihrt. Fiihrt man diese Rechnung fiir den Kopfteil durch, so 


erhilt man fiir den Aufgabeboden eine geringere Temperatur, als wenn man die Rechnung fiir 


den unteren Teil durchfiihrt. Das bedeutet, dass in diesem Falle kein gemeinsamer Aufgabeboden 


fiir die beiden Kolonnenteile existiert. 


Minimum variable reflux for distillation of a six constant. Constant relative volatilities are 


component system was calculated in a previous assumed [1, 3). 


communication [4]. The calculation was _per- In the present article the liquid and vapour 


formed by a method which reduces the caleula- rates of the components, and the temperature at 


tion of a distillation in which the reflux varies various trays, are calculated for the separation 


from tray to tray to one in which the reflux is at minimum reflux. The absorption factor of the 
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heavy key component at the feed tray 
is determined for minimum reflux con- 70007; 
ditions. 


Definitions of Ap(m) and Sp(n) 
The method used has been set forth 5000 
elsewhere [2]. It is based upon : 
(a) the introduction of “ generalized — 
quantities” which enable us to treat var- ” 
iable reflux mathematically as constant 3000 


reflux, and on, | 
(b) the calculation of the absorption 


factor product of a reference component, 


here the heavy key component, in the 10008 


stripping section and of the stripping +92 | 
fac luct in the ecti 42500 
actor product in the top section. Oo5 © 6 20 25 36 35 4 @ 
1 

(a) It has been shown elsewhere [2, 4] /k 
that the expressions obtained for constant Fic. 1. Molar enthalpy of vapour and liquid in the column 
reflux remain valid in the case of vari- as a function of K and of 1/K respectively. A is the 
able reflux if the actual liquid and \ apour equilibrium constant of the reference component, C and R 


refer to condenser and reboiler respectively. The pinches 


rates are replaced by “generalized rates’’. 


1 are indicated in the curve for the liquid enthalpy. 
Io use generalized rates the molal en- 


thalpy of the total liquid rate H, must 
be represented by has its own approximating straight line and, 
H, = +H, /K (1) therefore, its own set of coetficients A,. (Fig. 1). 
For the calculation of minimum heat flux, 
and the molal enthalpy of the total vapour rate which quantity, if the reflux is variable, corres- 


H, must be written as : ponds to minimum reflux in the case of a constant 


H, = + /K (2) 


reflux, only the column section between the 


two pinches is of importance. 


The H’s in the right hand sides of equation (1) 
and (2) are constants. The generalized component (b) The absorption factor product is defined as 


rates are defined in both sections as : 


Ap(m) = A,. A,. 4y...... (6) 


= A, L; (3) 
Vio = AV, (4) with A, = LV, K, the absorption factor of the 
reference component at the kh" trayt 
The factors A, are constants and In the top section S, A,* V, Ay Ly is 
A, = — + (5) the stripping factor of the reference component 


at the A tray and the stripping factor product is : 


Equation (1) requires that the curve representing 
H, against 1. A must be approximated by straight Sp(n) = 8S, .S,.S,..... S,, (7) 


lines. It is natural to divide the column into three 


regions, separated by the pinches.* Each region + A list of notations is given at the end of the article. 

The trays are numbered upwards starting with m <— 0 

*The pinch of a column section is the infinitely long for the reboiler and m V 1 at the feed tray. They 

part of the section in which the liquid and vapour rates are numbered downwards in the top section with n 0 
are constant from tray to tray. for the condenser and n N + 1 for the feed tray. 
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Table 1. Separation scheme of six-component systemt 


Liquid feed Vapour feed F 


Component 


B; D a; 


Methane 1 | 0-010 0-250 0-260 sas 0-260 | 20-6 


Ethane 2 0-012 0-078 0-090 0-090 5-07 
Propane 3 | 0-070 0-180 0-250 0-004 0-246 2-06 
Butane 4 0-075 0-095 0-170 0-167 0-003 1-00 


Pentane 5 | 0-072 0-038 0-110 0-429 

Hexane 6 0-095 0-025 0-120 0-120 . 0-206 
0-334 0-666 1-000 0-401 0-599 


+ The relative volatilities are average values estimated from different sources for a temperature of 80°C and a pressure 
of 20 atm. 


2. The stripping section between pinch and reboiler O 04 O08 2 


The absorption factor product Ap(m) can be P| \ 
expressed as 
J M5 ++ 4 
Ap(m) = 2; >," (8) 
The reciprocal quantities ¢;' are the poles and M10 jh 


the coeflicients are the corresponding residues of 


the analytic function : 


V (1 Ay t) 


The summations must be carried out over all | 
components [2]. In the case of minimum reflux 
only the components which are specified in the «! iw | 
residue are present in the section between the | v\ 
pinch and the reboiler; the poles ¢,' can be S'y 
calculated straight forward from the denominator o OT 02 03 0465 06 67 06 


when the value of L,,, is known. In the present 


(b) 


Fic. 2. Liquid rates L,;,, in the stripping section. The 


numbers refer to the components ; 
{the light key-component. / 
varies, whereas the generalized liquid rate L., remains 


case four poles are to be considered, situated at 
t d ‘i 8, 4, 5, 6 in order of decreasing 
volatility, 
1 and 2 are unspecified in the residue. L,,) 


h is the heavy key, 


is the total liquid rate, which 


is the smallest value. Components 


constant. 


follows from : 


Ly = Hy, + 2 
boiler assuming equilibrium vaporization are : 
(Ref. [4], equation (12), p. 81) S, «; (11) 
In the previous article H, was found to be 


The generalized vapour rates from the re- 


B, (10) 


(0) 


105, thus L,,, = — 105 + 5563 — 5458 [4]. 
The resulting ¢,’s are given in Table 2. The 
residues ¢; follow from : 
Vie %;) 


lim (1 — ¢,t) A(t) = 
= lim ( t) A(t) 


1 
1,000, 


and 


in which A, is determined by the equation : 
l= 2 Kya, B; (12) 


or Ky, = 1-623, and S, is the stripping factor of 
the reference component in the reboiler or 
K Vio/ =; B;. 


The results are given in Table 2. 
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Table 2 


17700 


15-42 2-822 1 84-76 


& 


334-2 120100 T4585 OSS415 
5817 1427 56539 1057 
3053 0-20263 953-35 O-O554 OoT47S 


The absorption factor products can now directly The accuracy of the computations has been 
be computed by inserting the values of the ¢,s tested by computing LT V,,,.. which must equal 
and ¢.’s into formula (8), Law 2402. This condition is 


Since, satistied within 0-1 per cent for each tray from 


Ap(m) ‘ 

A. (13) tom 13. The equilibrium constant of 

tn the reference component K, ,, is computed from: 

A m (I A tm A 2; 17) 


the absorption factor of com pone nt 2 on trav m From the known graphical relation bet we en the 


temperature ¢ (degrees Centigrade) and kK, om 


can be determined. For the computation of 


the temperatures have been derived. Some 


the actual liquid and Vapour rates we use the Sa 
law of conservation of matter : results are given in Table 3. 


It can be proved that the absorption factor of 


B, Vow (15) the reference component in the pinch equals 


4, [2 
and the equilibrium condition ; (2). 


I 


‘im (16) 3. The top-sec tion between pinch and condenser 


The B;s are known, the A..’s have been This part of the column can be treated in a 


calculated, and the V's are easily found from quite analogous way. The stripping factor 
the V... since V 1" 1 In this way the product can be written as : 


“0 


actual liquid and vapour rates can be computed 


Spi n) Cy (18) 
step by step from the first tray upwards until 
at m 13 the values at the pinches are sufficiently (The suflix ¢ will be used to denote a quantity 
approximated (equation (51) in this article). pertaining to the top-section.) 


Table 3. Quantities pertaining to the stripping section below the pinch 


Ap(m) = | K 


m m m (m) 


l O-0549 1-33 134-6 24014 
2 1-239 1-212 1-102 O-7559 120 125-5 2401-3 
2-559 1510 1-155 7226 144) 115-3 240161 
6 7 1-771 1-084 00-6463 106-6 2400-8 
27-25 O-ST4 100-7 2400-6 
10 110-8 2-033 O-S40 


e740 2072 9719 05693 O-s24 97-0 


l 
2 
“ 105 
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The ¢,' are the poles and the ¢, are the stripping factor product (14). Equilibrium con- 


corresponding residues of the analytic function: densation is assumed in the condenser thus : 


) (a) 


S(t) 


in which A, follows from the equation : 


One has: 


and 
The expression K, 2, D,; is the ab- 
Vi) = Hp — 2, +H, 2) D, (21) 
. : sorption factor of the reference component in 
which gives V,,, = 3827 + 1739 = 5566 [4]. the condenser. 


From 


| 


N-§ N-9) 

) 

O =~ 6-4 0-6 8 Fe) 


Fig. 3. Vapour rates V,;,, in the top section. The numbers 


refer to the components. VV, is the total vapour rate, 


which varies, whereas the generalized vapour rate Ve) 


remains constant. 


The most important quantities are given in 
Table 4. 
The stripping factor is computed from the lic. 4. Temperature distribution in the column 


Table 4. Generalized distillates, roots and other data for the top section 


Comp A D t x;] D; $1; L «0 


1 10375 2698 2256 0-090236 710-1 0-14866 
2 2184 196-6 + 22-323 0214427 710-4 0-08613 
1878 486-6 0-603719 1215-8 0-71129 


11-43 + 7-761 100722 61-09 0-05409 


100017 


| | ; 
| 
4 6 } 1-2 Sin x, S,, (24) 
~ - 
P in-l “in i 
and 
Vin Li. (26) 
> \ 
\ 
\ 
| | | 
4 3811 
1739 2198 
LC 
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Table 5. Quantities pertaining to the top section between pinch and condenser 


L 


Sp(n) 
1-000 
0-5158 
0-31938 
01504 
0-0909 
0-06985 
0-0627 
0-06044 
0-06126 


0-8282 
00-9805 
00-9231 
0-8268 
(7404 
7077 
0-6835 
00-6805 


0-5157 
00-6192 
7072 
0- 8032 
00-8993 
0-9606 
10005 
1-0006 


0-599 

14274 
15795 
15639 
14741 
1-3828 
13237 
1-2851 
1-2801 


0-255 
0-354 
0-389 
0-417 
0-451 
0-487 
0-514 
0-532 
0-535 


2198-3 
2198-1 
2197-7 
2196-2 
2197-9 
2198-0 
2197°5 
2198-0 
2198°5 


the actual liquid and vapour rates can be cal- 


culated. Finally the temperatures are derived 


from A, ,, K, , being given by 


Some results are given in Table 5. 
It can be proved that the stripping factor of 


the reference component in the pinch equals 


du (2). 


4. The absorption factor A, at the feed tray 


The numbers M, N, the 
containing the feed tray are all infinite in the 
case of We still 
however, the feed tray by M + 1 as belonging to 
the stripping section and by N 


m and n in region 


minimum reflux. indicate, 
1 as belonging 
to the top section, but introduce the symbol « 
to denote a tray in the stripping section defined 
by m = M — x» and in the top section » defined 
byn=WN 
calculations in this region are much more com- 


vy is used to denote the trays. The 


plicated owing to a discontinuity in the total 
generalized liquid and vapour rates at the feed 
tray and owing to the presence of the unspecified 
components. There is also a discontinuity in the 
absorption factor. The actual absorption factor 
A, at the feed tray differs from Ay, as cal- 
culated from the absorption factor product in 
the bottom section. Using the index u to denote 
the unspecified components in the bottom section 


(u = 1 or 2) we can write : 


L ium +1) V um) = @ (28) 


2 


because in the case of minimum reflux the un- 
specified residues are zero. Below the feed tray 
one has: 


L : A. V, = A 


(um) um (um) V cm) (29) 


where A,, is the absorption factor for the reference 
component as calculated from the absorption 
factor product in the bottom section. On the 


feed tray one has : 


Dig = Ap Vig (30) 


As has been stated above A, is not equal to 
Ay, , Ap(M+1)/Ap(M) but one can prove 
that : 

APM +1) 


81 
Ap M) V (31) 


(ft) 


From these four equations it follows that : 


+1) 


For the specified components we have : 


m 


[ %,"* Ap (k) + 


k=0 

| /Ap(m) (33) 

Combining equations (32) and (33) we can de- 
rive : 


Ap(m) ip ( (b), 
(M—m+1) 


iuN +1) Ly 


+ { Aplk) + Vio) (34) 
k=O 


30 


| Ve | Ke (°C) 
48-8 
53-6 
4 | 57-5 
6 | | 61-5 
| 65-7 
10 | 68-3 
4 | 70-7 
| | 
VOL. 
: 
Ls 


The sum over i is to be taken over all the 
specified components [2]. In these expressions 
Ap(m) means the absorption factor product of 
the reference component in the bottom section, 
in which the contribution of the unspecified 
components is taken into account. Actually 
equation (8) can only be used to calculate the 
absorption factor product Ap,(m) in which the 
contribution of the unspecified components is 
ignored since the quantities ¢; in equation (8) 
cannot easily be calculated for an unspecified 
component. 

The relation between Ap(m) and Ap,(m) is : 


Ap(m) = Ap,(m) + Ap(M +1) 2, ¢,%, (M-m+1) (95) 
in which the constants c,, are : 
Vp) Viunst) Vin 


mA) 


(36) 
Li) - u 

Because there are an infinite number of trays 
between the region under consideration and the re- 
boiler, we must take the limiting value of equation 
(36) for m + o which makes c, independent of 
the tray number m : 


| (b) J (uN+1) (t) 


37 
Li) 2) ') 7) 
From equation (35) it follows that : 
Ap(M +1) = (38) 
and 
Ap,(M +4 
Ap(M) = Ap,(M) + %' (39) 
=a Cu 


From equation (8) we find : 


+ 3) 
m>o Ap,(M) 


where ¢, is the largest of the ¢; occurring in 
equation (8). Finally, with the use of equation 
(27) we find: 

Vi») $3 


(1 $y %, ') (41) 


= $5 (40) 


Similar reasoning for the trays above the feed 


tray gives the equations : 


Sp(n) = Sp,(n) + Sp(N +1) * (42) 
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Here the double index tu refers to the un- 
specified components in the top section (tu - 
5 or 6); also the relationship : 

— 2; Dey — 
and 

Lit) dis 
Kquations (37), (41), (43) and (44) enable us to 

compute A, by a repetitive procedure. The zeroeth 


approximation consists in taking D,,, instead of 


(u) 
the still unknown V,,, y,,) in equation (37). With 
these approximate values of c,, the first approxima- 
tions of Ly, 47.) are computed. Here it may be 
mentioned that Ly, .7,,) is the liquid rate of a 
component, unspecified in the top section, but 
specified in the stripping section; V,,, y,,) is the 
vapour rate of a component, unspecified in the 
bottom, but specified in the top section. 

For the components, specified in the bottom 


section we have : 


Ly M +1) 
1—Z,¢ 
4 be Cu 
B.) l—«@, % 45 
¢, (1 — $5 


and for components, specified in the top section : 


View 


With the improved values for Lyn, 7,,) and 
Viunsi) We get the next approximations of c¢, 
and ¢,, and so on. To judge the progress of the 


approximation we calculate at each step the A, 


tu 


according to equations (41) or (44). The results 
are given in Table 6. 

Thus the two values of the absorption factor 
A,, calculated from the top section and from the 
stripping section respectively, are equal. 


5. The liquid and vapour rates in the region con- 
taining the feed tray 

From the basic equations (28) and (29) and with 
the aid of the equation (35) for the absorption 
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Table 


720° 


2605 

2 504 
3 

4 9506 OS1-6 
5 


from 
eq. (44) 


O-05646 


1-120 
O-06550 1-130 

O- 2325 O5162 1-131 
O-06 560 1-151 

2323 1-151 


factor product we can express the generalized 


liquid rate of the components unspecified in the 


bottom section at the tray m l M l Tv 
in the corresponding liquid rate from the feed 
tray 
(6,2, °°) .L uM (47) 


The L,,, 47. ,,'s are computed from Table 6 with 
equation (30). From equation (28) it can be de- 
rived that the liquid rates of the components 


specitied in the bottom section are given by 


(48) 


The equations also hold for the actual (non 


weneralized) rates. Similarly we have for the top 


section : 
=m Cm (I Dis 
‘ 
(Dy, %,,) (40) 
and 
iN 
l 
— ty Cony (4 
Pry Lon (50) 
5 
| 


If we take the limiting value of equation (48) 
lor « oo and of equation (50) for v a we find 


the expressions for the rates at the pinches : 


B 


(51) 


4 


i, pinel x; by 


1" (52) 


i, pinch 
l %; Pig 


The absorption factor Ay enw becomes : 


(ds *)*] 


(33) 


in the bottom section. In the top section we 


have : 


(54) 


The equilibrium constant of the reference 
component Ay (heavy key) is calculated from 
equations (17) and (27). From A, the temperature 
IS derived. 


6. What happens at a reflux lower than the minimum 
reflux ? 

The absorption factor of the heavy key com- 
ponent at the feed tray has also been determined 
at a slightly lower heat input corresponding to a 
reflux below the minimum value. In this case one 
obtains two different values of Ay according to 
whether the calculation is performed in the top 
section or in the stripping section. The stripping 
section leads to a lower absorption factor than the 
top section; this result means that no common 
feed tray exists. In the case of minimum reflux 
the*net upward heat flux from the reboiler is 


i, 105 and the net upward heat flux from 
the condenser is H,, + 8827. The results for 
Hy, 520 and H,, — 3412 are shown in Table 7, 


It is not necessary to compute at each step of 


approximation both sets of constants ¢, and ¢,, 


Step 

eq. (41) 

VOL 
Sway 
Ou 
l Ps bed 
l e, [1 (Py } 
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Absorption factor below minimum reflux 


| 
Step | Tes | Mewes) 


1 1035 2582 9553 
2 2565 9558 
3 1018 2504 


| 

| | 
‘u=2 Cts from 
eq. (41) 


Ip from 
eq. (44) 


2275 0-532 0-9166 12551 

OOT1ISI 0-2637 0-5692 1-031 1-118 
166 O-2506 00-5652 1-033 1-133 
O2505 0-5651 1-133 


Scale enlarged 


| 
~ 
M-2 M-1 M M41 W P N-2 N-3 
+ 
Tray number 
Fic. 5. £ and » near the feed tray 
Ly, Ap(M a) Lig l v) 


At the feed tray the relation ¢ 


equivalent to the requirement that the equilibrium 


7 must apply. This is 


constant of the reference component K at the feed tray 
must have the same value when calculated from the bottom 
section (£) as from the top section (»). The full drawn 
curve pertains to minimum variable reflux ¢ 7 at 
point 4. The broken curve pertains to a slightly lower 
reflux ratio and K remains lower in the top section (point 
B) than in the bottom section (point C). Thus no common 


feed tray exists. 


respectively and both values of A,. It has only 
been done in this case to demonstrate cl arly the 
rapid convergence of the calculation process for 
the generalized rates L,5 Van 
Vien 41» for the constants c, and ¢,,, and for the 
absorption factor A ». (Compare Table 6). 


Ap(m 1) 


In Fig. 5 the functions 
Vi. Sp(n) 
Ly Sp(n + 1) 


Ap(m) 
have been plotted as a function 
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of m M pw and n N I respectively. It 
has been shown previously that at the feed tray 
for which m M l and n N 1 these two 
expressions must be equal [2]. They are denoted 
by € and by » respectively when referring to the 
feedtray. Thus & = » must apply at the feed tray. 

This requirement has a direct physical in- 
terpretation. In the case of a constant molal 
reflux the meanings of 
Ap M 1) L, Ap(M).Since Ap(M +-1) 
Ap M) L,/K i.e., 


the reference component for the feed tray but cal- 


and of » are obvious : 
the absorption factor of 
culated ignoring thi discontinuity in liquid and in 


1K. Similarly LAK 
but now the calculation has been performed from 


vapour rate, one has é 


the top section. Thus the relation é 7 Means 
that one must arrive at the same value for the 
equilibrium constant of the reference component 
at the feed tray by performing the calculation from 
the Stripping section as from the top section. If 
the two values are not equal no feed tray is present. 

In the case of a variable reflux a special proof 
is needed (but the proof is omitted here) to show 
that : 

Vis Vs Lis) 


and 


Vw mer Ly, 
means that the requirement & 


and this relation 
7 at the feed tray 
still means the equality of the equilibrium 
constants calculated from either section. 

In Fig. 5 one sees that below minimum reflux 
» iS larger than & or expressed in words: the 
temperature in the top section always remains 
below the lowest temperature in the stripping 


section, 
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NOTATION 


1 function of a vuriable ¢ defined by the | Pe Pe total liquid rate from tray m or nm res- 
relation 1 1+ Ap(i)t ipi2) t* pectively 
ip(m) L, total liquid rate in top section 
Ip absorption factor of the reference com- Ly total generalized liquid rate in top 
ponent on the feed tray section 
1,, = absorption factor of the reference com- / subscript indicating light key 
ponent A,, on trav number factor defined by the relation 
, :bsorption factor of component i on tray The component ratesare transformed into 
generalized rates by multiplication by 
ip(m) absorption factor product for the reference I; 
component at tray number m M molecular weight 
ip(m) = A, Ag.... Am m index for tray number in bottom section 
4, = Telative volatility of component ¢ in wu — index used to denote trays near the feed 
respect to the = reference component trav in the bottom section m VW mm 
x AA n — index for tray number in top section 
B, residue rate of component 4 ' index used to denote trays near the feed VOL 
R generalized residue rate of component 1, tray in the top section n \ \ 5 
defined as B 1, B; Ry reflux ratio in condenser 
Oo subseript referring to the bottom section, Rp vapour reflux ratio in reboiler 
i constant occurring in the expression of &.., &,, stripping factor of the reference com 
the absorption factor product Ap(m) ponent on tray m or n respectively 
a he oa Sp(n) stripping factor product for the reference 
D, = distillate rate of component t component at tray number n. Sp(n) 
generalized distillate rat f component 
i, defined as Di 1 Dy Sin stripping factor of component t on tray 
F; feed of component i, F; B; 
a reciprocal root of denomimator in the pH sommation over all components 
ut quantity belongs to the top section 
Hp net upward heat in bottom section ete.) 
ttn tu’ 
Hp net upward heat Sux in top section t mathematical variable used in the ecal- 


i, molal liquid heat content 
culation of the absorption factor product 
constants occurring in the « \pression 
molal liquid heat content as a linear 
function of 1/K " subscript to indicate unspecified com- 


il molal vapour heat content ponents 
i, 0 i, 1) constants occurring in the expression of Vey) total generalized vapour rate in bottom 
molal vapour heat content as a linear section 
function of A Fas Vn total vapour rate from tray m or n 
h subscript to indicate heavy ke respectively 
i subscript to indicate a specified com- j sil r be vapour rate of component i from tray 
ponent m or n respectively 
K equilibrium constant of the reference 
on) V tim generalized vapour rate of component i 
comp me nt (tm) (in 
from tray m or n respectively 
K,, = value of K on tray m : 
K.. equilibrium constant of component i on Vin yeneralized total vapour rate in top 
section : 
tray 
Le total generalized liquid rate in bottom Vin, Ap(M +1 p) 
section Lip) ip( M p) 
L,.., liquid rate of component i from tray m 
im? “‘in ppt N v) 
or n respectively (t) 
generalized liquid rate of component ( OFA r) 


“*( om) 


at the feed tray £ = » L/K 


from tray m or m respectively 
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deal steady-state stagewise countercurrent distillation, liquid 


The concept of a carrier” 


a cifference equation whi h describes all four processes This 


equation may be solved by the transformation methods de veloped for distillation processes 


The assumptions made in the derivation of 


worked out 


Résumé l-auteur déerit une théori generale 


the equation are discussed and an example is 


relative aux opérations fondamentales suivantes 


eflectuces a contre courant. con ces CLage par ctage, dans un état de régime idéal : distillation 
extraction liquide liquide uiserption et absorption Le concept d'un vehiculeur est 
developpeé et utilisé pour obtenir u equation differentielle qui rend compte des procéedés. Cette 


equation peut étre resolue par les 


distillation. Lo auteur discute les | 


Zusammenfassung cis 
Destillation 
mitgeteilt, ks wird cic 


fliissiy-fliissig 


Vorste ul 


pothéses faites en dérivant les équations, et donne un « sample. 


stationire stufenweise, im Gegenstrom 
rption 


emes Triigers 


ethodes de transformation dévek ppees pour les procédeés ce 


verlaufende 
Adsorption, wird eine allgemeine Theoric 


entwickelt und zur Ableitung einer 


Differenzeny lk ihuny benutzt welche alk Vier Proz besehrenbt ie se Glek huny kann mit 


Hilfe von Transformationsmet hod velost werden, wie sie fiir die Destillation entwickelt sind 


Die bei der Ableitung der Gleichung gemachten Annahmen werden diskutiert und ein Be ispie' 


durchgerechnet. 


Because the raw materials for. and the products 


SECTION INTRODUCTION 


from, industrial chemical reactions are scldom 


pure substances, s« paration processes are of 


interest to the chemical enginecr, Important 


among such processes are those in which matter 


diffuses 


mutual 


between two immiscible phases in 


contact, e.g. distillation, liquid liquid 
extraction, adsorption, and absorption processes. 

In the design of such processes it is nece ssary 
to calculate the number of ideal stages, or the 
number of transfer units, required for the separa- 
tion, and often, as a preliminary to these cal- 
aleulated. 


All these calculations become tedious when the 


culations, the minimum reflux ratio is 


number of stages or transfer units is large and 
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when the mixture to be s parated contains more 
than two compone nts. 
For distillation 


elegant and less tedious methods of calculation 


processes, howeve r. some 


have been reported [1, 2): similar methods do 
not appear to be available for the other processes, 
This fact is rather surprising when the close 
similarity of the four processes is considered, and 


it prompted the work reported in this paper. 


Secrion I] 


In order to ce sign any of the processes con- 
sidered here the equilibrium data for the appro- 
priate two-phase system is required. Such two- 
phase systems in thermodynamic equilibrium are 
characterized by a similarity which makes it 


pp. 226 t0 263. Pree Lid. 
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possible to describe them all in the same terms, 
i.e. both phases in such a system are saturated. 
To make this clear consider first the two-phase 


system for a distillation process. At thermo- 


dynamic equilibrium addition of heat to the 


liquid phase would cause the formation of some 


vapour, and similarly removal of heat from the 
vapour phase would cause the formation of liquid ; 
with 


in this sense both phases are saturated 


heat. 


For a liquid liquid extraction process, conside r 


the system consisting of two liquid phases in 
Addition of solvent 
the 


thermodynamic equilibrium, 
to the 
formation of 


solvent—lean phase would cause 


some solvent-rich phase, 


similarly the removal of solvent from the solvent 
rich phase would cause the formation of some 
solvent-lean phase. In this sense both phases are 


saturated with solvent. The analogy between 


specific enthalpy in distillation solvent 
content in liquid-liquid extraction is well known 
and quoted in several textbooks ; [8]. 

For an adsorption process consider the system 
consisting of the fluid and adsorbed phases in 


Addition of ad 


sorbent to the thuid phase would cause the forma 


thermodynamic equilibrium. 
tion of some adsorbed phase, and similarly the 
removal of adsorbent from the adsorbed phase 
would result in the formation of some fluid; in 
this sense both phases are saturated with ad- 
sorbent. 

Finally, consider an 


absorpt lon process in 


which only one substance, A, has, when pure, a 
vapour pressure greater than the working pressure 
of the 


the vapour pressure of the liquid is equal to the 


process. At thermodynamic equilibrium 


working pressure, and addition of A to this 


liquid would result in the formation of some 


Table 


vapour, similarly the removal of A from the 
vapour phase would cause the formation of some 
liquid. 

Thus for each of the processes considered there 
is a special substance, e.g. heat for distillation, 
characterised by two properties ; firstly this 
substance is necessary for the production of two 
phases, and secondly both phases in the process 
are saturated with it. It is convenient to refer 
to this substance by the same name for all the 
called the ** 


It is also convenient to refer to the carrier—rich 


processes ; it will be carrier.” 
and carrier-lean phases by the same names for 


all the 


carrier lean 


pre W'CSSES, ; Carrier rich and 


phase 


phase’ are rather cumbersome 


terms. therefor -phase and L-phase ” 


‘ 


The terms “ carrier,”’ 


defined 


re spectively will be used. 


phase and phase are for 


each process in Tabk 


Sectrion III 


Consider a portion of an ideal countercurrent 


cascade, consisting of several 


consecutively 
Let l 


the flow of j™ comporn nt im the L-phase leaving 


numbered stages (see Fig. 1). ., denote 


the s + 1™ stage (moles per unit time), and let v 
the flow of the y= 
phase leaving the s™ stage (moles per unit time). 


If one assumes that no accumulation of the j™ 


j 
component in the V- 


denote 


component occurs in the cascade, then application 


of the conservation principle gives : 


u; (1) 


Here u; is a constant, the net upflow of the j™ 
the 


carricr, 


component. Suppose mixture contains n 


the 


equation (1) for all the components In the mixture 


components plus then summing 


Process Carrier 


Heat 
Adsorbent 
Solvent 


Distillation 
Adsorption 
Liquid liquid extraction 
Absorption 


working pressure 


Substance with vapour pressure 


V-Phase L- Phase 


Vapour 
Adsorbed 
Solvent—rich 


Liquid 
Fluid 
Solvent—lean 


Vapour Liquid 


| 
| 
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except the carrier, and defining V-phase and A sequence of stages for which the u, and K are 


L-phase flow rates : constant is called a “ cascade section.” 


Knowledge of the equilibrium data for the 


2 v;,, and L | one obtains, to f 
ivs - i system enables one to write further equations, 


| L, VoL which are assumed to have the form : 


«C.y;., (3a) 
(For simplicity Z is abbreviated to in this 
— Sea (3b) 
miper). ” ” 
pal / (4) 
If one defines mole fractions 
Phe x, €;, and the are all constants in- 
Zo, Jel depend nt of composition. The are separation 


; factors (relative volatilities for a distillation 
then one may write equation (1) as 


process). Putting 2; =y; =1 in equations 
Vy... (3a) and (8b) it is seen that and ec; are the 
- carrier contents of saturated phases containing 
I} u terms change only whe re material Is th 
only the 7 component and the carrier. 
dded to or removed from the cascade . 


The correspondence between equations (3a), 


Now consider the carrier, writing ¢ for the (3b), and (4), and the real systems which they are 


ontent of the |-phase leaving the s“ supposed to describe, is discussed in the next 
1.¢ unit ol carrie! r mole ot the other n Section. 
nonent the mixture), and writing c,, , for Using equations (3a) and (3b) we can write 
ontent of tlhe i phase le iving the equation (2) as: 
|" stage (units of carrier per mole of the other 0, Te] 
componel in the mixture); assuming no U;., =" 
mulation of the carrier in the cascade. 


conservation principle gives 


Rewriting this equation one obtains : 


(2) 


and using equation (1) to substitute for the 


A. the net upflow of carrier « hang only when 
bracket on the right hand side of this equation, 


is added to or removed from the cascade. 


one obtains : 


aC; U;, to) 


The right hand side of this equation is a constant 


for a cascade section. Now put 


(C 


' 


¢,) = 6, and b, b, B, 


then equation (5) may be written : 


‘ 

“'v;., B; = a constant (for a cascade section). (5a) 

Stage:s 

Similarly it can be shown that : 

21, ,., B, = another constant (for a cascade 
Stage:s-2 section). (5b) 
Since B; is a number, v, , B; is a flow rate 
MG. 1. Flows in part of a cascade section (B; moles per unit time), which will be called a 
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— 


“ pseudo-molar flow rate.” Pseudo—molar phase 


flow rates will be defined as : 
we 
Zv;,,B,andL= Zl,, 
The subscript s has been dropped from V and 
L because they are independent of s (see equations 
(5a) and (5b)). Similarly pseudo-mole fractions 


may be defined as : 


B,v;.,/2 B,v,., 
B, 
Equation (1) may be rewritten as : 
me l B 
u. B. 


When expressed in pseudo-molar quantities 


this equation becomes simply 


Vy;.. — (6) 


Here B,. 


Equation (4) may be transformed into pseudo 
molar fractions without any change in the separa 


tion factors, viz.: 


(4a) 


s 
Using equation (4a) to substitute for Yj. in equa- 


tion (6) one obtains : 


If equations (3a), (3b) and (4) apply, equation 
(7) deseribes all the four processes considered in 
this paper, and it can be solved by the methods 
referred to in the INTRODUCTION. 

So far the numbering of the components has 
not been specified ; the usual convention will be 
adopted here, i.e. components will be numbered 
in order of increasing separation factor. 


Section IV 
The correspondence between equations (3a), 
(3b) and (4), and the real processes which they 
are supposed to describe, will be examined in 
this section. However, discussion will be limited 
to those processes for which it would be desirable 


to use the elegant methods of calculation referred 
to in the introduction, i.e. processes which require 
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a large number of stages. A large number of 
stages is required when the separation factors are 
small : the methods referred to in the Inrropuc- 


TION are applicable when the separation factors are 


constants. Thus only systems for which the 
separation factors are small and constant will be 
considered. Such systems will generally be 


ideal” in the thermodynamic sense. 

Distillation processes. Mixtures which display 
small and constant relative volatilities generally 
contain only substances having little tendency 
to molecular association or dissociation. In such 
systems the heats of mixing in the liquid and 
vapour phases will be small compared with the 
latent heat of vaporization; for such a system 
equations (8a) and (3b) will decribe the enthalpy 
of the liquid and vapour phases quite well. In 
any case equations (8a) and (3b) will provide better 
approximations than the assumption of constant 
molar latent heat which is usually made. 

The 


pre ViIous 


pseudo-molar quantities defined in the 
the 
which have previously been used for distillation 
calculations [8, 4, 5). 


Liquid liquid extraction 


section are similar to heat units 


processes. Two-phase 


result of excessive non- 


the 


liquid systems are the 
Nevertheless, if 


Se parated consists of similar substances with no 


ideality. mixture to be 


tendency to association or dissociation in the 


solvent, then equation 3a) should describe the 
solvent content of the solvent-rich and solvent 


Little data is available 


for multi-component this 
but the mixture used in the example is well 


lean phases quite well. 
mixtures of nature, 
described by equations (3a), (3b) and (4a). 

If the solvent 
then the equations obtained in section IIT would 


were a mixture of substances, 
only be valid if this solvent could be treated as 
one substance, i.e. if the molecular weights of the 
substances in the solvent were very similar, or 
if the composition of the solvent did not alter 
with position in the cascade. 

Adsorption processes. These processes are run 
so as to minimize attrition of the absorbent: 
thus the adsorbent content of the fluid phase 
should be negligible, i.e. all the “c,’ terms of 
equation (3b) will be zero. 


It is usual to assume that dynamic equilibrium 
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is achieved between the fluid and the adsorbed 
phases : that monolayer 
adsorption occurs then equations (3a) and (4) 
describe the system, as shown below. 


if one also assumes 


Let @, be the 


covered by the i component, 


fraction of adsorbent surface 


m; be moles of the i component per unit 


area of surface, and A the surface area per 


unit weight of adsorbent. 
Then : 
) 
0; =m; 0, 
and, if dynamic equilibrium is achieved : 


The number of moles of j*' 
the 


The numberof molesof jt” 
rly ad 


per unit tine 


component leaving component ads on 


surface per unit time the surfaces 


That is : 
Am; -(rate constant) 


(] 2 0.) (rate constant) 


Simplifying : 


m, 0. constant . 2 


Writing a; for the constant and substitutine 


for m;@; in the expression for y, written above 
one obtains : 


Equation (4) is obtained from this equation if 
the right hand sic Is divided throughout by a, 
Now 


and a,/a, Is put equal to x,. when a : 


then and (say), so that 
1/Am, and, forming we find 
1/A om; (since 2 0, 
but 1/4 Im, 
so TC, y; = C, ice. equation (3a), 


Thus equations (3a) and (4) deseribe th system 
exactly. 


Absorption processes. Consider an 
to 
Henry's law applies, so that : 


ponent mixture, each 


component of which 


D. A. More 
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Putting «, = g,/g, and ap 


P/g, we obtain : 
Lp 


Suppose that the absorption process runs at 


such a temperature that only g,,, > P, ie. the 
n I component is the carrier. The C, and 
the ¢; are found by considering the set of binary 
mixtures obtained by taking the components 
& (n 1), m. separately with 
the n jth component. For each of these 
mixtures one has : 


Solving the last two equations for aw, and x, 
one obtains : 

Using this and a similar expression for C, it is 
easy to show that equations (3a), (3b) and (4) 


describe the system ¢ xactly viz.: 


- 
“ 


For such systems then, calculation for ann + 1 
component mixture can be reduced to calculations 
for an n component mixture to which equations 
(3a), (3b) and (4a) apply ‘ xactly. 

When more 


pre ssure 


than one substance has a Vapour 
greater than the working pressure this 
method can only be applied if these substances 
can be treated as one substance, which is the case 
the 


(4) the molecular weights are equal, and (c) the 


if (a) the g. are 


equal for substances, 
ratio of the amounts of these substances in each 
phase does not vary with position in the cascade. 

Since the x, ¢,, and C; are almost independent 
of temperature in such a system, equation (7) will 
apply to adiabatic absorption columns. 

For a stripping process, in which only one 
substance has a vapour pressure less than the 
working pressure, a similar set of equations can 
be obtained ; remarks similar to those given in 


(a), (4) and (c) above also apply. 
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AN 
To work out an example, UNpERWoon’s trans- 
the 
resultant equations apply to all the processes 


SECTION V. EXAMPLE 


formation will be applied to equation (7) ; 


considered in this paper. An outline of the deriva- 
tion of these equations is given next. 

The composition of an n component mixture 
may be represented by a point in an n dimensional 
space, the co-ordinates of the point represent the 
mole fractions of the components of the mixture ; 
possible compositions all lie on the figure described 


The 


by # 


i 


pseudo-mole 


section II] are obtained by a transformation of 


this co-ordinate system, VIZ. : 


UnpDERWoop’s transformation is given by 


If the wu, are all positive, then : 


¢; 


if the u, are all negative, then : 


/ 


Cases exist in which all the u, are neither all 
positive nor all negative, but they will not be 
considered here. 


When all the ¢, the A, 
tained by solution of the set of equations, 


are found can be ob- 


2 A; a; /(%,; ¢;) 1 1, 2, 3, 


The number of ideal stages required for a given 
separation is obtained quite simply by solution 


of the equation : 
(¢, s) (9) 


When a liquid composition lies on one or more 
of the figures described by 0 an infinite 
number of stages is usually required in the cascade, 
the condition corresponding to minimum reflux. 
Minimum reflux conditions in distillation cascades 
have been thoroughly discussed in the literature 


[6, 7], and N. L. Frankuin and R. W. Goutcuer 


fractions of 


hope to publish another paper discussing single 
sectioned cascades. 

For distillation processes the wetness—fraction 
q has been defined because its use simplifies 
the 
described by equation (7) a saturation fraction 


certain calculations. For general cascade 
analogous to the wetness fraction can be defined ; 
q, the saturation fraction of a mixture is defined 


as follow ge 


The amount of carrier which must be added to con- 
vert the whole mixture into a carrier-rich phase. 
added to a 


carrier lean phase, containing the same weights of 


The amount of carrier which must be 


the other substances as the mixture, to convert this 


carrier-lean phase into a carrier—rich phase. 


This is the same as the wetness fraction for 
a distillation process. 

Example. McCase and Smirn [8] give the 
following example. A feed of 40 per cent methyl 
cyclohexane and 60 per cent n heptane is to be 
separated into products containing 1 per cent 


and 98 per cent of methy leyclohexane in a liquid 


liquid extraction cascade with reflux at both ends. 


Aniline is to be the solvent. 

Equilibrium data for the system, published by 
VARTERESSIAN and Fenske [9], together with 
other data calculated therefrom, is given in 
Table 2. All this data is required before the 
Jaxecke diagram used by McCare and 
can be plotted. 

From Table 2 


considerably in a rather unsystematic manner 


it can be seen that a, varies 
(probably because of difficulties in analysis) ; 
the arithmetic mean of « has therefore been used. 

C was plotted against y, and ¢ was plotted 
fitted to 
the intercepts made by these 


against a, straight lines were these 


graphs by eye; 


1, Ys 0, 
= 1 have been taken as ¢,, C,, and C, 


straight lines and r, 0, Be and 


respectively. 


0-075, Cy = 0-125 

c, and ec, are so small in comparison with C, 
and C, that they have been ignored in the 
calculation of B, and B,. Thus 6, 15 
b, 5, therefore B, 1 and B, 1 3. 


and 
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where d; is a solution of : j 
a, u;/(x%; — ff (8) 
< 
J 
= 
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Table 2 


Carrier—lean layer Carrier_rich layer Calculated values 


(wt. (wt. 


MCH n heptane aniline M.CH n-heptane aniline 


0-000 O-000 15-18 


9-2 83:1 7-7 Os 60 03-2 00-0996 1176 O-O8S4 13-70 1-82 
18-6 73-4 2-7 5-3 92-0 2021 O-3375 11-50 201 
24 69-8 30 | O-37038 11-4 1-86 
33°58 57 O-3608 0-5055 0-09040 1-74 


OSS 00-6727 O-0080 8-09 201 
OO04 O- 7666 7-33 1-49 


67-2 1! 2:1 86-6 OT Oil 651 1-92 
716 12-7 85:7 O-SS81 O-1185 202 
73 5 1160 


O-1274 


Now the feed and product specifications can be ratio R — 7-7. one first solves equation (8a) for 
transformed into pseudo-mole fractions two values of é obtaining : 
B, 2,.¢/(B,2,.¢ + Bato. p) = 0-9913 and 4, = 1-67. 


1. 0-6 0-6 + (1/8)-.0-4 0.82 


These values of 46 are used in equation (9) to 


Similarly, xy. , 0-18,.2,_,) — 0-058 and obtain the number of stages required between 


the feed and product compositions. Writing 


Ts. y 0-0034, 


equation (9) in full one obtains ; 


Using the definition of q given above, 2, 2,. $1) 2%», p/(ae 


The minimum reflux ratio will be caleulated %: (% hy 
first. The reflux ratio R LN L). Dividing % — be) p/\% $2) 
‘ 
equation (8) by 2 u,, one obtains 
& G; Di \F ? ony Noting that A, / Ae cancels from both sides of 
VV L) R l (Sa) the equation, and that the large fraction on the 


right hand side of the equation is equal to one, 


This equation gives the minimum reflux ratio 


since both the numerator and the denominator 


when one uses the value of 4, obtained by solution 


of : 


are equal to R + 1, this equation ts solved quite 


easily for ¢, the number of ideal stages required. 


p/(a; — >) = 1 — 


t 11 stages, 


Using the values of 2; », x; and g» obtained 


above one finds ¢,=— 1-62 and R 5-8 A similar calculation gives the number of 


(Implicit in this calculation is the assumption stages required in the other section of the cascade. 


that the reflux stream to the cascade is a carrier 


t = 18 stages. 


lean phase). 
To calculate the number of stages required for McCarr and Ssrru find that R 5-5 and 


min 


the M.C.H. end of the cascade when the reflux the total number of stages to be 29 by the 
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McCabe-Theile the 


Ponchon-—Savarit construction. 


construction, and 26 by 


The differences in the answers obtained by the 


three methods are most probably due to slightly 


different use of the somewhat scattered equilib- 
rium data. 


Mr. 
this work for publication. 
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NOTATION 
component flow rate in the V’-phase 
component flow rate in the L-phase 
net upflow of component 
}-phase flow rate 
L-phase flow rate 
mole fraction of a component in the V-phase 
mole fraction of a component in the L-phase 
carrier content of the }-phase 
carrier content of the L-phase 
net upflow of carrier 


n = the 
excluding the carrier 


number of components in the mixture, 


= separation factor for a component 


fraction of adsorbent surface covered by the 
jth component 


moles of i h 


component per unit area of surface 
surface area per unit weight of adsorbent 

a constant 

Henry's law constant for the i h component 
the vapour pressure of a component 

working pressure for an absorption process 

the transformed pseudo mole fraction 
parameters for UNprERWoopbs transformation 
the saturation fraction of a mixture 

product rates 

the 
the number of stages required in a cascade 


reflux ratio 


section 


Superscripts on V, L, y, a, and u :— 


denotes ordinary molar quantity 
denotes molar quantity on a carrier-free basis ; 
no superscript denotes a pseudo—molar quantity 


for these symbols only 


Subscripts 


k 


& 


i n 
is written for = 


iel 


throughout this paper 
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Abstract— An investigation of the r quired stirrer speed and stirrer dimensions for the suspending 


of solid partie les in a low viscosits liquid has been carried out The results are given im dimen 


sionless form. By means of these empirical correlations the design for this operation can be 


carried out in many simple cases kor more complicated systems it will probably be better to 


carry out trial experiments. Rules for sealing 


the condition is chosen that no solid matter 


homogeneity of the suspension throughout 


ire discussed. As a crite non for complete suspension 
rests on the bottom of the stirred vessel the 


the whole vessel is not especially examined 


The results are compared with data from literature. Theoretical consideration for the me« hanism 


underlying the lifting of the particles from 


mechanism will be continued. 


Résumé—Des recherches ont éte effectuées sur la vitesse et les dimensions de l'agitateur exigées 


bottom are not given, but research into this 


pour la SUSPcTision cle particules solides dans un Dicquriele a basse Viscosite Les résultats en sont 
donnés comme corrélation entre des nombres sans dimensions \ DVaide de ces corrélations 
empiriques, il sera possible, dans un grand nombre de eas simples, d@etablir le plan d'un appareil 


pour cette operation. Pour les svstémes plus compliqués il vauda probablement mieux effectuer 


des essais sur model Les régles « 


eritere d'une complete na chotsi 


sur le fond du récipient agité. L*h 


ernant la 


yrandissement a Péchelle sont discut‘es. Comme 
la condition qu'aucune matiére solide ne reste 


gencite de la suspension a tout endroit dans le ré& ipient 


n'a pas fait Pobjet d'un examen spécial. Les résultats en sont compares avec les données fournies 


par la littérature Des observations the 


des particules du fond ne sont pas données 


continuces 


Zusammenfassung Zur Suspendierung ve 
Viskositét wurde eine Untersuchung tiber 
Abmessungen des Rithrers anyeste! 
Mit Hilfe dieser empirischen Bezichungen ist 


rijues sur le mécanisme qui est a la base du souléveme nt 
i 


mais les recherches sur ce mécanisme sé ront 


Festteilchen in einer Fliissigkeit mit geringer 
cli benotipte Ruhrerye schwindigkeit und 
Die Ergebnisse sind in dimensionsloser Form anyegeben 


die Vorausberechnuny dieser Operation in einigen 


einfachen Fillen méglich, Bei | plizierteren Systemen wird es wahrscheinlich besser sein 
Vorversuche anzustellen Auel le Regeln zur Massstabsvergrésserung werden diskutic rt 
Als Kriterium fiir vollstindige Suspension wurde dic Bedingiing gewihlt, dass kein Feststoff 
auf dem Boden des Riihrkessels zuriickbleibt die Homogenitiét der Suspension im ganzen 
Kessel wurde nicht besonders ge prift Die Ergenbisse werde mit Literaturwerten ve rulichen 


Theoretische Betrachtungen fiber det 


Mechanismius dles ms «der leik hee nm vom Boden 


werden nicht mityeteilt, aber die Untersu hungen tiber diesen Mechanismus sollen fortgesetzt 


werden. 


INTRODUCTION 
One of the applications of rotating stirrers is the 
suspension of a solid in a liquid. The require 
ments for the homogeneity of the suspension 
vary according to the purpose of the operation, 
Sometimes a high che gree of homogeneity Is 


required, e.g. when the stirred vessel is used as 
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a mixer from which the suspension must be fed 
continuously into a chemical reactor: however. 
in many cases, it will be suflicient if the liqnuid 
phase is well mixed, and all the particles of the 
solid phase are in suspension, so that the whole 
solid surface is well exposed to the liquid. As 


examples the batchwise dissolution of a solid in 


| 
voL 
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a liquid, or a reaction vessel with a suspended 
catalyst may be mentioned ; in these cases it ts 
often immaterial whether the upper layer contains 
only few or even no suspended particles, 

The experiments described in this paper were 
carried out in order to find the conditions for 
complete suspension of all the solid particles, 
without obtaining a quantitative description of 
the distribution of the particles over the height 
of the vessel. At the point of complete suspension, 
this distribution Is fairly homoge neous whe th the 
particles are small and light, whereas coarst and 
heavy particles remain in the lower part of the 
vessel, the upper part being a layer of clean 
liquuid, 

The investigation was restricted to fully battled 
vessels, It is not certain whether the use of 
baffles is always desirable for this operation, 
but for various reasons they ar used In many 
cases, especially with high speed stirrers. 

Previous investigations have furnished only 
very limited quantitative information regarding 
the stirring conditions required for complet: 
suspension, HirsekorRN and MILLER t dealt 
with the laminar region, where Re 12. They 
give no correlations, but only scaling-up rules. 
In the present work the Reynolds number ts 
always above 3,000, Hixson and co-workers 
[2, 3, 4] published a series of papers on the 
mixing of solids, and used the rate of dissolution 
of the solid as a measure of the intensity of the 
agitation. Most of their experiments were carried 
out in vessels without baffles, with low. stirrer 
speed, and relatively large stirrers. Wurre and 
Summerrorp et al. [5) studied the distribution 
of sand under the influence of a simple paddl 
agitator in one unbaflled tank. RAGHAVENDRA 
Rao and Munerst [6] investigated the same 
problem in small vessels of 3 1. capacity. These 
authors do not give definite general rules for 
determining the properties required in the stirrer 
system for complete suspension. 


KNEULI [7]. in his investigation of the rate of 


solution of crystals, stresses the importance of 


the point at which all particles are just in suspen- 
sion. In many cases it is most eflicient to design 
a vessel for this condition, because with an 


increase of stirrer speed above this point the 
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mass transfer from the particles to the liquid 
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increases only very slightly. KNeuLe published 
a relation indicating the power per unit volume 
required to reach and maintain this complete 
suspension. In this relation a constant appears 
the value of which is not given, and which 
depends on the type of stirrer. At the end of this 
paper KNEULe’s relation is compared with the 
results of the present investigation. 


2. EXPERIMENTAL 

The effect of rotating stirrers on the suspension 
of solids was studied in transparent vessels 
of glass and Perspex which were fitted with four 
baffles along the walls. Only visual observations 
were made, no samples being taken from the 
suspension vessel; this made the examination of 
a great number of combinations feasible. To 
ascertain whether the solids were completely 
suspended, the bottom of the vessel was observed; 
this could be done easily by illuminating the 
contents and looking through an inclined mirror 
which was placed underneath the vessel. 

At low stirrer speeds, there is a layer of solids 
at the bottom, the particles being at rest, but 
when the speed is increased the solid is partly 
entrained by the liquid. At still higher speed, 
the bottom is partly swept clear by the liquid 
flow, but there are places where some solid 
remains piled up. The position of these last 
heaps of solid depends on the character of the 
flow in the lower part of the vessel and on the 
form of the bottom. At a certain stirrer speed 
the last solids are lifted from the bottom. 

It was diflicult to determine this point exactly 
and objectively, because even at high stirrer 
speeds streaks of solid particles are always 
visible at the bottom. At the border between 
incomplete and complet suspension there are 
particles which settle temporarily at the bottom 
and remain for a short time in a fixed position 
relative to each other. When such a small pile 
remained at rest longer than 1 or 2 seconds 
before being broken up the suspension was 
judged incomplete. When no deposits remained 
on the bottom for more than Isec_ the 
suspension was considered complete. Following 


this rule the minimum stirrer speed for complete 
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suspension can be determined, by the same ob- Table 1. Dimensions of the vessels 


server, with a reproducibility of 2 to 3 per cent. 


This minimum stirrer speed was determined Vessel diameter (m) | 0-154 | 0-192 | 0-24 | 0-29 0-45 0-60 
Volume (1.) 29 55 11 19 170 


in a great number of cases, varying the vessel 


dimensions, the type of stirrer, the dimensions 


and the position of the stirrer. the amount, 


Table 2. Types and dimensions of stirrers 


the density and the particle size of the solid, 


and the density and viscosity of the liquid, 


Type of stirrer Diameter (m) 
The cy lindrical vessels we re fille d in veneral 7 7 
Paddles, D/W 2 0-08 0-16 
to a height equal to the vessel diameter, 
Paddles, D/W 0-06 OO8 O16 (0-224 
although the liquid height has only a smal Flat blade turbines | 0-06 0-08 0-12 | 0-20 
influence. The dimensions of the vessels Vaned disks 0-06 0-08 10-10 (6-12 | 0-16 | 
are listed in Table 1. In the majority of the Propellers 0-05 OOT | O115 


experiments flat bottoms were used. Two 


of the vessels of 0-192 and 0-29 m diameter, 
could be fitted with flat bottoms, dished Table 3. Density and particle size of the solids 


bottoms (radius vessel diameter) and 
conical bottoms (120°). Along the side of Solid Cy (kg/m) Particle size () 
the vessels four baffles were mounted. with 
widths 0-1 of the tank diameter. 


The experiments were performed with chloride 2160 125-150 150-250 250-350 


Sand 2600 125-150 250 350 710-850 
Sodium 


propeller stirrers with flat blades, with two- 


bladed paddle stirrers with ratios of dia- 


meter to width of DW t and 2, flat Only two solids were used, sand and sodium 
blade turbines and vaned disks. The dimen- chloride. From both solids thre« particle 
sions are given in Fig. 1 and Table 2. The size fractions were obtained by sieving (see 
stirrers were mounted centrally, with vertical Table 3). The amount of solids in the slurry 
shafts, at various distances from the bottom. was varied from 0-5 to 20 per cent by weight. 


D 
Paddle st rrery = Paddle stirrer =4 


Propeller 


Fic. 1. Five types of stirrers used in the experiments. 
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Properties of the liquids 


Carbon tetrachlorick 


Potassium carbonate 


The liquids were water. 
a solution of potassium car 
bonate in water, and a light oil. The 
physical prope rties 
listed in Table 4. 


at room te mperature are 


EXPERIMENTAL 


The following factors ce 


solid particles will be in suspension ; 


vessel diameter 


~ 


stirrer diameter 
stirrer speed 
. particle size 
. density of the solid 
density of the liquid 
kinematic viscosity 
weight of solid per weight of 
liquid times 100 
, acceleration due to gravity 
From these 10 quantities, a set of seven dimen- 
sionless combinations can be formed : 
T/D, T/C, 


Reynolds number 


geometrical ratios 


Froude number 
density ratio 


» percentage of solids 


A disadvantage of this set is that most of thes: 
dimensionless 
dependently, thus complicating the analysis 
the experimental data. It is possible, however. 


to adapt the set of dimensionless numbers to 


experimental 
each quantity that is varied in the 
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only with quantities which are constant. To 
simplify the analysis further, at first the experi- 
ments with sand and water were considered 
separately. Then there remain six quantities 
Which can be varied independently, =m 7 €& 
D.n, x and B. the other four quantities remaining 
constant. From dimensional theory it follows 
that from experiments on these quantities the 
relation between all dimensionless numbers can 
be established. except for the influence of the 
density ratio {s/C,. The dimensionless numbers 
which can be used for the analysis are in this 


Case : 


and B (1) 


It will be seen that each of these groups can 
be varied independently in the experiments, 
without changing the values of v and g. When 
an empirical relationship between these groups 
is found in this way, theoretically it is always 
possible to convert this into a relationship 
between another set of dimensionless groups 
should this be more convenient. 

From over 1.000 experiments with sand in 
water, with different types of stirrers, it was 
found that the required stirrer speed for complete 
suspension could be expressed very well as a 
product of powers of 7", D, x and B, but the 
dependency on C could not be expressed as a 


power function. 


In these experiments it was natural to 
consider the stirrer speed as the dependent 
variable, since this could be adjusted most 
easily. However, there is no fundamental 
reason for choosing this quantity as the 
dependent one. When designing an agitator 
for this operation in a given case one will 
rather choose n beforehand and consider 
D as the variable which is to be adjusted 
to obtain the desired effect. 


From a graphical analysis of the experimental 
data the following relation was found : 


n kK, B for constant (2) 


Because the influence of C could not be 
expressed in the same manner it appeared 


Liquid (kg/m") » (eP) v (m? /sec) 
Acetone 0-31 0-39 
1600 10 065 10°6 
solution lth 35 
Oil 840 0-3 ® 
i 
; 
VOL. 
as 5. 1s 
] O58 : 
termine whether the 
(111) 
, distance between stirr r and 
bottom m 
combining 
experiment 
| | 


Ta, 


to be more appropriate to use the ratio T/C 
as a constant parameter in this equation 
than the form C (g/v*)’. This could be done 
easily because the value of C could be 
adjusted in a very simple way so that the 
ratio T/C could be chosen at will, indepen- 
dent of the value of the other variables. 


From the dimensional considerations given 
above it follows that the general relations 


between the dimensionless forms will b« 
[a e/g B (3) 


Here the constant kK, can be inferred from the 
value of K, by substituting the numerical values 
of g and v for water. Rearranging and sub- 


stituting the symbol f for a hb: 
n DY Dit B 


The values of the exponents ¢ and ¢ were 
found to be constant. viz. « 0-20 and ¢ 0-13, 

The exponents a and 6 varied appreciably for 
the different stirrer types, but their sum f was 
fairly constant, lying in all cases between the 
extremes 0-78 and 0-94. From this it 
follows that the exponent of v in (4) must lie 
between 0-05 and 0-16 and the exponent of g 
between 0-42 and 0-47. These values represent 
the extremes from a great number of estimates 
of the exponents on T and D from experimental 
series in which these dimensions were varied, 
keeping vt, Band T/C constant. 

In this manner the influence of viscosity and 
gravity could be established without varving 
their values. However, for the purpose of 
determining the effect of the density ratio, it 
is necessary to examine other solid-liquid com- 
binations. When this is done, the viscosity varies 
simultaneously, thus allowing for checking the 
relation given above between » and vp. 

The experiments with other liquids were made 
only in the three smaller vessels, again with all 
types of stirrers and with all possible combinations 
of vessel and stirrer diameter. With the solids 
and liquids listed in Tables 3 and 4 it was possible 
to study eight combinations; in these the 
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viscosity varied between 0.39 1° and 
11-1 10° m?/see (0-39 and 11-1 eSt), and the 
density difference between 720 and 1810 kg ‘mm’. 
From these experiments it was found that : 


This formula, which was found independently 
from the previous one, is in very good agreement 
with the results on sand in water. The exponent 
of v is just halfway in the range which was found 
earlier and, moreover, from the value of 0-10 it 
follows that the exponent on g must be 0-45 
[see formula (4)]. Thus the variables g and 
At/t, have the same exponent, and can be 
combined. This purely empirical result is there- 
fore in excellent agreement with what would 
have been expected, Further, the value of f 
appears to be 


for constant T/C (5) 


The values of A and a depend on the type of 
stirrer and on the position of the stirrer, that is 
on the value of 7 ice In dimensionless form the 


relation becomes : 


85 


(9 AL cy” 


= K(T/DS (6) 


The left hand side of this equation can, if 
desired, be split up into more conventional 


LrOUps 
(n (fn? dD g At)’ 7) ry’? 


Here the first group is a Reynolds number 
and the second a modified Froude number. 
Finally all the rimental values were 
plotted in vraphs, one for each type of stirrer, 
with 7 Don the abscissa, and the dimensionless 
form : 
yl (yg Al ry" (4) 
on the ordinate axis. The ratio T/C was con- 
sidered as a parameter. Figures 2 to 6 show 
these curves for the various stirrers : owing to 
the numerous points these have not been 
included. In general the deviations are only 


slight, rarely higher than 10 per cent, 


| 
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Fic. 4. Correlation fi ‘ol ‘te st sing 
«a. 2. Corre lation for complete Susp nsion with dD W » ‘ orretatio reo nplete Isp nsion using i 
7 turbine with six flat blades : here Ne 6-2 

and Ne aD. 


20 ry 
| 


T 0 
O Fic. 5. Correlation for complete suspension using a 
: vaned disk ; here Ne = 4-6. 
Fic. 3. Correlation for omplete suspension with D/W 
and Ne 2-5. 


20 


Referring to Fig. 6, the line for the pro- 


pelle r shows a break which corresponds with 


10 
a change in the flow pattern. All propellers : 
were rotated in such a direction that the 8 
liquid flow was downward. In comparing a 6 


series of propellers with increasing diameter 


in a vessel of constant dimensions, that is. 4 
decreasing ratio 7/D, it is found that for : 
D < 0-45 T the solid particles move radially D 


outward along the bottom, and the last piles 


Fic. 6. Correlation for complete suspension using «a 
propeller ; here Ne = 0-5. 


of particles occur at the periphery. 
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ru. 


solids 
the 
pattern ts 


the 


into 


0-45 T 
up 


when D 
move inward and centrally 


However, 


propeller. In this case the flow 


somewhat analogous to that of a radial flow 


stirrer. There is no doubt that the shape 
of the propeller blades will influence this 
phenomenon, but probably all types of 


propellers will exhibit this sudden decrease 


in efliciency at low values of TD. A similar 


effect 


with axial turbines i.e., paddles with inclined 


was observed in some experiments 


blades. 
The influence of the shape of the bottom 
When using paddles, th 


case of a 


is small. solid is 
more easily suspended in the 
dished bottom, especially with large values 
of T/D. Equally, conical bottoms sometimes 
but 
relatively small. When propellers are used 
the different the 


practically equivalent. The homoge 


offer an advantage, only when C 1s 


shapes ot bottom are 
neity of 
the suspension is not influenced noticeably 


by the shape of the bottom. 


4. MEASUREMENTS ON Pires 
Duranp and Conpo.ios [8] investigated the 
minimum liquid velocity, necessary to prevent 
deposition of solids in horizontal pipes for 
hydraulic transport ; they give the relation : 

v,/[2 D, = F, (8) 
Here F, is a dimensionless number, which 


depends on the particle size and the concentra- 
tion of the solids, v, is the “limiting velocity 
of the 
pipe. The value of F, is plotted in their paper as 


for deposition,” and D, is the diameter 


a function of z for different values of the con- 
When the data of thes 


compared with the results of the present 


centration. ithors are 
investi 
gation the linear velocity in pipes can be re garded 
as somewhat analogous to the stirrer speed in 
an agitated vessel. Then the following conclusions 
can be drawn : 

For hydraulic transport, v, is proportional 
to (g 


stirrer speed is proportional to (g A fC) 


required 


and for suspension the 


Duranp and Conpo.ios do not mention the 
influence of the liquid viscosity. Although they 


N. 
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used only water as a liquid, and therefore did 
not vary the viscosity, the Reynolds number was 
a large In the experiments 


varied over range. 


on suspension in agitated vessels, the influence of 
the viscosity proved to be very small, n ~ v'. 
The 


carried out over a large range of particle sizes, 


hydraulic transport experiments were 
from 20 « up to 10cm. The dependence of v,, on 
function. 


200 and 


r cannot be expressed as a power 


However, in the region of 2 between 
800 » and B between 2 per cent and 15 per cent 
it can be approximated by v, ~ 2”. 

In the experiments on suspension this relation- 


ship is described by i~a 

For hydraulic transport the influence of the 
to be ~ for x 
between 200 and 800 « and B between 5 per cent 
the 


concentration appears 


and 30 per cent. In stirred vessel it was 
found that n ~ B’™, 
On the 


in these different fields agree remarkably well. 


whole the results of the experiments 


However, the influence of the overall dimensions 


of the apparatus is different in the two cases, 


but it is diflicult to make an exact compar!- 


son. When sealing up the agitated vessel, nD?’* 


must be kept constant. The peripheral speed of 
the impeller will therefore increase with D’", 
When sealing up a transport pipe the velocity 


must be increased proportional to the square 


root of the seale factor. No explanation of this 


disagreement has been found. 


Regarding the homogeneity of the mixture, 


Duranp and Conpo.ios give the following 


classification : 
with 


Particles a diameter of less than DO 


form a homogeneous mixture even when the 
movement of the liquid ts very slight. 
Particles from 50 to 200 do not form a 


perfectly uniform mixture, but when they are 
completely suspended, they are distributed over 
the whole height of the liquid ; this distribution 
approaches uniformity at higher speeds. 

Particles above 2mm move only by successive 
jumps along the bottom. They remain in the 
lower part of the liquuid, 

Particles from 0-2 to 2mm form a transition 


category. 
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The same phenomena are observed in agitated 
vessels. Particles with a diameter of 800 » do 
not reach the upper layer of the liquid, particles 
under 200 , form a fairly homogeneous suspension 
when they are completely suspended. Probably 
with complete suspension, the homogeneity in 
agitated vessels is better than in transport pipes, 
because the stirrer induces strong vertical currents, 
as compared with the flow which is mainly 
horizontal in transport pipes. 


5. ScaLING Up 


An important question is the applicability of 


the results of the small scale experiments described 
to the design of technical scale apparatus, A 


direct application will nearly always be an 


extrapolation, because the interdependency of 


the appropriate dimensionless forms has been 
investigated only in a limited region. This can 
most clearly when examining the 
(1). 

dimensionless forms 


Table 5. 


be shown 


set of forms In the experiments, these 


were between the limits 


shown in 
Table 5. 


Dimensionless forms and the limits involved 


Form Limits 


31 to 600 
15 to 600 


D 10 to 200 
0-028 to 1-5 


B 0-5 to 20 


These quantities were varied independently, 
with the exception that nearly 
T/C <7. 
within the region defined by these inequalities, 


always 


For systems 


the correlation is established fairly well. However 
in many cases for technical scale apparatus the 
value of the first three dimensionless forms will 
be appreciably higher, so that the experiments 
cannot be looked upon as truly representative. 
Theoretically, in order to apply the results to 
large vessels, it would have been necessary to 
use liquids of a much lower viscosity and, even- 
tually, solids of finer size fractions. Nevertheless, 


in order to get an idea of the rules for scaling-up, 
and to compare them with statements put for- 
ward in the relevant literature, we shall consider 
an extrapolation to larger values of these dimen- 
sionless numbers as legitimate. 

The only rule, known for the scaling-up of an 
agitator which must keep a solid in suspension is 
the rule of equal power per unit volume. This 
1937 in a 
rather general form. KNevu Le [7] used the same 


(9) in 


rule has been given by Bicue 


rule but extended its significance by relating the 
required power per unit volume to some of the 
properties of the liquid and the solid (see below). 
When considering this rule it is not clear whether 
the authors mean to say that it can be applied 
only to geometrically similar systems, or whether 
it has a more general application. It is possible 
now to compare this statement with the results 
of our experiments. 

Assuming geometrical similarity between model 
and technical apparatus, and equal liquid and 
solid properties, that is also equal particle size 
in both cases, the stirrer speed for complete 
suspension [see (6)] is related to a linear dimension 
by the formula : 


— constant. 


region : 


Further, in the turbulent 


P ~n® D°. From this it is found, by elimina- 
tion of n : 

P ~ D’*. Thus for the power per unit volume 
we find, as V~ D’: 

P/ V D 40-55 (9) 
for the same effect. 

Therefore, the required power per unit volume 
decreases when increasing the dimensions of the 
apparatus. The rule of equal power per unit 
volume is, therefore, on the safe side. 

When designing the larger apparatus it will 
often be convenient to choose one out of certain 
fixed stirrer speeds, especially when using high 
speed stirrers. It would be possible, then, to 
fix the stirrer speed beforehand, and to adjust 
the impeller diameter so as to get the required 
unit volume, thus vitiating the 


eT 
I I 


geometrical similitude. Thus it is interesting to 
consider whether the required power depends on 


the geometrical ratio 7'/D. 
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When v, x, Cs, Band T/C are fixed, our 
correlation reads : 
nD?™ — A,(T/D)*. From this it is found that : 

The power per unit volume will thus de pe nd on 
the ratio TD, unless the exponent 3a 2.45 
vanishes, which occurs for a O-S2. This ts 
true only for propellers, whereas for the paddle s 
the value of a is about 1-3, and for the turbines 
approximately 1-5. The conclusion drawn ts 
that for radial flow mp llers it is more economical 
to choose a large impeller and a low speed, but 
that for the propeller it matters little which value 
of the ratio 7T'/D is chosen, so that these are 
more suitable for use at high speeds. 

Summing up the conclusions of this section a 
confirmation is found of the validity of the 
scaling-up rule based on the power per unit 
volume, with two provisos: it is conservative 
because in larger apparatus the power can 
probably be lower, and secondly, for radial flow 
impellers it applies only when there is geometrical 
similitude, but for propellers it is more generally 
valid, in that the ratio of tank diameter to 
impeller diameter can be varied at the same time, 


6. KNeure’s Resutts 

KNeuLe [7] has undertaken an investigation 
into the rate of solution of crystals in stirred 
vessels. He concludes that it is profitabl to 
design the agitator so that all the particles ar 
just in suspension. An increase in stirrer speed 
above this point does not pay, as from this pot 
onwards the mass transfer from solid to liquid 
increases only very slowly with the intensity 
of the agitation. 

The condition for complete suspension is given 
in his paper in the form of a formula for the 
required power per unit volume as a function of 
liquid and solid properties, In the symbols used 
in this paper, his formula becomes : 

V 


os 


P= Ay (11) 


Ss 


where A, is a dimensionless constant coetlicient, 


depending only on the geometry of the agitator 


system. The value of A, is not given in KNEULE’s 


paper. 
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To compare this with our results, the following 


expression can be derived from (6) : 


os A O15 


3 


p) pow | p) 


The dependence of P on (5 and ¢, thus differs 
by a factor ACT®’, This factor is, 
however, approximately constant in the normally 
encountered region. In fact, for values of (<¢/C, 
between 1-02 and 4-0 it varies between 1-36 and 
1-70, with a minimum at C./¢, 1-43. In our 
expernments the density ratio was varied between 
1-35 and 3-2, and therefore it is impossible, on 
the basis of our experiments, to decide between 
these two « xpressions, 

The dependence of P on the particle size @ 
differs by a factor"! 
of the experiments ts also in fairly good agreement. 

The dependence of P on B is given by B’* and 
B’™ respectively. As B was varied over a 40-fold 


. which, for the limited region 


range im our experiments this discrepancy Is 
rather large. 
The ae ney ons and is not yiven by 


ULE so that a comparison cannot be made. 


NOTATION 


B weight of the solids in suspension, per 
weight of liquid, times 100 per cent 
( distance between the stirrer and the 
bottom of the vessel m 
dD stirrer diameter m 
P power input to the stirrer Nm /see 
T vessel diameter m 
volume of the vessel 
" vertical dimension of paddle stirrer m 
g acceleration due to gravity m_ sec” 
" stirrer spec d rev, sec 
s dimensionless form, from equation (7) 
’ particle size of the solid m 
kinematic viscosity of the liquid mm? /see 
ty, = density of the liquid kg/m* 
ls density of the solid 


Ne _ Newton number 


Re Reynolds number 
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Abstract—A mathematical solution to the velocity and concentration distributions of an 
incompressible, isothermal, planar laminar jet issuing from a long rectangular slot into quiescent 
surroundings of the same density is presented. The features of the solution are : concentration 
and velocity distributions remain finite back to the orifice : the orifice is finite in width « the 


change in concentration distribution, which is flat at the orifice, and the change in the velocity 


distribution from the Poiseuillian form at the orifice are taken into account 

The present solution for velocity distribution is a modification of that of Okane [8, 9] who 
solved for the velocity distribution in a viscous jet for flow from a nozzle which provides a uniform 
velocity profile at the throat. Previous analyses of the planar, laminar jet assumed the jet to 
originate from a line source. 

AnprRADE [1] reported experimental velocity distribution measurements of a laminar jet 
issuing from a nozzle under conditions similar to those specified by the theoretical studies 
of Okabe. The results of these two investigations are compared here after extension of numerical 


values of the essential parameters beyond those originally published. 


Résumé — |. auteur présente une solution mathématique de la répartition de la vitesse et de la 


concentration d'un jet plan laminaire isotherme incompressible arrivant par une fente rectangulaire 


dans un milieu immobile de méme densité. Les carac téristiques de la solution sont : Répartition de 


la vitesse et de la concentration limitée en arriére de Vorifice ; Orifice limité en largeur ; Il est 


tenu compte de la variation de répartition de la concentration, uniforme A Vorifice, et de la 


variation de répartition de la vitesse a partir de la forme de Poiseuille A lorifice. 


La solution préconisée ici pour la répartition de la vitesse est une modification de celle de 


Okabe (8, 9) qui donnait la répartition de vitesse dans un jet visqueux a travers un ajutage 


produisant un profil de vitesse uniforme a la buse. Les études précédentes dun jet plan laminaire 


supposaient celui-ci issu d'une source linéaire. Andrade (1) experimente la répartition des vitesses 


d'un jet laminaire dune tubulure dans des conditions analogues a celles spécifiées par Okabe 


dans ses études theoriques, Lauteur compare ici les résultats de ces 2 études, aprés extension 


des valeurs des paramétres principaux au dela des valeurs publiées précédemment. 


Zusammenfassung Fiir die Geschwindigkeits- und Konzentrationsverteilungen in einem 
inkompressiblen, isothermen, ebenen, laminaren Strahl, der aus cinem langen rechteckigen 


Schlitz in eine ruhende Umgebung gleicher Dichte tritt, wird die mathematische LAsung mitge- 


teilt. Die Merkmale der Losung sind : Konzentrations- und Geschwindigkeitsy erteilungen bleiben 
endlich hinter dem Spalt ; die Miindung hat eine endliche Breite ; die Anderung in der Konzen- 
trationsverteilung, die an der Miindung gleichférmig ist, und die Anderung in der Geschwindig- 


keitsverteilung vom Poissculle-Profil an der Mindung werden beriicksix htigt. 
Die mitgeteilte Losung fir die Geschwindigkeitsverteilung ist eine Modifikation der Lésung 


von Okabe (8, 9), der die Geschwindigkeitsverteilung in einem zihen Strahl berechnete. der aus 


einer Diise mit gleichfOrmiger Geschwindigkeit austrat. Friihere Untersuchungen betrafen 


den ebenen, laminaren Strahl, der aus einer Linienquelle stammte. 


Andrade (1) berichtete tiber gemessene Geschwindigkeitsverteilungen an einem laminaren 


Strahl, der aus einer Diise unter Ahnlichen Bedingungen austrat, wie sie Okabe seinen theoretischen 


Untersuchungen zugrunde legte. Die Ergebnisse dieser beiden Arbeiten sind hier zum Vergleich 
herangezogen, nachdem die hauptsichlichen Parameter tiber ihren urspriinglichen Bereich 


hinaus ermittelt wurden. 


°This research was supported by the United States Air Force under Contract No. AF 33 (038)-23976 monitored by 
the Office of Scientific Research. This work was performed by the author while Research Assistant at the James 


Forrestal Research Center. 
TPresent address; KE. 1. du Pont de Nemours and Company, Photo Products Department, Parlin, New Jersey. 
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INTRODUCTION 
Few theoretical and fewer experimental studies 
of laminar incompre ssible, isothermal jets have 
appeared in the literature. 

ScuiicutTinG [12] set up the force balance on 
a unit mass of fluid in differential equation form 
and obtained an exact solution for u (2, y) and 
v (x, y) for the circular laminar jet and an approxt- 
mate numerical solution for the planar laminar 
jet. Bickiey [4] subsequently found that the 
force balance as expressed by Scuticurine for 
the planar jet could be integrated in closed form 
and he obtained an exact solution. In_ the 
solutions by these two authors the jet issues 
at infinite linear velocity from a point or line 
source. 

ANDRADE [1] and ANDRapE and Tsien [2] 
using a liquid jet attempted to verify experi 
mentally the work of and BickLey. 
The measurements which were obtained for the 
u (a, y) velocity distribution for the planar and 
circular jets were adjusted to the theory by the 
introduction of a vertical distance, X,. as the 
effective origin of the jet. , was related to the 
orifice width and Reynolds number by a numerical 
constant. The numerical constant, in turn, 
was evaluated from the experimental data. 

OKABE [7, 8, 9] applied two integral equations 
in two alternative solutions of the force balance 
on the unit mass of fluid for circular and 
planar jets. In his solutions, the velocity 
approaches the actual discharge velocity at the 
orifice; the orifice is finite in width; and the 
velocity distribution at the orifice is flat, so 
that the orifice must be a converging-diverging 
nozzle or a conduit with incompletely developed 
laminar flow. The solutions do not apply to the 
case of discharge from a slot with parabolic 
velocity distribution at the orifice. Unlike the 
exact solutions of ScuuicuTinc and BIcKLey, 
the unknown function appears under the integral 
sign and numerical computation is required. The 
force balance was satisfied exactly at the centre- 
line and edge of the jet and an error analysis [9] 
showed the deviations of the force balance from 
zero for a numerical case. This approach was 
initially Yamapa and Oxape [16] for the case of 
the laminar wake beyond a flat plate, ; a velocity 
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distribution was assumed which satisfied the 
momentum flow equation. 

LANDAU [6] was apparently the first author to 
obtain a solution in closed form for the force 
balance of the circular laminar jet, issuing from 
a point source, which included pressure drop. 
The solution gives infinite linear velocity at the 
origin. Yarsryvev [17] obtained a more general 
solution in which LANpDAu’s solution falls out 
as a particular case, and in which the planar 
jet also is a particular solution with the origin 
treated as a line source with the velocity issuing 
therefrom at an infinite rate. 

Squire [14, 15] also obtained an exact solution 
for the case of the circular laminar jet for the 
u (a, y) and v (x, y) distributions, including 
pressure drop. This solution also states that the 
jet emerges from a poimt source at infinite 
velocity. 

Rumer [10] continued the work of Lanpau 
and Yarseyev. He pointed out that since 
LANDAU’s solution yi lded zero mass flow at the 
origin, the solution was only a first approxi- 
mation. Lanpavu’s and Yarseyvev’s functional 
forms for velocity and pressure were used but 
a second component function was added, the 
solution to which contained an arbitrary constant. 
This constant made it possible to transpose the 
origin of the jet at distances of the order of the 
size of the radius of the tube. The author stated 
that the solution has meaning as long as the tube 
is thin, ie., 7 r,. Where r, is the radius of the 
tube and r is the distance from the origin. 

Cnou [5], using BickLey’s formulation of the 
differential force balance and the equation of 
diffusion of Chapman and Cowling, derived a 
general expression of the equation of continuity, 
in differential equation form, which describes the 
concentration distribution in the planar, laminar 
gas jet. The particular differential equation for 
the case of the planar jet was set up. If this 
equation is integrated, the result apparently is 
a concentration distribution that is an odd and 
not an even function about the axis of the jet. 

Rumer [11] using the velocity distribution 
obtained by Lanpau for the circular laminar jet, 
solved the mass balance and obtained the first 
recorded instance of an exact solution to the 


, 
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concentration distribution in the circular laminar 
jet field. His solution is valid for any value of the 
Schmidt number but the jet is assumed to issue 
from a point at infinite concentration. 

Beran [3] has presented a proof that states 
that the stream function & rf (@) used in above 
mentioned solutions is incompatible with the 
Navier-Stokes 


circular laminar jet. 


complete equations for the 
Beran’s proof says that 
there is no non-trivial solution of the form 
ys rf(@) that Navier-Stokes 
equations and their boundary conditions. 


satisties the 


Shadow photographs of the flow pattern of a 
laminar jet issuing from a rectangular slot are 
presented by SINGER {13}. 


THEORY 


The ultimate objective of the present work is 
to describe the concentration field in the planar, 
laminar jet. To do this, it is necessary to solve 
balance. Since the balance 


the mass mass 


Force balance on unit mass of fluid : 


at pomt (7, y): 


at pomt (7, 0): 


Mass balance on unit mass of fluid : 
(p) 
at (7, y) 


at point (7, 


Momentum flow equation ; 


Mass flow equation : 


1 
A 3 


a 
0 


Fic. 1. 
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Basic relationships and balances pertaining to planar, laminar jet 


contains velocity terms, it is essential first to 
solve the force balance. This general approach 
was followed by Cuou [5] and Rumer {11}. 

Known physical events associated with the jet 
field are (1) the force balance on a unit mass of 
fluid, (2) the mass balance on a unit mass of 
fluid, (3) the momentum flow equation, and 
(4) the mass flow equation. The momentum flow 
equation re lates to the total force « xerted by the 
flowing fluid across any cross-section of the jet; 
the mass flow equation states that the total 
flow of the primary fluid in the jet (material 
issuing from the orifice) is constant. Figure 1 
gives a synoptic view of the events in terms of 
the theory as deve loped here. 

The present theory develops the four above 
mentioned facts into a unified system which 
relates 


velocity, concentration and 


pressure 
distributions and effective width of the jet. 
The force and mass balance Ss are used to connect 


the w-co-ordinate distance to the dependent 


VOL 
du y* P da du dy vey 
de Uy wat 
oy” AP oy dy" 
du U 
A o o 
' OP, 
a 15 U2 » | 
ve 
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variable A, which appears in both the velocity of the primary jet decreases with increasing 
and concentration distributions. The mass flow distance from the orifice, (3) the concentration 
equation serves to determine the effective width 
of the jet, and the momentum flow equation ts 
used to determine the effective field in) which Velocity distributions 
there are pressure gradients. OKABE’s approach in planar, laminar 
woater-in-water jet 
to the velocity problem is used: he considered at room.temperature 


flow from a nozzle and the present work ts 
concerned with flow from a rectangular orifice. Orifice width =0°3 mm 
The difference in ve locity distributions for these X=t3cm Re=6°4 


two geometries is shown in Fig. 2. -=- Uniform velocity 
at orifice 


Parabolic velocity 
tration distributions are assumed and written at orifice 


The general forms of the velocity and concen- 


down in Fig. 3. The stipulated velocity distribu- 
tion states that (1) at the orifice the velocity 


profile is parabolic with the jet issuing from a 


rectangular orifice, (2) the velocity decreases 


with inereasing distance from the orifice, (3) 


the velocity gradient in the y-direction is zero 
at y = 0 and at the edge of the jet, and (4) the 
derivative of the velocity gradient, d* u ‘dy*, is 


negative at y 0 and zero at the edge of the jet. 


The concentration distribution written” in 


such a way that (1) at the orifice the concen- -0°5 
Radial distance ,cm 


tration profile is uniform with the jet issuing 


from a rectangular orifice, (2) the concentration Fic 


Velocity equation ; 


Where 
normalized velocity in X-direction at point (», A) 
velocity at (», A) in X-direction 
maximum velocity at orifice 
normalized radial position, y/A, where \ = half-width of jet at X 
\ is a function of the Reynolds No. and orifice half-width, a 


normalized concentration at point (», A) 
concentration of primary fluid in the jet, at (», A). Primary fluid is the fluid that issues 
from the oritice 


density of the jet 


Fic. 3. Equations of velocity and concentration in planar, laminar jet as a function of co-ordinate position in the jet. 
Definition of terms 
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gradient in the y-direction is zero at y = 0 and at 
the edge of the jet, and (4) the derivative of the 
concentration gradient, d* c/dy*, is negative at 


y = 0 and zero at the edge of the jet. 


Velocity distributions 
in planar ,laminar 
water-in-woter jet 

>» at room temperature 


Orifice width=0O-3 mm 
Re =6°4 


X=2-3cm o 


-O5 
Radid distance ,crm 


Fic. 4 


In the present solution for the velocity and 
the jet field, the 
force and mass balances are exactly satistied at the 
The deviation of 


these balances from zero at points between the 


concentration distributions in 
centre-line and edge of the jet. 


centre-line and edge is determined for a numerical 
case. 

Also, ake 
solution to velocity distribution [7, 8, 9] 
ANDRADI 


values of A versus for OKABE’s 


are 
’s data 


extended to A 13 to verify 


shows agreement 


v aRe 


for the planar jet. Figure 
between theory and experiment and A v. 
values are listed in Table 1. 

The method applied to this hydrodynamic 
problem is similar, in a sense, to the methods of 
chemical _ kinetics. Empirically determined 
chemical rates are used as the framework for the 
the 


equations. The assumed velocity distribution is 


mechanisms which are deduced from rate 


a mechanistic one which is applied to the general 


force balance in the hydrodynamic field. This 
method of approach is general in hydrodynamic 
problems connected with fields of flow. The basic 
relationships pertaining to the present solution to 
the planar, laminar jet are illustrated in Figs. 1 
and 3. 

By taking the pressure drop term into account 
the shear stress variation starting directly at the 
for. Velocity gradients, 


orifice is accounted 


ata 
oY 
and 6 for jets issuing from a nozzle and rectangular 
slot. The 


follows. 


0, 0-1, 0-3 and 0-6 are shown in Figs. 5 


mathematical development as 


city gradient, - 98 
versus 

Jistance,7, for 

or, laminar jet with 
ow 


O3 O4 O6 OS OY 
n= ‘3 


Fic. 6. 


Force balance 
the field of the 
laminar, incompressible, isothermal jet issuing 
(Fig. 7) 


stagnant region of the same fluid at steady 


Consider two-dimensional 


from a long rectangular slot into a 
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state. The Navier-Stokes force balance, including Let 
pressure drop, in the direction is : (6) 
2 ox A 
Mr dy dy? p ww Here +, = the shear stress at the wall within the 


rectangular slot from which the jet 


Let the velocity u be expressed in terms of 


issues 
A and y, where A is some function of z and where \pP 
A, the half-breadth of the jet, is also a function No 7 at which > goes to zero; ny = f(A) 
of z. Using the general form originated by = 
Oxane [7, 8, 9] define u such that : 
© om Within the slot u (1 (7) 


Velocity and concentration profiles of the 


plane jet 2 U 
Hence, where A =a, the half- 


Ie breadth of the orifice. Making use of equation (2), 
its first derivative with r, , its 


respect to 
oa 


78 
second derivative with respect to y, and 
Profiles in duct oy” 
ft jot near equation (6), equation (1) becomes, at » = 0, 
a in pt for from 
<e 


. 


r - aRe 
Fic. 7. J (1 — + 6A(1 — 


0 
A/a and Re aU». 


v/aRe versus A is given in Table 1. 


Here 


Mass flow equation 


The total flow of fluid issuing from the orifice 


where a, 8, y are functions of A 
is constant over any cross-section of the jet : 


n=y 
eme “ | p, udy — constant (9) 
U = maximum velocity within the slot 0 
At the orifice : 


a \ 
At » 1 @ - = @ (3) p, udy pUa p, udy at any x (10) 
on a 3 
0 


and it follows that = 
Rewriting, 


+ BB + by (4) where ¢ = p,/p 


where Define ¢ such that 
A =e(1 —e)*/A 
3 

(1 — € a e(1 


Solution of equation (4) gives « 9A, 


B 


9A, 


JOL. JN 
; 
1958 
Velocity profiles i : 
vie 
(l—»)\3 3 
a+ Bo 3A 
Qa BB + 4y 3A 2. 
i= (5) + an* + By* + yn*. (12) 
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giving 


3A 


Substituting and @ from equations (2, 12) 3 47 


into equation (11) gives : 


where 


2 \ 724.9 \3 


36 


6430 


wt 
— 
— 
~ 


13 \3 1905 ” \ 1. A? 


A versus 5 values are tabulated in Table 1. 


3141 «+ 19840 
3 


Material balance 


62 
. 140 A 108 q lhe material balance over a unit cube of fluid 
83 80 in the jet is, neglecting longitudinal diffusion, 
3339 19440 \" 

(p, U) (p, Vv) d 


Aty=0, (17) 


S230 


By use of the equation of continuity, 


or oY 


and where 
equation (16) is simplified to, at » — 0 


6 
2 ) ) 
dy Fie”; Fo=1—2—A . 20 
0 n= 


At »=1 3 25 5 
‘ — 101 + 419 ” 
0 (13) d 36 
F. 25 A — 108 A? 
B 9 5 
A (14) — 1270 
y | Is 27 
3 sf. 989 — 566 — \* 1206 At 
0 
3 5 
I 31 A 136 
‘1447 12 3 
VOL 
8 
3 61 195 
140 
235 os 37 
216 648 ot 
Is 
.. 
L 3342 \® 4. 6582 oy 
108 B24 
and equation (16) becomes 
127 
oY da dy? 
52 r° 
108 9 27 0 (19) 
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dy dy Integrating with respect to x: 
Substituting @,, and , expressions into 
ou A A Yo Yo 


equation (20) as derived from equations (2, 12) utdy — | utdy l P dy l P dy (23) 
and letting the Schmidt Number of the jet equal } ; p. p 
vie 8). 
unity, D 1, yields equation (8) Let x, = 0 and #, =<. 
After the numerical relationship between 2 and A 
has been established by use of equation (8), the ‘u2dy Pdy 4 l (24) 
concentration and velocity distributions are 15 
established for any value of x. 
Where P, the pressure at a, 0. 


Momentum flow equation Find P and P, in terms of A and » making use of 


The momentum flow equation is obtained by 


multiplication of equation (1) by dy and integra- (25) 
tion between the limits y = 0 and y = A, giving: p p 
\ \ \ dx (1 6A(1 €) 

undy t dy u | d dy (21) dy o- (1 e)* aRe 
2dx, oY pls oa 

0 0 0 0 Equation (6) becomes : 
Term Term 2 Term 3 Term 4 


Consider Term 3: ‘ 


mA (I €) 


(26) 
du | p da (1 6 A(l €) 
d | | and 
oY oY | OY | 
0 yes Pa w Ne 
du of(l —e mA 
0 at y= 0 and y A dP Zo l dA (27) 
3e)+ 6A (1 €) 


\ P 


% 8 ; 0 Substitution of the expression for P, as given 
at 
+4 4 by equation (27), into equation (24), and division 


. - , of equation (24) by U* vields the following form 
Consider Term 2: 


for the momentum flow equation : 


\ \ \ \ 
dy v du uv | u dy 5 
Pe 
pU? 
. ou AY 
1 dy | dy me 
| dy 2 | | 
0 —8e) + 6A(1 
ou 
By the equation of continuity, 
so that 4 4 2{ et e } 
u dy u~ dy (I €) 
(28 
8 P ) 


a 


Therefore, equation (21) simplifies to 15 p 
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1745 


889 .. 
A* 


6 


14 
26 125,041 \4 1,500,263 


1354 
540,000 


3.600.000 


G, = 14448 


Error analysis 


An error analysis was made of the force and 
mass balances, equations (1) and (16). The error 
in each balance is expressed in terms of F, 


and F,, defined as : 


f(A, 0) 


vy WP 
A?U 


(29) 
al 


where, by the equation of continuity, 


0) 
v= A dy 
or) 


” 


Fic, 5. 


F, = = 0 (31) 


The magnitudes of F, and F,, were determined 
from several numerical examples, and the value 
of v computed in each example. The computations 
were performed on an IBM 
machine. F,, F, and v are plotted for a typical 


epe accounting 
case at various values of A in Figs. 8 and 9. F, 
and F, are of the same order of magnitude as 
the 
v increases with » and remains positive to » = 1. 
That is, this solution of the laminar jet problem 
Actually, in a 


individual force and mass flow terms and 


states that there is no inflow. 
greater part of the actual jet there is a net 
outflow which drops to zero, followed by a net 
inflow which increases then drops to zero at a 
the Thus, the 


be considered as taking 


large distance from centreline. 
present solution may 
place in that region defined by v 0 to v equals 


a maximum. The reason for this result for v 


does not change 


) to (+) 


That it does not change 


derives from the fact that 
sign. It should change in sign from ( 
at large values of ». 
sign follows from the original definition of @ of 
equation (2). 
Both F, and F, go to 
» 0, so that F, and F, maxima are some- 
Values for F, in the 


ANDRADE’S 


zero aS > and 


as 
where at 0 z< O. 


range corresponding to water-in- 


water jet are already quite small in the z experi- 
mental range. 


in ¥-coordinote direction, | 
versus 


aterc! distance from jet centerline 


9 
— 2 
5 20 
G. - A” 
4 
160 
G, AS 
81 
45 
G, = - 
128 
G vol 
4 
$125 
G 195 
324 
du AU 
| 
of terms of force and moss | = 
C balances, & and Fo, | Ast 
0d versu 
loteral distance from jet centerline a Velocity 
~ -80 a=0°015 cm ; ‘ 
Py = POIXIO’ dyn/or* | 
OT OF OS OF 05 OF O8 OD TO 
Fic. 9. 
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CONCLUSIONS 

The merit of the present solution to the con- 
centration field in the planar, laminar jet, which 
follows from and depends on the solution for the 
velocity field, lies in the fact that the velocity 
profile and nature of the outflow from the orifice 
are taken into account, including the width of 
the orifice. Other solutions to the jet problem 
have yielded zero mass flow issuing from the 
orifice, with infinite linear velocity and with the 
orifice being either a point or a line. 

OKkaBeE’s method is fruitful in that it specifies 
the nature of the orifice and outflow conditions 


and thus is amenable to experimental verification. 


Table 1. 
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Tabulation of 5 and x/a Re values 


However, the force and mass balances are 
satisfied exactly only at the centre-line and at 
the defined edge of the jet. Also, the form chosen 
for @ yields negative values of )u/dr at all 


values of y and no inflow at the outer edge. 
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Jet with no pressure drop 


A ria Re 


0-1 1-3242 O-1014 
O-2 1-8901 0-2331 
0-3 2-863 04450 
4-422 O-8026 
0-5 OST7 1-4546 
0-6 10-628 2-550 
O-7 16-209 4407 
24:31 7-836 
0-9 35-79 13-443 
10 51-76 22-65 
11 73-52 37-42 
1-2 102-8 60-62 


13 141-4 96°56 


Jet with pressure drop 


1-000 0-0 
0-1 1-1088 0-08842 
0-2 14584 0-1676 
0-3 2-108 00-2899 
3-185 04754 
0-5 4-890 00-8032 
0-6 7-541 1-3606 
O-7 11-428 2-340 
Os 16-09 3-980 
0-9 24-67 6-687 
1’) 11-07 
1-1 51-43 18-31 
1-2 70-00 29-50 
13 07-08 16-48 


half-width of rectangular jet orifice 
D coeflicient of molecular diffusion 
J = function of 
Fy sum of all terms in the elementary force balance 
rr. sum of a'l terms in the elementary mass balan e 


P pressure 


P, — ambient pressure 
Po pressure at the jet orifice 
distance from orifice 
Re = Reynolds number, aU’ /v, unless otherwise indicated 


Sch Schmidt number, »/D 
U = maximum linear velocity of jet at orifice 
u = linear jet velocity in z direction 

v = linear jet velocity in y direction 


NOTATION 


a rv co-ordinate direction of jet, parallel to jet outflow 
issuing from orifice; a () at orifice and increases 
with increasing distance from orifice 

y =y co-ordinate direction of jet, perpendicular to 
# and defining the lateral distance from the 


centre-line y 0 


Greek : 
9A 
Y 3A 
A/a 
A = half-width of jet at A 
eu 
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Ternary liquid equilibria 


M. Raga Rao. M. Ramamurty and C. Venkata Rao 


Department of Chemical Technology, Andhra University, Waltair, India 
(Received 1 October 1957) 


Abstract—Ternary liquid equilibrium data at 3 are presented for the system 
water and methyl isobutyl carbinol with formic, acetic, propionic and butyric aci 

as solutes. The tie-line data have been correlated by the method of Hand. For 

to favour the aqueous layer, whereas propionic and butyric acids are observ 

solvent layer. With acetic acid, however, the phenomenon of solutropy is obser 

10 per cent acid composition, although the distribution ratio is closer to unity at all acid 
trations. The zone of heterogeneity decreased from formic to acetic acid, and 
increased with increasing molecular weight of the acid. The behaviour of thes 
regards the area of heterogeneity and distribution relationships has been exp! 

of the relative intluences of the alkyl and carboxyl groups of the fatty acids toward 


and water constituents of the ternary system. 


Résumé—-Des résultats d’équilibres liquides ternaires sont présentés. [ls concer 
contenant de Teau et du méthyl-isobutyl carbinol en présence de solutes 
formique, acétique, propionique et butyrique. Les résultats ont été representes 
de Hand. L’acide formique favorise la couche aqueuse alors que les acides f 
sont plus solubles dans le solvent. Cependant avee lacide acétique a 10 

le phénoméne de solutropie quoique le coeflicient de distribution soit plus 
toutes les concentrations en acide. La zone dhétérogéncité décroit de 
acétique pour augmenter ensuite avee le poids moléculaire de lacide 


systemes du point de vue de la surface dhetéerogencité et des rel 


tio 
a été expliqué en vertu des influences relatives des groupements alcoyles et ¢ 


gras sur le solvant et eau du systéme ternair 


Zusammenfassung —Fiir die Systeme Wasser und Methyl-Isobutyl-Carbinol, enthalt 
Ameisen-, Essig, Proprion— oder Buttersiure als Gelostes, werden terniire Gleichgwi htsdaten be 
30°C mitgeteilt. Die Konodenwerte wurden nach der Methode von /iand korreliert. Es wurde 
gefunden, dass Ameisensiiure die wiisserige Schicht, Propion- und Buttersiiure die Losemittel- 
schicht bevorzugen. Bei Essigsiture wurde Solutropie bei ungefiihr 10°, Séuregehalt beobachtet, 
obwohl das Verteilungsverhiltnis bei allen Sdiurekonzentrationen niher an Eins lag. Die Z 
der Heterogenitit fiel von Ameisen— zu Essigsiiure ab und stieg dann mit steigendem Molekul: 
wicht der Séure. Das Verhalten dieser Systeme in Bezug auf die Heterogenitiitsilache und die 
Verteilungsbezichungen wurde erklirt in Ausdriicken der relativen Einfliisse der Alkyl— und 
Carboxylgruppen der Fettsiuren auf die LOsemittel- und Wasserbestandteile des terniren Systems, 


THE increasing importance of liquid extraction higher acids between similar water-solvent 
as a means of separating liquid components The present work is a part of the 
when other vaporization methods are impractical general programme of work under * F 

has enhanced the need for reliable liquid equili- extraction and investigations of liquid p! 
brium data. During recent years considerable equilibria” and includes the following s 
interest has arisen in the separation of fatty (1) Formic acid-water-methyl isobutyl carbu 
acids by liquid extraction. Systematic investiga- (2) Acetic acid-water-methy! isobutyl car! 
tions on the mutual solubility and equilibrium (3) Propionic acid-water-methy! isobuty! 
distribution of the lower members of fatty acids binol. 


have been carried out to predict the behaviour of (4) Butyric acid-water-methy! isobutyl carbinol. 
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Materials used 

Formic acid, Merck, 99-8 per cent pure. 

Acetic acid Merck, 100 per cent pure. 

Propionic acid Merck, 99-9 per cent pure. 
Butyric acid, Naarden, 99.8 per cent pure. 
Methyl isobutyl carbinol, the technical grade, 
supplied by Burmah-Shell Oil Co. This was 
fractionated, the fraction boiling from 131-0 to 
181-5°C. being collected and used. 


DD,” = 0-8017 g/ml nj’ = 1-4080. 


Distilled water, free from carbon dioxide, was 
always used. 
EXPERIMENTAL PROCEDURE 

In this investigation the method of Orumer 
et al. [1] was followed in the determination of the 
ternary liquid equilibria. The mutual solubility 
and distribution data were obtained at 30°C 
using a constant temperature bath of precision 

0-1°C. The solubility data were determined 
by the turbidity method which has been detailed 
in the earlier papers [2, 3}. 


Table 1. Mutual solubility and tie-line data for the system formic acid-water-methyl isobutyl carbinol 
(MJI.B.C.) at 30°C, 


Mutual solubility data 
(Weight per cent) 


Formic acid 


Water M.LB.C. 


90-0 1-3 
760 15 
0-1 2-1 
8-0 3-8 
39-0 7 
82-1 12:1 
30-8 13-8 
28-0 
26:8 19-5 
23-3 28-0 
20-5 39-3 
16-4 50-3 
140 60-4 
9-9 726 
81 82-0 
77 86-1 


Tie-line data 


(wt. per cent Solvent layer) 


(wt. per cent Water layer) 


Formic acid W ater MLL BC, Formic acid Water 


89-0 
85-1 
82-5 
79-5 
76-5 
73:5 
73-0 
65°5 
64:5 


6-0 

13-9 

21-5 

29-2 

38-2 

41-7 

48-2 

54-6 
| (Plait point) 


VOL. 
S.no. 
1 &-7 
2 22-5 
3 37-8 
+ 45-2 
5 54-0 
6 55-8 
7 55-4 
8 54-2 
53-7 
10 48-7 
11 40-2 
12 33-3 
13 25-6 
14 17-5 
15 | 9-9 
16 6-2 
| 
| 
M.LB.C. 
| 3-7 73 | 92-6 14 
| 78 84-6 | 15 
| 9-0 8:5 | 76-9 | 1-6 
11-5 9-0 68-8 20 
13-7 9-8 59-4 2-4 
| 10-5 47°58 40 
22-0 | 12-5 37-2 8-2 ; 
22-9 12-6 
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To obtain the tie-line data mixtures of the 
three components of known composition within 
the two-phase region were prepared in separating 
funnels, shaken vigorously for about 15 minutes, 
and allowed to settle in the constant temperature 
bath for three hours with intermittent shaking. 
The phases were then separated, and the acid in 
each layer was estimated by titration with 
standard alkali. The complete composition of 
the conjugate layers was then evaluated from the 
mutual solubility data. 


RESULTS AND CORRELATIONS 
Woter 
The mutual solubility and equilibrium distribu- py, 4, Ternary liquid equilibrium diagram for fatty 


tion data for the four systems are recorded in acid-water-methyl isobutyl carbinol systems at 30°C. 


Table 2. Mutual solubility and tie-line data for the system acetic acid-water-methyl isobutyl carbinol 
(M.1.B.C.) at 30°C, 
Mutual solubility data 
(Weight per cent) 


Acetic acid Water 


49 
17-0 
32-9 
34:3 
37-6 
37-7 
37-6 
36°8 
35-4 
28-8 
21-4 
18-4 

5-7 


93-8 
81-2 
68-6 
63-0 
60-8 
520 
45-7 
39-9 
33-5 
27-3 
25-1 
16-9 
12-6 
10°8 

738 


1-3 
18 
2-9 
41 
49 


Tie-line data 


(wt. per cent Solvent layer) 


Acetic acid | W ater 


7-5 88-8 
| 
9-2 80-1 
10-7 | 730 
13-0 65-0 
318 31-0 


MLB. | 


(wt. per cent Water layer) 


Acetic acid 


M.I.B.C. 


3-5 
5-8 
10°8 
17-4 
24-1 

(Plait point) 


13 
14 
17 
19 
2-6 


asia 
| frm aN 
{ 
JOL. 
S.no. 
2 
3 
4 
6 10-4 
7 16-6 
8 22-5 
29-7 
10 37-3 
11 | | 40-4 
12 54-3 
13 66-0 
14 | 70-8 
15 86-5 | 
| | | 
| 
S.no. = —— = 
| | | | 
2 | 6-0 | 92-8 | 
3 10-7 87°5 
4 16-8 | 80-7 
5 | 22-0 | 73-3 
87-2 | | 
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Tables 1-4. For comparison the saturation 60 
isotherms of all the systems are drawn in Fig. 1, 1 
and the equilibrium distribution curves are shown 


CFormic acid 
© Acetic acid 

* Propionic 
e Butyric acid 


in Fig. 2. 


The plait points of the systems were estimated 
by the method of Treypat et al. [4] and their 


compositions are recorded in Tables 1-4 under 


the tie-line data. The estimated plait points are 


also shown plotted on the respective saturation 


8 


isotherms in Fig. 1. 
In the systems studied formic acid favoured 


wt. % acid in layer 
w 
oO 


the water phase in preference to the solvent 


phase, whereas both propionic and butyric acids a 


preterre d the solvent lay r. case ot acetic 
wt. % acid in woter layer 


acid, however, the distribution ratio ts closer to 


unity, the acid favouring the solvent layer in Fic. 2. Equilibrium distribution diagram. 


Table 8. Mutual solubility and tie-line data for the system propionic acid-water- methyl isobutyl carbinol 
(MJI.B.C.) at 30°C, 


Mutual solubility data 
(Weight per cent) 


propionu acid Water MUTI 


10-4 
21-1 


Tie-line data 


(wt. per cent Solvent layer) (wt. per cent Water layer) 


Propionic acid Water MLB. Propionic acid Water M.LB.C, 
8-1 7-4 84-5 2% 96-4 13 
2 14-8 8-4 76-8 49 93-7 14 
3 21-4 0-6 69-0 78 90-7 15 
‘ 28-2 110 60-8 11-3 86-8 1-9 
36-3 49-7 14-0 | (Plait point) 


VOL, 
= 
| 

2 763 26 
3 26-9 
32-5 581 
5 45-6 17-6 
6 36-8 22-2 
7 42:1 32-1 25-8 
42-2 27-2 30-6 
41-8 22-3 35-9 
10 19-5 
1! 35-3 14-0 
lz 28-9 10-6 60°5 
13 21-3 8-7 70-0 

12-1 73 80-6 

| 
268 


JOL. 


Ternary liquid equilibria 


low acid concentrations, and the aqueous layer decreased from formic to acetic acid. and ther 
in high concentrations, thus giving rise to after it increased with increasing molecular 


“ Solutropy.” The area of heterogeneity weight of the fatty acid 


Table 4. Mutual solubility and tie-line data for the system butyric acid-water—methyl isobutyl carbinol 
(MJ.B.C.) at 30°C, 


Mutual solubility data 
(Weight per cent) 


Butyric acid Water M.1.B.C. 


13-1 85-8 11 
1-2 


2:1 


8 
2 
10 58-5 
52-6 12-9 34-5 
S-4 60-3 
23-0 76 69-4 
17-7 7-4 749 
11-6 7-3 


Tie-line data 


(wt. per cent Solvent layer) (wt. per cent Water layer) 


Butyric acid Water Butyric acid Water 


6-9 83-4 98-0 1-2 
2 ™ 7 760 1-7 97-0 1-3 
3 24-9 76 67-5 2-6 96-05 1-35 
4 32-4 8-3 04-7 1-4 
5 43-6 10-4 46-0 92-05 
6 510 12:5 36-5 73 1-2 1-5 
45-0 51-6 3-4 (Plait point) 


Table 5 


Constants of hand correlation 


System 


Formic acid—water—M.1.B.C. 0-25 0-62 


Acetic acid—water—M.1.B.C, 1-05 1-00 
Propionic acid—water-M.1.B.C, 3-55 1-00 


Butyric acid—water—M.1.B.C, 


1 
2 
3 34-0 63-9 = 
+ 50-4 3-6 
5 55-9 37-9 6-2 
S.no,. - 
K n 
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Tables 1-4. For comparison the saturation 60 


isotherms of all the systems are drawn in Fig. 1, 


and the equilibrium distribution curves are shown 


oFormic acid 

| © Acetic acid 
* Propionic acid™ 
Butyric acid 


in Fig. 2. 
The plait points of the systems were estimated 
by the method of Treypat et al. [4] and their 


> 


compositions are recorded in Tables 1-4 under 
the tie-line data. The estimated plait points are 


also shown plotted on the respective saturation 


8 


isotherms in Fig. 1. 


In the systems studied formic acid favoured 


wt. %acid in M.1.B.C. layer 
w 


the water phase in preference to the solvent 


phase, whereas both propionic and butyric acids 
0 20 30 40 50 60 


preterre d the solvent lave r. In the case of acetic 
wt. % acid in water layer 


acid, however, the distribution ratio ts closer to 


unity, the acid favouring the solvent layer in Fic, 2. Equilibrium distribution diagram. 


Table 3. Mutual solubility and tie-line data for the system propionic acid water- methyl isobutyl carbinol 
(MJI.B.C.) at 30°C. 
Mutual solubility data 

(Weight per cent) 


propiont acid Water MLB. 


10-4 
21-1 


Tie-line data 


(wt. per cent Solvent layer) (wt. per cent Water layer) 
_ — _ - 
Propionic acid W ater B.C. Propionic acid W ater M.L.B.C, 
| 7-4 84:5 2-3 
2 14-8 8-4 76-8 49 93-7 1 
3 214 0-6 69-0 78 
28-2 110 60-8 113 86-8 


36-3 49-7 140 (Plait point) 


VOL 
1 14 
2 76-3 246 
3 26-9 
32-5 
5 i764 
6 4140) 36-8 22-2 
7 42:1 32:1 25-8 
s 42-2 27-2 30-6 
41-8 22-3 35-9 
10 19-5 10-3 
1! 35-3 140 
lz 28-9 10-6 60-5 
13 21-3 &-7 
12:1 73 80-6 
: 
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low acid concentrations, and the aqueous layer decreased from formic to acetic acid, and ther 
in high concentrations, thus giving rise to after it increased with increasing molecular 


“ Solutropy.” The area of heterogeneity weight of the fatty acid 


Table 4. Mutual solubility and tie-line data for the system butyric acid-water-methyl isobutyl carbinol 
(MJ.B.C.) at 30°C, 
Mutual solubility data 

(Weight per cent) 


S.no. Butyric acid Water M.1.B.C. 


13-1 
23-0 
340 
wo 


59 37-0 6-2 


Ge 


59-8 30-4 
7 61-1 25-5 13-4 
61-2 21-5 17-3 
60-2 18-7 21:1 
10 58-5 25-0 
52-6 12-9 34-5 
12 43°5 
31-3 8-4 60-3 
23-40 76 69-4 
74 749 
11-6 7-3 81-1 


Tie-line data 


(wt. per cent Solvent layer) (wt. per cent Water layer) 


Butyric acid Water M.L.B.C. Butyric acid Water 


6-9 
17” 7 760 1-7 97-4) 1-3 
76 67-5 26 96-05 1-35 
3-9 04-7 1-4 
43-6 10-4 10 5-6 92-05 145 
6 510 36-5 73 91-2 15 
45-0 51-6 34 (Plait point) 


Table 


Constants of hand correlation 


System 


Formic acid—water-M.1.B.C. 5 
Acetic acid-—water—M.1.B.C. 1-05 1-00 
Propionic 3-55 14) 


Butyric acid—water—M.1.B.C, 


1-1 
2 75-8 1-2 
3 63-9 2-1 3 

4 50-4 3-6 

5 

VOL. 
105R 
S.no,. 
K n 
14-80 1-00 
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The tie-line data are examined by the correla- 
tion method of Hanp [5] and was fitted satis- 
factorily by the following equation : 


Xx 


cs 


X 


© Formic acid 

© Acetic acid 
* Propionic acid _| 
Butyric acid 


16 00 


3-6 
Xow / Xww 


Fic. 3. Hand correlation, 
The constants K and n were evaluated for the 
Except 


0-25, 


four systems and are listed in Table 5. 
0-62 and A 
the value of n is unity for the other acid systems, 


in the A 


value from formic acid to the higher fatty acids. 


with formic acid, where » 


and further there is a gradual increas« 


Disct SSION 


An examination of the saturation isotherms 
and the distribution curves for the 
systems studied that the 


heterogeneity is larger, when the solute greatly 


four acid 


indicates zone of 


prefers one of the non-consolute liquids, and ts a 


minimum when the distribution ratio is close to 
unity. The behaviour of these systems can be 
explained as being due to the relative influences 
of the alkyl and carboxy! groups in the fatty acid 


on the other two components of the system. 


With formic acid the effect of the carboxylic 
acid group is predominant, and so the immisci- 
bility zone is increased as the acid favours the 
water phase in preference to the solvent phase. 
For higher members in the homologous series of 
fatty 
group is gradually reduced, with the result that 


acids the effect of the carboxylic acid 


when the solute is acetic acid both the methyl 
and carboxyl groups influence the behaviour of 
the acid towards water and the solvent to nearly 
the same degree. Consequently the immiscibility 
zone is reduced as the distribution ratio closely 
approximates to unity. The increase in the area 
of heterogeneity and increased distribution of the 
acid in favour of the solvent phase from acetic 
acid to propionic acid, and from the latter to 
butyric acid, may likewise be attributed to the 
increasing influence of the alkyl group in the 
fatty acid towards water and the organic solvent. 
This is in confirmity with our earlier observations 
on fatty 
fatty 

It can therefore be concluded that for a given 
the 
increasing molecular 


acid water-ester systems [6, 7] and 
acid-water-hydrocarbon systems [8]. 

combination, as solute is 
the 


weight in the homologous series of fatty acids the 


water. solvent 
changed in order of 
acid increasingly favours the solvent phase, the 
area of heterogeneity is larger as the solute 
the 
liquids, and is a minimum when the distribution 


increasingly prefers one of non- consolute 


ratio is closer to unity. 


NOTATION 
solute 
solvent 
water 
weight fraction of the solute in solvent layer 
weight fraction of the solute in water layer 
weight fraction of solvent in solvent laver 
weight fraction of water in water layer 
constants in Hand equation 
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Mass transfer from spheres in various regular packings to a 
flowing fluid 


D. Tuorenes, Jr.* and H. Kramers 
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Received 18 October 1957 


Abstract——Mass transfer coefficients have been determined for fluids flowing around single 
spheres, which were part of a regular arrangement of similar, but inactive spheres. Eight different 


kinds of packing were used. The experiments were carried out with solid spheres dissolving 
into a stream of water, and with porous spheres soaked with a liquid, which evaporated into a 


stream of gas. The results are tentatively explained on the assumption that the total mass 


transfer from the sphere surface is the sum of three different contributions : one for laminar 


convective transfer, one for turbulent convection transfer and one for diffusion in the stagnant 


regions near the contact points with adjacent spheres. The latter contribution would only be 


important for gas flow at not too high Reynolds numbers (Re < 500). 


Résumé —Des coefficients de transfert de masse ont été déterminés pour des fluides s’écoulant 


autour de sphéres uniques, qui faisaient partie d'un remplissage régulier de sphéres identiques, 


mais inactives. Huit différents types de garnissage ont été utilisés. Les expériences ont été 


réalisées avec des spheres solides se dissolvant dans un courant d'eau, et avec des spheres poreuses 


mouillées avec un liquide qui s'évaporait dans un courant gazeux. Les résultats ont été expliqués 


provisoirement A partir de Thypothése basée sur le fait que le transfert de masse global de la 


surface de la sphére est la somme de trois différentes contributions : transfert par convection 


laminaire, transfert par convection turbulente et transfert par diffusion dans la courbe immobile 


se trouvant prés des points de contact avee les sphéres adjacentes. La derniére contribution 


serait seulement importante pour écoulement de gaz a nombre de Reynolds peu élevé (Re < 500). 


Zusammenfassung—Fiir die UmstrOmung von Einzelkugeln, die Bestandteile eines regel- 
miissigen Haufwerks von dliniichen, aber inaktiven Kugeln waren, wurden Stoffiibergangszahlen 


bestimmt. Hierbei wurden acht verschiedene Packungsarten angewendet. Die Versuche wurden 


durchgefiihrt mit festen Kugeln, die sich in einem Wasserstrom lésten, und mit porésen Kugeln 


vollgesogen mit einer Fliissigkeit, die in einen Gasstrom verdampfte. Die Ergebnisse wurden 


versuchsweise unter der Annahme erklirt, dass die gesamte Stoffiibertragung von der Kugelo- 
berfliiche die Summe dreier verschiedener Beitriige ist : die laminare konvektive Cbertragung, die 
turbulente konvektive Ubertragung und die Diffusion in den stagnierenden Bereichen in der 
Nihe der Bertihrungspunkte mit den Nachbarkugeln. Der letzte Beitrag ist nur fiir Gasstrémung 
500) von Bedeutung. 


bei nicht zu hohen Reynolds-Zahlen (Re - 


1 INTRODUCTION beds has been a subject of many investigations. 


A LARGE number of chemical processes, especially In Fic. 1 results of previous mass _ transfer 


catalytic processes, are carried out continuously measurements have been collected as a plot of 
in “packed bed” reactors. As in any other Sh/Sc'/* vs. Ret. These results were obtained 
heterogeneous process, the rate of conversion is for purely physical mass transfer, viz. dissolution 


not only determined by reaction kinetics, but and evaporation, in packed beds of known poro- 
also by the mass transfer to and from the phase sity and specific surface area, consisting of 
particles that did not differ greatly from spheres. 


boundary. Therefore mass transfer in packed 


*Present address : Central Laboratory of Staatmijnen, Geleen, Netherlands, 
+See notation at the end of this paper. 
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PHOENES, 


10° 10 


10 


Mass transfer measurements im pr ked beds, Here 1 repre sents isobutyl alcohol water 16), 2 n-butyl 


iter 16), 3 succinic acid-n-butvl hol [3], 4 salievlie acid benzene [3], 5 succinic acid — acetone [3], 
i eld — givcerine-water mixture (10), 7 benzoie acid — water [10], 8 naphthol water [11], 9 water — air 
10 water air (cylinders) [5], 11 water [21], 12 various liquid-gas systems [7], 13 hydrogen 


peroxide — helium [14], and 14 naphthalene — air (pellets) [1]. 


ntal results can be summarized number as a variable. He found a much larger 


the following rough correlation spre ad. 
en Sh, Re and Sc: The lines in Fig. 1 represent the “* best lines " 
Sih 0.9 Rel/? Se (1) drawn through the points of each series of 
measurements ; the points themselves have been 
: K 5,000 nd I a 10,000) omitted in order to get a clearer picture. From 
ition of 25 per cent the data given by the various authors it follows, 
s ith have found that Sh for given however, that the experimental error and the 
Ke depended on the particle size. spread in each series of measurements are . 
( or t interpretation of their generally much smaller than the spread between 
mensionk variabl ild have the various “ best lines.” This indicates that the 
wcount, This means that either discrepancies between the various investigations 
of different sizes were not exactly are to a large extent of a systematic nature. 
pe, O1 hat an uneven flow distri For a discussion of the possible systematic 
d mixing effects in the bed affected th errors, the usual method for calculating the mass 
n of these measurements. These — transfer coeflicient from the measurements referred 


famed by Horr [8], Resxick and to will be briefly outlined. The differential 


\\ 13) and part of the results of Bar-ILan equation for mass transfer in a packed bed can 


have not been included in be written as follows :— 


vo F dC k(C* C)S F dz (2) 


GAMSON [4], who in an earlier review collected 


ewer measurements, found a smaller spread than Here, a uniform flow velocity over the entire 


t mentioned for equation (1). ExGun [2] cross-section of the bed has been assumed, and 


plotted the quantity J = 6Sh Re against Re, transport of matter by diffusion and mixing in 


nd in this way did not consider the Schmidt the direction of flow has been neglected with 


10° 
s 
10 
3 
BAA 
Ca 
10 
VOL. 
| 
Fig. 1. 
||| 
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respect to the convective transport. With the 
additional assumptions that vp, S and C* do 
not vary in the direction of flow, equation (2) 
can be integrated to give : 


(3) 


The value of k, derived in this way from average 
inlet and outlet concentrations (C, and (, 
respectively) can be regarded as a mean effective 
mass transfer coefficient for the entire bed. 
For purposes of correlation, k and vy, are combined 
with the appropriate physical properties of the 
fluid system and with a characteristic linear 
dimension, yielding relationships of the type of 
equation (1). 

In this entire procedure the following factors 
may cause systematic deviations in the final 


result of a measurement. 


(a) Bed inhomogene ity 


A bed consisting of particles, which are not 
very small compared to the dimensions of the 
tube or the container that holds it, cannot be 
considered to be randomly packed. The sections 
where the fluid enters and leaves the bed have a 
structure which is different from the main body. 
The region of a bed near the walls has a porosity 
which is larger than the average porosity of the 
bed. Accordingly, the flow velocity near the 
walls may be considerably higher than the average 
velocity ; a maximum velocity ratio of 2 has 
been observed by Scuwartz and SMITH [15] for 
a ratio of 12 between the tube and particle 
diameters. Thus in practice one of the conditions 
for the validity of equation (2) ts only approxi- 
mately fulfilled. As a consequence, a k-value 
caleulated from equation (3) may differ appreci- 
ably from the actual mean value for & in the bulk 
of the bed, in particular in short beds where 
transverse mixing has insuflicient opportunity to 
level out radial concentration gradients. 


(b) diffusion 

The value of the axial (eddy) diffusivity in flow 
through packed beds, as recently published by 
McHenry and Wineto [12], is such that for 
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high axial concentration gradients (i.e. when 
saturation of the fluid occurs over a short distance) 


its contribution to the total transport in axial 


direction may not be neglected. 


(c) Bed height 

It can be shown that, for mass transfer between 
a packed bed and a gas flow, saturation is nearly 
complete when the gas has passed a number 
of particles of the order of 10. Since for the 
application of equation (3) the value of (C* — C,) 
must be measurable with sufficient accuracy, very 
short beds have to be used. In that case the wall 
effect and the axial diffusion reduce equation (2) 
to a definition rather than an adequate descrip- 
tion of the mass transfer process*. In our opinion 
this is the main reason for the large discrepancies 
between various results for mass transfer of gases 
in the literature. 

*\ method for eliminating the complications arising 
from (a), (b) and (c) for mass transfer to gases is the use 
of a dilution method, The recent paper by Bar-ILan 
and Resnick [1] on the rate of evaporation of naphthalene 
pellets dispersed in a bed of relatively many inert pellets 


provides an interesting example of this method. 


(d) Ewperimental errors 

Out of the many pitfalls for the experimenter, 
only one difficulty is mentioned here. In deter- 
mining the mass transfer rate of evaporating 
substances it is generally assumed that the 
temperature of the surface is near the wet bulb 
temperature. It should be measured with great 
care because of the strong dependence of the 
vapour pressure on temperature. However, this 
measurement may involve serious errors due to 
heat flow through the measuring device, generally 


a thermocouple. 


(e) Phys ical properties 

Of these. the diffusion coefficients of the 
evaporated or dissolved material in the fluid 
medium is often not known with sufficient 
accuracy. Therefore a systematic error may 
easily be introduced by representing experimental 
results as a relationship between the usual 


dimensionless groups. 


v 
k= 
SL | C2 | 
OL. | 
1958 | 
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2. Tare PrINcIPLE oF THE MEASUREMENTS 


It follows from the above considerations that 
one must distinguish between the mass transfer 
coeflicient for one particle and the effective mean 
mass transfer coeflicient for a bed of particles. 
The latter is of practical interest, but the former 
is a better measure of what actually occurs 
between the particle surface and the fluid. The 
effects of other phenomena such as non-uniform 
velocity profiles and mixing, which influence the 
total mass transfer rate, can be studied separately 
[20]. 

The subject of this investigation is the rate of 
mass transfer between a fluid and the surface 
of one particle in a packed bed. If such a bed 
were a random packing of particles, statistically 
large numbers of measurements would have to be 
carried out with one particle in many different 
geometrical configurations. In order to avoid 
this complication only regular packings of spheres 
were considered. Because the mass transfer rate 
was measured for one central sphere in regularly 
packed surroundings, the volume of the packing 
did not need to be very large in comparison with 
the volume of one sphere. Therefore the packings 
that were used consisted only of the central 
sphere and a sufficient number of other spheres 
and parts of spheres in order to make a represen- 
tative packing and to fill up the prismatic tubes 
in which the packings were mounted. 

The experiments were carried out with solid 
spheres dissolving into a stream of water, and 
with porous spheres soaked with liquid, which 
evaporated into a stream of gas. In this way a 
wide range of Schmidt numbers could be covered. 
In all cases the spheres surrounding the central 
sphere were of inert, that is, non-dissolving or 


non-evaporating, material, 


3. EXPERIMENTAL DeratLs 
The packings 
The packings of spheres that were used have 
been listed in Table 1, together with some of their 
geometrical properties. Figure 2 shows a number 
of cross-sections through the packings and the 


way in which they were arranged in the prismatic 


tubes. 
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Table 1. Regular packings used 


Symmetry axis 
Packing | Geometrical form | Porosity | parallel with 
direction of flow 


la 4 
Dense cubic 0-260 

lb 3 

2 Dense hexagonal 0-260 3 

3 Orthorhombic 6 

4a 4 
Cubic 0-476 

4b 3 

5 Body-centred cubie 0-320 4 

6 | O48 4 


The packings nos, la and 1b are identical and 
differ only in orientation ; the same is true for 
ta and 4b. 

Packing No. 5, the body-centred cubic lattice, 
is in structure rather similar to a random packing 
of spheres. In a horizontal layer the spheres do 
not touch one another, while each sphere of the 
next layer lies in a “ pit” formed by four spheres. 

Packing No. 6 is not a regular packing. It 
contains 9 spheres in a_ body-centred cubic 
arrangement, without other spheres filling up the 
tube. 


Experiments with evaporating liquids 

The apparatus is schematically represented in 
Fig. 3. Evaporation occurred at the surface of a 
porous sphere made of filter paper fibres, which 
were hardened by means of a synthetic glue. 
It was connected by means of a capillary tube 
with a liquid reservoir, a test tube; the liquid 
level in the reservoir was always adjusted so 
that the sphere surface was completely wetted, 
without drops falling from the sphere. If the 
adjustment was good (roughly at the same height 
as the lowest point of the sphere), a change in 
level of a few millimetres up or down did not 
influence the rate of evaporation, under fixed 
conditions, 
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Fic. 2. Cross-sections through models of regular sphere packings. Each cross-section is made through the centres 
of the spheres in one layer. Each model consisted of three layers, except No. 4b, which consisted of five layers. 


The rate of evaporation from the sphere was 
determined by measuring the rate of flow of 
liquid from the reservoir to the sphere. To this 
end a small air bubble was inserted in the capillary 
which had a calibrated cross-section, and the 


rate of travel of the bubble was measured. The 
air bubbles were collected in the vertical side 
tube from which they were removed from time 
to time. 

Temperatures were measured by means of 
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Fic. 3. The essentials of the apparatus for gas flow, 


thermocouples (the wire diameter being 0-05 mm) 
at two points at the surface of the sphere. The 
temperatures of the air stream above and under 
the packing were measured with mercury thermo- 
meters. 

The liquids used for evaporation wer p-xylene 
and n-propyl alcohol. Several gases were used as 
a fluid: air, carbon dioxide, hydrogen and 
Freon F-12 (CCI1,F,). The Schmidt numbers 
of the chosen systems lie between 0-84 and 2-06. 

Mass transfer rates were determined over a 
wide range of Re numbers. The velocity of the 
gas flow was measured by means of a series of 
three rotameters, calibrated for air and hvdrogen. 

The diffusion coeflicients were calculated using 
the formula of Wike and Ler [22], and have 
been collected in Table 2. The Vapour pressures 
of the liquids were taken from the work of 
Stuckey and Sayior [18] and of [19). 


Table 2. Diffusii ities and Schmidt numbers used 


D at 20 ¢ Se 

Systems (10°° see!) werec) 
p-Xylene, air 7:35 2-406 
p-Xylene, CO, 6-43 1-24 
p-Xylene, Hy, 30-8 2-8* 
p-Xylene, Cl, Fy 2-93 
n-Propyl ale., air 10-2 1-14 
Benzoic acid, water 000078 STOO-450 


*Average between surface and bulk gas conditions 
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Experiments with dissolving solids 

The apparatus was essentially the same as 
that used for the experiments with evaporating 
liquids. The capillary was omitted and the porous 
sphere was replaced by a massive sphere of 
benzoic acid. These spheres were dried and 
weighed before and after each experiment, which 
consisted of an exposure of the sphere to a con- 
stant flow of water for a known period. 

Because of the negligible temperature effect only 
the temperature of the outflowing water was 
measured, The Schmidt number was varied from 
$50 to 3,700 by working at different temperatures. 

The diffusion coefficients of benzoic acid were 
taken from the work of Suerwoop and Liwrox 
[17] (see Table 2), and the solubilities from 
SEIDELL [16]. 


Discussion or tue RESULTS 


The results of the mass transfer measurements 
have been collected in Fig. 4 to 13. In order to 
make them clearer only the “ best lines” have 
been drawn. Each line has been drawn through 
at least 7 points, the average number of points 
per line being 20. The mean deviation in all 
series was less than 5 per cent except in the 
series with benzoic acid dissolving in water in 
packing number 4b; here a mean deviation of 
13 per cent occurred due to poor reproducibility 
with this complicated arrangement*, 

In Figs. 4 to 7 the data have been combined 
for given systems, the numbers referring to the 
packings ; in Figs. 8 to 13 for given packings ; 
here PX means p-xylene, PA n-propyl alcohol 
and BZ benzoic acid. From the latter graphs it 
can be seen that the influence of Re on Sh/Se' 
still depends on the value of Se; the a priori 
and very widespread assumption that Sh is 
proportional to Sc'/* at constant value of Re is 
only a rough approximation, particularly for 
gases at small Re numbers. That different 
packings give different results is only natural, 
as no packing parameters have been, nor could 
be considered as independent variables. The 


definitions of Sh and Re have only been chosen 


*The complete data have been published elsewhere 
[20] and are available on request. 


VOL 
~ 
195 


Mass transfer from spheres in various regular packings to a flowing fluid 


Fic. 4. 


in such a way that somewhat comparable results 
might be expected. They are based upon the 
hydraulic diameter ed (1 «) and the average 


axial velocity of flow in the bed vy ‘«. 


Liquids 


The results of the mass transfer measurements 
with benzoic acid spheres, partially dissolving 
in a stream of water (BZ-H,O), are the least 
was 


Satisfactory agreement 


plotting Sh/Sc'/* vs. Re, if Sc was 


complicated. 
obtained by 


10° 


Benzoic acid 


water. 


varied by a factor of 8. Thus for every packing 


the points form one single line ; these lines are 


presented in Fig. 4. 

The packing No. la shows peculiar results : 
1000 Re < 2000 and a 
comparatively steep slope for Re < 1000. These 


irregularities were reproducible but could not 


an irregular turn at 


be explained. Packing 1b, which has the same 


structure and only differs in orientation, did not 
show this peculiarity. The other lines were all 


smooth and roughly parallel. As the slopes 


Fic, 5, p-Xylene — air. 
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Fic. 6. p-Nylene 


increase with increasing Re one might express 
an average value of Sh as the sum of a laminar 
and a turbulent term both containing the factor 
Se' 


h 


Thus for packing number 5 : 
(1-26 Re'/* Se'/*) + (0-09 Se'/®) (4) 


These two the 
basis of the consideration that the flow through 


terms have been chosen on 


a packed bed is never completely laminar 1 
completely turbulent. Certain of t 


ior 


parts he 


carbon dioxide. 


surface of the particles may be exposed to 


turbulence, while in narrow and in converging 


With 


increasing Re there is an increase in the relative 


channels the flow may still be laminar. 


importance of turbulence. 


Gases 


The experiments with gas flowing through 
the bed show much more complicated results. 


In this case the value of Se was varied by using 


Fic. 7. p-Xylene — hydrogen. 
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10 


Fic. 8. Packing No. la. 


different combinations of gases and evaporating 
liquids. Contrary to the experiments with dissolv- 
ing solids the data obtained with gases did not 
form one single line for each packing. For each 


system a line was obtained which 


differed in slope from the other lines. To avoid 


separate 


confusion different graphs have been drawn for 
each packing (Fig. 8 to 13). In these graphs the 
lines from the experiments with liquid flow 
(BZ-H,O) are shown for comparison. 

The results obtained using the packing with 


h 
10°) 


| 
| 
| 
| 
| 


the highest fraction of voids (No. 6) are the least 
the fall together 
For the other packings the lines fall 


complicated : lines roughly 
(Fig. 13). 
further apart. The largest spread between the 
various lines occurs with packing la (Fig. 8). 
For lower values of Re higher values of Sh /Sc'/* 
are obtained in gas than in liquid flow. 

In comparing mass transfer in liquids with 
mass transfer in gases there is one important 
difference, the value of Sc. The fact that Sc is 


so much lower for gases than for liquids means 


Fic. 9. Packing No. 1b. 
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Fic. 10. 


that the relative of diffusion 


compared to convection ts much vreater In gases, 


importance as 


The proportionality of Sh and Se' roughly 
this fact 
diffusion occur simultaneously. But in a packed 


accounts for when convection and 


bed only a part of the particle surface is « x posed 
to 
considerable part of the surface of each particle 


simultaneous diffusion and convection. <A 


is shielded by other particles from intensive 


convection. In flowing liquids this part of the 


2 


10 


Fic. 11. 


Packing No. 4a, 


In 


flowing gases, however, the diffusion coefficient 


surface hardly contributes to mass transfer. 


is so much higher that even on the shielded part 
of the surface a certain amount of mass transfer 


takes 


stagnant 


place diffusion through an almost 
fluid. to the total 


average value of Sh for the whok particle would 


by 
This contribution 


be ck arly inal pence nt of Sc. Thus, in addition 
to the laminar and turbulent terms of the kind 


given in equation (4), a third term may be added 


Packing No. 4b. 
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which is approximately a constant. This explains 
why at low values of Re the measured values of 
Sh/Sc'/® 


accordingly lowest for liquids. 


are lower for higher values of Sc. and 


Re 


highest 


the 
for 


high values of 
Sh are 


liquids. This can be explained by assuming that 


However, for very 


measured values of 
the turbulent term contains Se with an « xponent 
higher than 1.3. 

The led to the 


following relationship for mass transfer to gases 


above considerations have 


and liquids in packing No. 5: 


(1-26 Re! ") 


laminar 


Sh (0-054 Sc” 


(turbulent 

(5) 
(stagnant 

The small positive dependence on Re of the 
third the 


of convection into the 


term indicates Increasing penetration 
quiet corners between 
the particles. Table 3 gives an illustration of the 
relative importance of the three terms of equation 
(5) for liquids and gases. The third term accounts 
for a considerable part of the mass transfer in 
vases, but only for a very small part of the mass 
transfer in liquids. For high values of Se equation 
(5) is practically identical with equation (4). 
As to the the the 


packing, it has already been suggested that for 


influence of structure of 


denser packings a larger part of the surface is 


10° 


Table 3. Sh calculated from equation (5) for 


Re 107 


Se 1-26 0-054 Re®® 0-83 Sh Sh 
Scl/8 

10-0 | 10-0 

76 76 


shielded from convection. Figs. 4 to 7 show that 


for these packings the measured values of Sh (Sc'/* 


are accordingly lowest. However there may be 


still another cause for this. Due to the definition 
of Sh and Re employed, packings are compared 
diameter. A denser 


according to hydraulic 


packing has, per unit of hydraulic diameter 
length, fewer particles than a less dense packing, 
and therefore the relative int nsity of convection 
IS smaller. 

The packings 3 and 4a are of particular interest. 
In them every sphere lies on top of a sphere in an 
adjacent layer. At the upstream and downstream 
contact points each particle is very effectively 
shielded ; this results in a rather large contribu- 
tion of the third term (compare equation (5)), 
also at high values of Re, as can be seen from 
Fig. 10. 

Finally a very rough average of all 438 mass 


Fic. 12. Packing No. 5. 
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mass transfer coefficient 


Sc = Schmidt number, Sc v/D 
Sh = Sherwood number, Sh = «kd/(1 — «) yD — 


divided by 
vy = kinematic viscosity 
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The temperature produced by heat of reaction in the interior of 
porous particles 


Cuartes D. Prater 


Socony Mobil Oil Company, Inc. Research and Development Laboratory Paulsboro, New Jersey, U.S.A 


( Received 9 September 1957) 


Abstract——Heat of reaction will cause a temperature difference to occur between the inside 


and boundary of a porous particle during an internal reaction. An equation is derived relating 


this temperature difference to the reactant concentration difference, and it is shown to be 


independent of the kinetics of the reaction and of the particle geometry. Examples are given of 


the applic ation of this equation to 


lohexane 


dehydrogenation and catalyst regeneration. 


Résumé—-La chaleur de réaction produira une difference de température entre la surface limite 


et lintérieur de particules poreuses, au cours de réactions au sein de cette particule, Une equation 


est obtenue établissant la relation 


entre cette difference de température et la difference de concen- 


tration des composés réagissant Elle est indépendante de la cinétique de la réaction et de la 


forme géomeétrique des particules. Les exemples dapplication de cette equation sont donnés 


pour la désydrogénation du cy /ohexane et la régenération de catalyseurs 


Zusammenfassung Reaktionswiirme erzeugt eine Temperaturdifferenz zwischen dem 


Inneren und der Aussenseite cines porésen Teilchens wihrend ciner inneren Reaktion. Es wird 


eine Gleichung abygeleitet, um diese Temperaturdifferenz auf die Konzentrationsdifferenz des 


Reaktanten zu bezichen und es wird gezeigt 


dass sie unabhingiy von der Reaktionskinetik 


und der Teilchengeometrie ist Fir die Anwendung dieser Gleichung werden als Beispicle dit 


Dehydrierung von Cyclohexan und die Regeneration eines Katalysators behandelt 


INTRODUCTION 


Many of the reactions which occur within porous 
catalytic solids have an appreciable heat of 
reaction. This heat of reaction will give rise to 
a temperature difference between the inside and 
boundary of the catalyst parti le while the 
reaction Is proceeding. An estimate of the siz 
of this difference is needed in Many catalytic 
investigations. A complete solution of the problem 
would give the temperature as a function of the 
spatial co-ordinates of the particle. To obtain 
this, a separate solution is required for each 
particle shape, boundary condition and kinetics. 
It is not possible to obtain such solutions in a 
closed form exct pt for certain spec ial cases and, 
in general, numerical methods must be used. 

It will be shown that a solution can be obtained 
which is valid for all kinetics and all particle 
geometries. This solution gives the temperature 
as a function of the concentration of reactant. 
Information has been sacrificed to obtain this 
generality, namely, the temperature as a function 
of co-ordinates. This solution is : 


| AHD, 1) 
c c 
s 


where 7 is the centigrade temperature at all 
points within the particle which have a concentra- 


tion of reactant ¢ (mole /em”*) 


7, is the temperature (°C) at the surface of the 
particle 

c, is the concentration of the reactant at the 
surface of the particle (mole cm”) 

AH is the heat of reaction (cal mole) (negative 
values mean exothermic reactions) 

dD, is the effective diffusivity of the particle 
sec) 

A is the thermal conductivity of the particle 
(cal, sec em “C) 

DamkOu.Ler [1] and [2] have investi- 

gated the thermal effects in spherical particles 

and have obtained the same equation for this 

special case. The derivation given below will 

show that the equation holds for particles of 

any shape. 


The maximum possible temperature difference 
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between the inside and boundary of the particle 
can be obtained at once from this equation. The 
maximum temperature for a given operating 
condition will occur when ¢ has its smallest 
possible value. This will be the concentration 
of reactant at thermodynamic equilibrium 
between reactants and products, Equation (1) 
is a useful tool in cataly tic research and operations 
since the maximum temperature can be used as a 
criterion to determine whether a more detailed 
analysis is desirable for a particular investigation. 
Equation (1) is derived as follows : In conven- 
tional porous catalyst particle s, diffusion through 
the pore structure of the catalyst is the only 
method operating in supplying reactants to and 
removing products from the region of an active 
site within the particle. Thus, the diffusion 
equation : 
D,V* 0 (2) 
dt 


must be satisfied at all points within the particle. 
In this equation dn dt is the rate of disappearance 
of reactant at a given point (moles ‘sec em*). The 
heat produced (or consumed) by the reaction must 
be removed (or supplied) to the active site by 
ordinary thermal conduction processes. Therefore 
the thermal equation : 


AH =0 (3) 


must be satisfied at all points within the particle. 
These two equations, with the proper boundary 
conditions, completely determine the tempera- 
tures within the cataly st particle. To obtain a 
complete solution for these two simultaneous 
differential equations, dn/dt must be known as 
a function of concentration of reactant and as 
a function of temperature. In addition, the 
geometry of the particle must be known. How- 
ever, the path to be followed here in obtaining 
our solution does not require this information. 

Multiplying equation (2) by AH / and 
dividing (3) by K, 


4 
A dt 


The temperature produced by heat of reaction in the interior of porous particles 


285 


dn AH 
dt K 


=0 (5) 


Let « = AHD,/K 


Since « is a constant it can be taken under the 
differentiation. This gives for equation (4) : 


Vrae = 0 (6) 


Equations (2) and (3), which have sources 
and sinks of different strengths, have been 
transformed into two equations (6) and (5) which 
have the same sources and sinks. This means 
that the difference (f) between «xc and 7 must 
satisfy Laplace’s equation (V* f = 0). The solu- 
tion for equations (5) and (6) is then: 

aC T +f (7) 

The usual boundary condition is that the 
temperature and concentration of reactant at 
the surface of the particle are constants. Let 
their values be 7, and c, respectively. Then, 
on the boundary, equation (7) becomes :— 


T,+/f, (8) 


where f, is a constant. Any function which 
satisties Laplace’s equation in a finite closed 
region and is constant on its boundaries is 
constant in the region. Therefore, f = f, for all 
points within the particle. Equation (7) now 
becomes 


T (9) 


Eliminating f, between equations (8) and (9) 


and using the definition of «, 


AH D, 


T —T,) = 
( K 


(10) 
As can be seen the derivation is independent of 


particle geometry and reaction kineties. 


APPLICATIONS TO CATALYST SYSTEM 
Two examples will be given to illustrate its 
application to practical catalyst systems. The 
first system to be considered will be the regenera- 
tion of silica-alumina catalyst used in Thermofor 
catalytic cracking of gas oils. This catalyst 
is in the form of silica-alumina beads containing 
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10 wt. per cent alumina and with an average 
particle radius of 0-12 in. The beads are gradually 
fouled with carbonaceous materials in the reactor 
and in 20-80 min require regeneration by burning 
in air, During burning, the temperature in the 
kiln may reach 760°C. It becomes important to 
know how much the temperature within the 
beads can exceed the external t mperature, 
because the beads can be damaged by tempera- 
tures much in excess of 760°C. 

Coke burning within catalyst particles differs 
from the stationary state problem discussed 
above in that the reaction zone moves with time 
into the catalyst. However, it can be shown 
that for most cases of interest equation (1) still 
applies. The conditions for this are (1) that the 
amount of oxygen required to fill the pores of 
the catalyst is small compared to that required 
to burn the coke and (2) that the time it takes 
for maximum possible transient te mperature 
differences to decay to a negligibk quantity is 
small compared to the time it takes to burn the 
carbon on the catalyst. These conditions are 
satisfied for the case to be discussed. The value 
of dD, for oxygen at 760°C, determined by the 
gas flow method discussed by Weisz and Prarer 
(3), is 4-6 10° cm*/sec. The value of K, 
taken from the studies of thermal conductivity 
of various catalysts by Seur [4], is 8-6 10* 
cal/em sec °C. The heat of burning (AH) is 
that observed in practical kiln operation and 
has a value of 83 10° cal /mole. Equation 
(1) then yields a maximum possible t mperature 
difference of 1-1°C. This temperature rise is 


{1} ier G. Z. Phys. Chem, 1943 A193 16. 
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negligible with respect to further consideration 
of thermal damage to the catalyst. 

The next reaction to be considered is the 
dehydrogenation of cyclohexane over a commer- 
cial platinum-alumina reforming catalyst. This 
is an endothermic reaction with a AH _ of 

52-54 10° cal mole at 400°C. The condi- 
tions are 400°C reaction temperature, 4:1 
hydrogen to hydrocarbon ratio and 25 atmos- 
pheres total pressure. The values of the constants 


for this catalyst are : 
D, = 15-6 10% em sec* (cyclohexane at 
C) 
K 5-3 10° cal ‘sec em °C, 


The value of the equilibrium concentration of 
cyclohexane is calculated from free energy data 
given by Rossini et al. [5]. The calculated 
maximum temperature drop possible is found to 
be 53°C which will be quite important in kinetic 
studies. Thus further study to determine the 
actual temperature distribution in the particles 
may be desirable. 

In the above discussion, the problem of the 
measurement of the particle boundary tempera- 
ture has not been considered explicitly. Conven- 
tionally, the temperature measuring device is 
located externally to the particle. The question 
then arises as to how well the temperature 
measured by the device represents the actual 
temperature of the particle boundary. A con- 
sideration of the particular conditions existing 
in a given experiment is needed to answer this 


question. 
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Abstract 
range of heat flux from zer that corresponding t boiling s determined by 


The aver we covered area of a h ating surt e¢ immersed in boili ter, over a 


the resistance between the heating wire and clar le immersed he boilir vater. 
Th 


the heat flux was increased since two adiacen i es W lel unite d the two nuclei would then 


random boiling 


form bubbles at the same stave in their growth 


Theory indicates that the contact ar: between bul ind I rf m hich it is 


VOL. yrowing vari approxit ‘ inearly wit! ‘ i that iter: tl , rea covered 

Q by one source is Ipproni itely hal rea by 

ras active agents that increase the bubble film rth per ‘ r nur heat sources to 


exist on the heating surfac« 


Résumé 


La surface mo nee couverte d'un rfu muffante im uce eau bouillante 


dans un intervalle de fl thermiqu llant r lant l lilition a lm. 


eté determines par la résistance entre k wll ! in 4 trode 1 i dans leau 


bouillante. 


L¢ébullition désordonnée a faib! jues a ét t remplacée par une 
emission svnchrot ilors que le ff t ilies adjacentes 


des 


theorie indique que la surface de cor t , irfa pl ir | iclle elle 


l ii se 
dévek ppe varie a peu pres li im entay i iu ia i nnesc huverte 
per une source unl jue envi iu la surfiae la source 
meme. Les agents tifs a rfc ii nt resistan de lilm des bulles permettent 


lexistence d'un plus grand nombre de sources thermi + sur la surface chauffante. 


Zusammenfassung—Der im Mittel mit Blaser ‘ kte Anteil einer in siedendes Wasser 


getauchten Heizfliche wurde j ; ict \ von Nul zum Einsetzen der 


Filmverdampfung bestir t em der W tand ; hen dem Heizdra ind einer zweiten, 
Wasser 


in das s! cle nde 


Das 


zufillige Sieden i gerings VW iis fl vurde schrittweise ersetzt durch eine 


svnchronisierte Abgabe, ji ler luss er it wurde kl zwei benachbarte Blasen 


sich vereinigen wollten und die beiden Kerne dann B! i leichen Stadium inhres Wachstums 


formten 


Die Theorie ergibt, dass die Beriihrur ihe z\ chen einer Blase und einer ebenen Wand, 


auf der die Blase wiichst, etwa linear t der Z nst und dass bei Wasser die von einer 
Quelle im Mittel bedeckte Plache w ahr die Hilfte eréssten, von dersel t Quelle in Anspruch 
genommenen Fliche betrii Obertlichenak Zusii die die Festig t der Blasenwand 
steigern, lassen die Existenz ciner grésseren Zahl von Wiirme juellen en der Heizfliche zu. 


INTRODUCTION LEIDENFROST (1756) [1]. He is said [2] to have 


Tue first recorded obs« rvations of so called tilm observed t ie be laviour of small drops of water 


boiling, that is, boiling which tak place when th when placed upon a heated plate. The effect 
heating surface is covered c mpletely with an which Lerpenrrost noted and which now bears 


insulating layer of vapour, are attributed to his name, together with that of film boiling, and 
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the spheroidal effect, is the non-equilibrium 
condition set up when a boiling liquid and a hot 
surface are brought together with their tem- 
peratures differing by more than a_ certain 
determinate amount. A thin layer of vapour is 
formed between the two which persists until the 
temperature difference has been suitably reduced. 
Thus drops of water when placed on a hot stove 
are observed to maintain their spherical shape, 
not wetting the surface, and to dance on the 
cushion of vapour which supports them. 

The whole field of boiling heat transfer was 
first subject to critical examination by NuKtyaMa 
[3]. He introduced the technique of observing th 
process taking place on an immersed, heated. 
electrical resistance wire and by means of this 
was able to get heater temperatures in water 
under atmospheric pressure ranging from 100 to 
1100 C. He observed a maximum pomnt on the 
curve of heat flux versus temperature difference, 
which corre sponds to the transition from nucleate 


to film boiling, This occurred at a comparatively 


low temperature difference, in the region of 


10— 50°C. Following this the boiling mechanism 
became confused before, at higher temperature 
differences, a stable film of vapour was fornix d. 
More recently this technique has been subject 
to further refinement, and in 1948 Farrer and 
Scoran [4] produced a more detailed study on 
the same lines. They studied the boiling of water 
at pressures ranging from 0 to 100 p.s.i.g. from 
ele ctrically heated wires of chromel A, chromel C, 


nickel and tungsten. Fie. 1 is typi al of the data 


Specific heat flux 


+O 10 WO 1000 10000 


° 
Temperature difference, 
Pp 


Fic. 1. A typical boiling curve 
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described by Farsper and Scorau and can be 


divided into five distinct zones : 


(a) Natural convection region where no 


vapour ts formed at the heating surface. 


(b) Surface boiling region where bubbles are 
formed on the heater but collapse as they 


rise through the liquid, 


(c) Nueleate boiling region where bubbles are 
formed on the heater and rise through the 
liquid at the same time increasing in 


volume. 

(d) Unstable region of film boiling with the 
vapour film on the hot surface rapidly 
changing in shape. 

(e) Stable film boiling with radiation making 
an increasingly larger contribution to the 


total heat transferred, 


Extremely high rates of heat transfer can be 
observed in the nucleate boiling region. 

Explanations of the mechanism of the transition 
from nucleate to film boiling have, in general 
been restricted to observations made in the course 
of experiments which were not specifically de- 
signed to elucidate this point. SARUKHANIAN [5] 
has stated that as the heat flux from a heated 
surface to a boiling liquid increases, so the number 
of points on the surface at which bubbles are formed 
(bubble sources) increase until a point is reached 
where thev coalesce to produce tha phenomena 
of film boiling. The recent work of Westwarrer 
and SANTANGELO [6] which constituted a photo- 
graphic study of the whole range of boiling 
phenomena points to the fact that the number of 
nuclei or bubble sources increase until the 
bubbles are in contact with each other. This 
corresponds to the maximum in the heat flux 
eurve and beyond it no nuclei are observed. 
JaAKon [7] is more cautious and in his book Heat 
Transfer part 1, states — ** whether the spheroidal 
effect is always initiated by coalescence or is 
spontaneously initiated will be diflicult to decide.” 

Ronsenow [8] has made a study of forced 
convection surface boiling and from observations 
of the ratio of liquid to vapour produced in his 
apparatus concludes that a film is formed when 


the composition of the steam water mixture 
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adjacent to the heating surface reaches a quality 
of 0-8. No investigations were made of the 
actual processes taking place at the surface of the 
heater. 

Vos and Van Srraten [9] have examined the 
various effects of the composition of boiling 
binary liquid mixtures on the maximum point 
in the heat flux curve. In particular, they were 
able to observe that for a given liquuid combina- 
tion, the largest maximum heat {lux was measured 
for the mixture in which the bubbles of vapour on 
becoming detached from the surface reached the 
smallest size. 

This paper is conce rned with the growth of 
Vapour bubbles in contact with a heating surface 
and the effects which result from this and the 
number of active bubble sources on the heater, 
especially at high heat fluxes. These factors have 
considerable bearing on the nature of the transi- 
tion from nucleate to film boiling and the study 
of them has led to the deduction of a mechanism 


for this pre wW'CSS, 


THEORETICAL TREATMENT 


Both theoretical and experimental studies have 
been made of the growth of single tsolated 
spherical bubbles in a superheated liquid. The 
experimental results due to DerGaraBeDian [10] 
agree very closely with the deduced expressions 
of both Forster and Zuper [11] and PLESSE’ 
and Zwick [12]. In Appe ndix La similar approach 
to that of Puesser and Zwick [12] has been 


vrowth of a single bubble in contact 


made to the 
with a heating surface, and an approximate 
expression deduced for the area of heating surface 
covered at the base of the bubble. 

Appendix II contains a development of this 
treatment to give the average area covered by a 
single source of bubbles. During nucleate boiling 
the heating surface is covered with.a number of 
such bubble sources from each of which bubbles 
originate. Each bubble commences growing 
from a vanishingly small radius and grows in 
contact with the heater covering an ever in- 
creasing area of it until it becomes large enough 
to leave and rise through the Supe rheated liquuid. 
As it becomes disengaged, a new bubble is formed. 


The average area covered at a source of bubbles 


” 


mechanism of the transition from nucleate 


to film boiling 


is shown in Appendix II to be one half of the 
maximum area covered by a single bubble. 

It is a well substantiated fact that as the heat 
flux increases, so the number of bubble sources 
on the heating surface increase [5, 6, 13]. Visual 
observations suggest that this increase is related 
ina linear fashion to the heat flux [13,1 4]. From 
this it is directly apparent that adjacent bubble 
sources will be closer together the higher the heat 
flux. Now in a pure liquid, two bubbles will 
coalesce easily when placed in contact with each 
other [15 18] provided they are held in contact 
for a short period [19]. At high heat fluxes, with a 
large number of active sources on the heater 
surface, a situation where two adjacent bubbles 
attached to the heater are close enough to 
coalesce is easily visualized. On combining to 
form a single bubble the two sources will be 
denuded at one and the same time, and thus 
future pairs of bubbles will be formed in phase 
with each other. 

When the heating surface has on it a sutlicient 
number of sources such that, if all the bubbles 
reach their maximum size at the same instant, 
the whole surface area of the heater will be 
covered with steam, then film boiling can com- 
mence, The maximum area covered is equal to 
the sum of the maximum areas of the individual 
bubbles on the heater and is equal to the total 
area of the heater. Thus the processes taking 
place are visualised as follows. 

\s the heat flux increases there is a gradual 
increase in the number of active bubble sources 
and this means the individual growing bubbles on 
the heater get close together. In favoured spots 
adjacent pairs will coalesce and begin to operate 
in phase with each other. At the point where 
film boiling starts there are suflicient bubbles on 
the heater which all reach their maximum size 
simultaneously. Hence in this region the total 
heater surface is covered with vapour columns in 
contact with each other whose walls will thin until 
coalescence takes place and the characteristic 
Vapour film is formed. Thus the onset of film 
boiling is governed by the condition determining 
the maximum area each bubble can cover and the 


number of active bubble sources on the hot 


surface, 
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As has been shown in Appendix IT the average used in place of the variable resistance. The 


area covered by a single source of bubbles is steam evolved was r turned to the system by 
approximately one half the maximur rea any means of a glass condenser. At one end of the 
single bubble can cove r. Hence as the max mum boiler a 7 in. dia. 7 mm thick Py rex window was 
heat flux is approached the average covered , fitted to fa tate observation of the processes 


of a heating surface should approach on lf the taking place, and also to give visual indication of 

total available area. he onset of film boiline in a pote ntially unstable 
ystem. 

EXPERIMENTAL VERIFICATION The method evolved for measuring the average 


To demonstrate this point a series CX irface area of e heater covered by steam 
ments were designed and carried out |; 1i4 Phe, bubbles involves me urement of the electrical 
were aimed at measuring the average co lareca resistance betwe the heating surface and a 
of a heating surface immersed in boil \ cr secondary = elec immersed in the b viling 


over a range of heat flux up to that corresponding ‘ j resistan has been found to vary 
to the onset of film boili directly with the covered a of heater in the 
The apparatus consisted essential of a particular apparat used [14]. In order to make 
eviindrical boiler, with a Pvrex elas | p tl resistance measurements as sensitive as 
in Which distilled water was boiled on a possib ll the connections to the heater inside 
the boner were carefully insulated by giving two 

cou waterprool paint, leaving the heating 

<= \\ 1 t electrical contact with the 

liquid, Further, small additions of hydrochloric 

if elev f le 1 were made to reduce the liquid resistance. 

Material: brass Resistance measurements were made with a 


pica sult of these experiments is given 


in Fig. 3. This shows a linear variation of the 


8'disc of 
pyrex 
gloss 
7mm thick 


8 % bolts on } | 


Fic. 2. Side view of boiler in the un! 


wire 0-061n. dia. and 2 in. lone re is water 


heated by passing an alternating cur ip to 4c / 
300 A at 50c¢/s. The heat generated in wire 

G 3 | 
was measured by taking read £ if r 


flowing and also the pot crop bet we wo 
| 
4 
J 


potential taps brazed to the wire at about 2 


spacing. K 


| | 
The heat flux was adjusted by means of a rheo / 


stat in the primary circuit of the current trans O 100 a 300 400 
former in the higher heat flux range whil oo olve Heot flux, W 
close control in the lower range, a Vari wi Fic. 3. Effect of Teepol on nucleate boiling 
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average area of heater covered vs. heat flux for 
the boiling of distilled water. It also indicates 
that, at the onset of film boiling, an average of 
about half the heating surface is covered as 
predicted theoretically. A mean figure obtained 
as a result of several experiments at both atmos- 
pheric and reduced pressure gives the maximum 
covered area as 53-5 per cent for distilled water 
[14]. These experiments also show that the linear 
relationship between the number of bubble 
sources and heat flux which is observed at low 
heat fluxes [13, 14] persists throughout the range 
of nucleate boiling examined. 

In order that these ideas could be further sub- 
stantiated, additional confirmatory experiments 
have been carried out. In these the covered area 
of the heater was measured over a range of heat 
flux under nucleate boiling using solutions con- 
taining Teepol as the boiling liquid. This surface 
active agent has the property of reducing the 
surface tension of the liquid, and also of increasing 
the bubble film strength. The results are shown in 
Fig. 3. As can be seen covered areas much in 
excess of 50 per cent were observed. This is due 
to the bubbles not coalescing easily as they are 
compressed together at high heat fluxes. This is 
thought to be due to the increased film strength. 
Hence more such sources can exist on the heater 
before the isolating films break and film boiling 
starts and, as observed, the maximum average 
covered area is greater than that measured for 
water alone. Further, at high fluxes, the 
cramming of bubbles causes a departure from 


linearity. 


CONCLUSIONS 

By utilizing the existing knowledge of the 
processes taking place during nucleate boiling, 
and combining this with a programme of experi- 
mental work designed to a particular end, an 
explanation of the processes taking place during 
high heat flux nucleate boiling in a pool of liquid 
has been arrived at. Further, a mechanism is 
proposed and substantiated for the way in which 
nucleate boiling is induced to start. 

These mechanisms suggest that a clearer under- 
standing of nucleate boiling can only come from a 
study of the causes of nucleation, and further 
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that the heat transferred in such processes is 
intimately connected with the formation and 
gowth of individual bubbles. 

This behaviour is not confined to boiling but 
can occur in the production of froth for floatation 
in small bubble air lifts. 
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APPENDIX I 


Theoretical Treatment of the growth of a Single Vapour 
Bubble on a Horizontal Surface 

The theoretical approaches to the growth of vapour 
bubbles made by both Forsrer and Zuper [11] and 
PLEsset and Zwick [12] have only, strictly speaking, been 
concerned with the growth of single isolated sph rical 
bubbles. It is therefore necessary to amend the treatment 
somewhat when the case of bubbles growing in contact 
with a plane surface is to be dealt with. 

The highly idealised system which is visualised for the 
purpose of this treatment, consists of a spherical bubble 
with a cap removed such that it is in contact with the 
heater surface at an angk to its surface of separation 
(see Fig. 4). That is, the plane of the cut is replaced by the 


U 


\ 3 


hic. 4. A single bubble on a plane heating surface. 


plane of the horizontal heating surface. The bubble is 


assumed to retain this shape until it becomes detached 
from the heater. The temperature distribution in the 
superheated liquid is such that close to the heater there is 
an extremely thin layer of liquid at a high temperature 
(Jakos [7]}) Fig. 5. Above this layer, over a height equal to 
the maximum bubble diameter, the liquid exhibits a 
smaller but constant superheat due to turbulent mixing. 
Thus the environment resembles closely that envisaged by 
Piesser and Zwick [12]. This being so, it can be treated 
in a similar fashion with modifications for the change in 
shape. 

The work of Ronsenow and CLark [8] has shown that 
the flow of heat from a hot surface to a boiling liquid takes 
place first into the liquid in contact with it. Later the 
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Fic. 5. Temperature distribution in nucleate boiling from 


a horizontal flat plate [7 


heat flows into the bubbles through the liquid vapour 


interface. Thus, in this re spect, the presence of the heat 


f the 


can be ignored, no « vaporation « 


heat 


ny surface Vapour 
taking 


bubbles. 


place due to conduction of through the 

The treatment of this problem which follows is based on 
the simplified physical deductions due to Pursser and 
Zwick [12 


original work 


This form of argument, in the case of the 


gives an approximate answer, differing by 


3\? 
only the constant factor | O77 from that 


calculated by a more rigorous mathematical treatment: 
it is therefore considered adequate 

Consider a spherical bubble growing in contact with a 
plane surface as in Fig. 4. The area of curved surface. 
taking the 


sphere, is given by: 


bubble as a spherical segment or capless 


Ilere A 
The solid, liquid, vapour contact angk mal 


By integrating we get : 


The area of curved liquid ‘vapour interface, 
The 


radius. 


cos 8). (1) 


! is the surface area from which evaporation can take 


place and replaces the total bubble area of t- used by 


Piesser. The volume V of such a bubble is viven by 


(2) 


This compares with the normal volume of 


3 


and 


a sphere of 


HASLAM 
Now, assuming [12] that the temperature drop takes 
place principally in a liquid layer of approximate thickness 
the diffusion length (at)} where a thermal 
The rate of heat flow 


given by 
diffusivity in the liquid 


dQ kATA 


dt (at)! 2 


Here hk 
thermal diffusivity, and AT 
At the same time the heat flow required for evaporation 


The liquid thermal conductivity, a The liquid 


The liquid superheat 


is: 


d 
(4) 
dt dt 

Here L The 


density 


The latent heat of the liquid, p 
and J 
equating equations (3) and (4): 


vapour 


The bubble volume 


(at)? 


Now 


3 cos 8) p, from (2) 


di 


rs 2 
at 


8 cos 8) 


since 8 does not vary with time and the vapour density 
can be assumed constant, as has been done by Piesser. 
Substituting in equation (5) we get : 
KAT 


(at)? 


3) 


> 
COs p) 


kAT 2 (1 


cos" 3 cos 8) 


(at)4 Lp, (2 


That is, for a particular superheat and contact angle : 


kAT 2 (1 


ad Lp, (2 


coos B) 
Ilere ¢ constant 


cos 


5 


This is a similar expression to that obtained by Puesswer, 


but in his case the factor C is given by : 


kAT 1 
ad Lp, 


Integrating equation (6) we get 
r=2cA (7) 


Now the area of the base of a bubble in contact with a 


plane surface is 7 r* sin® 8, and substituting from equation 


101-2 
Of} — + —{.| 
100-4! —_+—Water 4+Steam— 
j 
100-0 
VOL 
L (Vp...) (5) 
195 
d 
i 
Lp, wt? (2 — cos® p + 
P and dt 
dr 1 
( (6) 
dt 
— <= 
r 
l 
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(7) gives 4nc? 8, i.e., the contact area between the 
bubble and a plane surface on which it is growing varies 
approximately in a linear fashion with time. 

A similar expression for the surface area of a spherical 
The curve in Fig. 6 shows the square of the radius of a 
spherical bubble plotted against the time from the work 


of DERGARABEDIAN [10]. As can be seen, this is very close 


bubble and its dependence on the time ist 
a 


to being a straight line, and it gives added confirmation 


to the theory. 


15} 


2 
x 


Bubble rodius squared, 


| 
10 15 20 


Time, sec x10” 


Curve of (bubble radius)- vs. time. 


II 


Formation of Bubbles from Single Sources During Boiling 


In Appendix I it has been shown that, as particular 
bubbles grow on the heating surface, the surface area 
which they cover increases linearly with respect to time. 
Consider now the processes taking place and the mean area 
covered at a particular source of bubbles. 

The bubbles will each grow in contact with the heating 
surface until they reach some equilibrium size. Then they 
become detached and rise through the hot liquid. The 
factors governing the maximum size such bubbles can 
reach have been very fully examined by Frrrz [21) 
following on from the papers of Basurorrn and AbAmMs 
[22] and Wark [23]. The results of this last author can 
be expressed both in the form of an equation and graphic- 
ally. The equation for the maximum volume V’,, of such 


a vapour bubble is :- 


The mechanism of the transition from nucleate to film boiling 


2938 


a2 
00-1198 
£ (py Py) 
Here JV,, Maximum bubble volume, 8 Contact 
angle, o Liquid-vapour surface tension, p; Liquid 


density, and p, Vapour density. 

After the bubble has left its parent surface, a fresh one 
is formed at the same spot. WestTwarer and SANTANGELO 
6] have observed that in boiling methanol the bubbles 
form immediately after each other and the present work 
agrees with this conclusion. Jakop [20] however has 
observed from his photographic work that the formation of 
new bubbles is delayed by about half a period or 1/40 see 
after the previous bubbles departed. This can take place 
at very low heat fluxes. Then the flow of cool water into 
the place formerly occupied by a bubble will cause a 
momentary cooling: however the effect disappears as soon 
as full nucleate boiling sets in 

The average behaviour of a single source of bubbles can 
be described as indicated in the curve in Fig. 7; this is a 
plot of the area covered by the bubbles against time at 
any heat flux. This diagrammatic curve approaches the 


sawtooth form in Fig. 8. 


Time 


Fic. 7. Curve showing the area covered by regular bubbles 


at a single source. 


Time 


Fic. 8. Approximation to Fig. 7 used for calculation 


of average covered area. 


Considering Fig. 8, the average area covered over a 
whole number of cycles is the same as the average area 
covered during one cycle. Hence by inspection of Fig. 8, 
A 44m, where A Mean area with respect to time, 
and Am Maximum area. 

That is, the average area covered by one source is approxi- 


mately one half the maximum area covered by the same 


source, 


| 
13} 
12 
VOL. | 
~ 
] 958 re) j 
7 
| 9 
5) 
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| 
2} 
| 
9 
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NOTATION 


liquid thermal diffusivity r = bubble radius 
area of bubble liquid vapour interface time 

mean area covered by a bubbk liquid superheat 
maximum area covered by a bubble 
bubble volume 
liquid, vapour, solid contact angle , 
= maximum bubble volume 
constant 

liquid vapour surface tension 


liquid thermal conductivity 
latent heat of liquid vaporisation _ = liquid density 
heat flux P, = vapour density 


REFERENCES 
Leipenrrost J. De aquae communis nonnullis qualitatibus tractatus. Duisburg 1756. 
Daew T. B. and A.C. Trans. Amer. Inst. Chem. Engrs. 1937 33 449. 
Nukryama 8. J. Soc, Mech. Engrs. Japan 1934 37 367. 
Farper EF, A, and Scoran R. L. Trans. Amer. Soc. Mech. Engrs. 1948 70 369. 
SARUKHANIAN G, Chem. Ing. Tech. 1953 25 477. 
Wesrwarer J. W. and Sanrancero J. G, Industr. Engng. Chem. (Industr.) 1955 47 1605 
Jaxon M. Heat Transfer Part 1. John Wiley, New York. 1950. 
Rousenow W. M. and CiarK J. A. Trans. Amer. Soc. Mech. Engrs. 1951 73 609. 
Vos A. S. and Srrauen 8, J.D. Chem. Engng. Sci. 1956 5 50. 
DERGARAREDIAN P. J. Appl. Mech. 1953 75 537. 
Forster H, K. and Zuper N. J. Appl. Phys. 1954 25 474, 
Piesser M.S. and Zwick S. A.J. Appl. Phys. 1954 25 493. 
Jaxon M. Res, Publ. Illinois Inst. Tech. 1942 2 159. 
Hastam F, Ph.D. Thesis University of London 1956 
Maier C. G. U.S. Bureau of Mines Bull. No. 260 1927. 
Verscnoor H. Trans. Inst. Chem. Engrs. 1950 28 52. 
C, W, and J. M. Industr. Engng. Chem. (Industr.) 1931 23 1283. 
Fou.k C, W. and Berkiey J. Industr. Engng. Chem. (Industr.) 1948 35 1013. 
Darra R. L., Narrer D. H. and Newrrr D. M. Trans. Inst. Chem. Engrs. 1950 28 14. 
Jaxon M. Z. Ver. Dtsch. Ing. 1982 76 1161. 
Farrz W. Phys. Z. 1935 36 379. 
Basnrorru F. R. and Apams J. Capillary Action Cambridge 188%. 


Wark I. W. J. Phys. Chem. 1933 37 623. 


: 
A 
Am 
k 
L 
Q = 
(1) 
[2] 
(3) 
(4) 
[5] VOL 
[6] 
195 
[8] 
[10) 
[12] 
[14] 
[15] 
[16] 
[17] 
18] 
19) 
(20) 
[21] 
(23) 
294 


ineeru Science, 1958, Vol. pp. 295 to 408 Per mn Lid 


Chemical 


The mechanism of mass transfer of solutes across liquid-liquid 
interfaces—II] 


The transfer of uranyl! nitrate between solvent and aqueous phases 


J. B. 


Chemical Engineering Division Atomic Energy Research Establishment, Harwell 


LEwIs 


1957) 


(Received 9 August 


Abstract—~ Overall’ transfer coeflicients were determined for the transfer of uranyl nitrate 


hetween water and the three solvents — dibutoxy diethyl ether (D.B.C.), methyl butvl ketone 


(Hexone) and a 20 per cent solution of tri n-butyl sphat I. B.P.) in * odourless kerosene 


(OW benzene or carbon tetrachlorid In sore er nts with D.B.C. the aqueous phase was 
nitric acid in equilibrium with D.B. it 1-SN it itric acid These experimental coefficients 
were compared with values pr ted as f ! nterfacial resistance and using individual 


transfer coeflicients calculated from the corr 


Some svstems gave larger transfer rates than predicted and in all of these interfacial tur 


bulence was observed Tm the ther stems the initial lues of the calculated and observed 


coeflicients were usual i | reel but the observed rate tended t liminish with tin 
apparently due to the buil p of i facial rrier The for tion of this barrier appears 
to depend solely on the dur 1 of eriment, bei independent of uranyl nitrate con- 


centratior 


The transfer of nitric acid alone between water and D.B.C. was investigated. No interfacial 


resistance was found when the tr sfer was from water to D.B. the transfer in the other 


direction was accelerated by interfacial turbulence 


Résumeé— Des efficients « tr fert lobaux mt été déte inés pour le transfert du 
et ! solwar nts dibutoxvdiethvether (D.B.C.) 


nitrate 
methvlisobut vicétone (1 et ut butviphosphate (T.B.P.) dans du 


certaines 


en équilibre 


ont 


let 


ent ete 


rcs les vale ‘ COS ‘ re terre if fact et utilisant 


coeflicien 


celles prévues 


Quelque 


et dans tou ‘ un n ervet les autres svstémes les 


valeu la 


thesse ce tendait a di ‘ i était i apparemment a la formation 


dune barrits fou ition i ipparait n rele que de la durée 
de Ve perience elle est i ken lante de la wentratio du nitrate d ranyle et de toutes les 
autres variables experi ntales 


seul entre DBC, a été étudi Aucune résistance 
le 


Le transfert de lacide nitri 


transfert dans 


interfaciale n'a été trouveée, quand le transfert rit leiu de Teau vers D.B.C, 


witre direction ctait acceléeré par turbulence inter ciale. 


Zusammenfassung l wurden Stoffdurchgangskoefli nten bestimmt fiir den Ubergang 
von Uranylnitrat zwischen W er und den drei Lésemitteln Dibut vdiiithvlither (D.B.C.), 
Methylisobutylketon (Iexor nd einer 20° en Lésung Tri-n-butylphosphat [T.B.P.] in 

eruchsfreiem Keroset OK Ber | oder Tetrachlorkohlenstoff. In einigen Experimenten 
mit DBA bestand die wiissrige I © al 8 n-Salpetersiiure im Gleichgewicht mit D.B.C. in 
1.5 n-Salpetersiiure Diese experimentellen Koeflizienten wurden mit den ve rausberechneten 


verglichen, wobei kein Zwischenflichenwiderstand angenommen wurde und die einzelnen 


Obergangskoeflizienten aus der in Teil I [4] abgeleiteten Beziehung berechnet wurden. 
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* kérosét inod du beng du trachlorure de carbone Pe 
périences réalisées ec le DBA pl at était de Vacide nitrique 3 N Ei 
avec le D.B.C, dans lacide nitrion btenus exm 
les 
transfert indi ‘ caleul i partie de | n derive ns I (4). 

systémes dor td transfert plus grande 


inige Svstenn 


Zwischentlichenturbulenz 


standes verursacht wurde 


Verinderlichen 


Auch die Ubertra 
cin Zwischentlicher 


In Part I (4) the factors influencing the 
binary solvent—water svstems were 
ina specially designed transfer cell, 

correlation between the 
Re vnolds rs ind Vviscositt 
transfer coellicrents to calculated for 


to 


between water and a number of oreat 


greater than pre dicted. 
phenome non 


he te rogencous 
reaction at the 
st« p 


reasonably be « xpected., 


solvent interface must necessarily be 


If the unionized molecule is the sper 
then during trans 
» aqueous to solvent phase, 


equilibrium will beome 


of transfer will 


Ubertragungswerte als vorausgesagt: in all diesen 


obachtet 


getnessenen mit den berechneten Werten in guter Chere instimmung 


verkleinerten sich mit der Zeit is Vermutlich durch den 


von der 


tersiiure allein zwischen Wasser und DBA 
widerstand gefunden bei Ubertra 


die Ubertragung in der anderen Richtung wurde durch Zwiselh nflichenturbulenz erhdlht 


ndividual 


ration of 


estigated 
itisfactory 
ents and 
of the two 


maliy idual 


inv two 
thon was 
rall Coe 


solute 
solve nits. 
Ny ri 
ith 
nts were 


small 


The fast transfer rates wer found ¢ xperiment illy 


interfacial 
cs were 
che mical 


The transfer of uranyl! nitrate across thy water- 


ite d 


with a chemical interfacial ste p as the molecule is 


and oft n 


compl xed with the solvent, in the solvent phase. 


es a tually 


fer from, 


the aqueous phase 


interface 


and more unionized molecules will be continuously 
formed by association of the UO,* and ).” lons, 


influenced by 


the rate of association which is propo tional to 
[UO,] Né Either some type of three body 


collision must occur or, more likely, successive 


B. Lewis 


Uranylnitrat-Konzent ration 
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Fallen 
inderen Systemen waren dic anfanyglich 
Aber die berechneten Wert« 


Aufbau cines Zwischentlichenwider 


Die Bildung dieses Widerstandes scheint nur voa der Versuchsdauer 


und allen anderen 


wurde untersucht 


ung von Wasser zu 


reactions via an intermediate such as | 0, NO, )*, 
This 


possibly 


process will bo followed by a second and 


more important step, the combination 


ot the unionized molec ule s with solve nt mole cules 


to form a compl x If ions cross the interface 
then the combination will occur in the solvent 
phase . Transfer in the reverse direction will be 


accompanied by a dissociation step. 
If this 
it to be 


chemical step is sufficiently slow for 
| 


significant under conditions of high 


turbulence, then it may set an upper limit to 


thre of solvent extraction contactors. 


Although h has publishe il on the 


the rine 
dynamics of the distribution of ut inyl nitrate in 


solvent—-water sv stems, verv little work has been 


carried out, other than of a purely technical 
nature, on the rate of reaching equilibrium. 
Murpocn and Prarr have studied the transfer 


and the two 


solvents dibutoxy di thyl ether (dibutyl earbitol 


of uranyl nitrate between water 


or D.B.C.) and thyl iso-butyl ketone (hexone) 
ina wetted wall column [7]. They found evidence 
reaction in all cases and 
the 


of a slow interfacial 


derived approximate values for velocity 


constants. The reaction appeared to be of a 


pose third orde r. More recently SMITH, 
and Prart [8] have investigated the 
transfer of uranyl nitrate between water and 


hee xone im a packed column and compared the 


experimental results with those predicted from 
individual coeflicient data obtained from binary 
They that 


between the predicted 


ethyl acetate -water systems. found 


satisfactory agreement 
and the experimental values could only be ob- 
tained if allowance was made for the interfacial 
resistance, using the values obtained by Murpocu 
The work 
published is that of Haun [2] who studied the 


and Prarr. only other relevant 


al md war unabhingig 
VOL 
Q 
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== phase svstem im thus cell This corr 
wed, in Part 
eflicients for the tr 
These coeflicients were compared w 
mental values and it was found tha 
instances the actual transfer coelher 
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tur 6 vi tT low 
| 
er 
si 
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transfer of uranyl nitrate between water and 
tri n-butyl phosphate (T.B.P.). Both phases 
were stagnant so that the transfer to and from 
the interface took place by molecular diffusion. 
No interfacial barrier could be detected. 

In the present investigation the overall rate 
coetlicient for the transfer of uranyl nitrate 
between water and the three solvents, D.B.C.., 
hexone, and T.B.P. were determined. The T.B.P. 
was diluted to 20 per cent by volume, with 
benzene, carbon tetrachloride, or “ odourless 
kerosene.” A few transfer rates we re also measured 
in the water—D.B.C. svstem with nitric acid 
present in both phases in order to investigate 
the effect of the latter on the transfer. The 
transfer of nitric acid alone between water 
and D.B.C. was also investigated. The T.B.P. 
experiments were carried out at 25°C and most 
of the D.B.C. and hexone experiments at 20 ( 
since these were the temperatures at which most 
equilibria and extraction data have been re porte ad. 
A few experiments with D.B.C. and hexone wer 
done at 50°C, in an attempt to measure the 


activation energies of the processes, 


EXPERIMENTAI 


The transfer cell was the same as that used in 
the previous experiments and described in Part I, 
Each compartment had a volume, V. of 300 ml 
and the area, A, of the interface was 30 em. 
The two phases were stirred independently, 
the rates being expressed as Reynolds numbers 
L? Nv, L being the length of the 
stirrer, N the revolutions per second and v the 


defined as Re 


kinematic viscosity of the liquid. Electrodes 
were inserted into the upper (solvent) compart- 
ment and at first the concentration in the solvent 
phase was determined conductometrically. 
Although this method was sensitive and enabled 
continuous observation to be carried out it was 
abandoned as erroneous values were occasionally 
obtained due to the presence of moisture on the 
polythene’ electrode-holder. About six small 
samples (0-1 to 1-0 ml) were therefore taken at 
intervals from the two phases, these volumes 
being too small to affect appreciably the condi- 


tions in the cell. The uranium concentration was 
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determined colorimetrically using the method 
described by Murpocu and Prarr. 


Systems used 

The dibutyl carbitol was a commercial product 
containing about 4 per cent of organic impurities, 
mainly monobutyl earbitol and 8 butoxyethyl 
butyl acetal together with various vinyl ethers 
and organic peroxides. The crude solvent was 
heated under reflux with N sulphuric acid for 
6hr in order to decompose the acetal and the 
ethers to aldehyde s. which were removed by a 
subsequent evaporation for 30 min. The dibutyl 
carbitol was then neutralised with sodium 
carbonate solution. washed several times with 
water, dried over calcium chloride and fractionally 
distilled under reduced pressure. Monobutyl 
carbitol and most other impurities were removed 
by the distillation process. The main fraction was 
then washed with successive batches of alkaline 
permanganate solution until the permanganat« 
was no longer reduced. After water washing 
and drying, the dibutyl carbitol was fractionally 
distilled a second time. Finally, the pure solvent 
was percolated through a column of activated 
alumina to remove any peroxides and surface 
active agents. The D.B.C. was stored in the dark 
in contact with freshly reduced copper foil. 
Before use the solvent was given a second run 
through an alumina column. 

The methyl iso butyl ketone was a commercial 
product which was purified by distillation in a 
column packed with 15 ft of 1 in. Berl saddles. 

Commercial T.B.P. was available which was 
approximately decinormal in acidity due to the 
presence of free dibutyl hydrogen phosphate, 
(D.B.P.), Other major impurities present were 
n-butanol, tetrabutyl pyrophosphate, butyral- 
dehyde and the ammonium salt of D.B.P. Two 
methods of purification were employed. 

The first method consisted of refluxing the 
T.B.P. with four times its own volume of 5 per 
cent caustic soda. This hydrolysed the pyro- 
phosphate and extracted the dibutyl phosphate 
formed, together with the original D.B.P. and its 
ammonium salt, into the aqueous phase. During 
the refluxing about one-third of the aqueous 
phase was allowed to distil away, thus steam 


| | 
- 
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distilling away the butanol and butyraldehyvde. 
Finally, the T.B.P. was well washed with water. 
Solvent treated by this method is referred to below 
as steam stripped T.B.P. A major disadvantage 
of this method is that the high boiling impurities 
are not removed. 

The alternative method of purification was to 
wash the T.B.P. successively with 4 per cent 
caustic soda until all the acidity had been removed 
then to wash with water, dry and ft tionally 
distil under high vacuum in the presence of 
anhydrous sodium carbonate. The T.BLP. 


purified by this distillation technique gave rise to 


‘ interfacial turbulence ” in contact with aqueous 
uranyl nitrate unlike the steam stripped T.B.P. 
The T.B.P. solutions used contained 20 0.5 per 


cent by volume of T.B.P. in odourless kerosene. 
benzene or carbon tetrachloride. 

lhe odourless kerosene diluent was an aromatic 
free high boiling fraction (200-2506 with a 
density of 0-79. For use with steam-stripped 
T.B.P. the kerosene was merely percolated 
through an alumina column. but that used to 
dilute distilled T.B.P. was itself distilled under 
high vacuum. The effect of using distilled odour 
less kerosene with untreated T.B.P. was also 
investigated. 

Benzene and carbon tetrachloride of Analytical 
Reagent quality were distilled without fractiona 
tion before use. The benzene used in Expt. 653 
was super fractionated under a high—tlux ratio, 
in an Sft still of lin. diameter, packed with 
Fenske helices. No significant difference in the 


transfer coetlhicoents was found between t] 


this benzene and the original A.R. benzene. 

The uranyl nitrate was a commercial product 
purified by ether extraction, while the nitric 
acid was of A.R. quality. Demineralized water 
was used for the aqueous solutions. All solutions 
were percolated through an activated alumina 
column before use in order to remove any surface- 


active impurities, 


Equilibrium and other physical data 
Equilibrium data at 20° have been published 
by Murpocu and Prarr for the two svstems 


uranyl nitrate —dibutyl carbitol— water and 


uranyl! nitrate — hexone-water at 20°C. All other 


Lrewts 
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equilibrium data were determined experimentally 
using the same technique as adopted by these 
workers. The viscosities and densities of all 
phases were also determined. The distribution 
data and physical propertic s will be published in 
detail elsewhere. 

Interfacial tension measurements were carried 
out on all systems using the pendant drop 
method of ANpreas and Hauser (1). This 
technique also enabled the presence of interfacial 


turbulence to be detected. 


Rate determining processes 

The transfer of uranyl nitrate from a solvent 
to an aqueous phase involves transport up to the 
interface in the solvent phase, transfer across 
the interface itself and transport away from the 
interface in the aqueous phase. Since there is no 
accumulation at the interface the rates of these 
three steps will be equal and three transfer co- 
eflicients can be defined such that each rate equals 
the product of a transfer coetlicient and the 
appropriate concentration ** driving force.” It ts 
known that if the two phases are agitated there 
will be fairly uniform turbulence in the bulk of 
each phase but the turbulence is considerably 
diminished in the vicinity of the interface. 

The conventional theory used to calculate the 
rate of mass transfer across a_ liquid—liquid 
interface assumes that two laminar layers exist, 
one mn each phase, adjacent to the interface in 
which turbulence is completely absent and mass 
transfer takes place sole ly by molecular diffusion. 
Since this is a much slower process than the eddy 
diffusion which is obtained where the flow is 
turbulent the major resistance in transfer to the 
interface in each phase lies in transport across 
these boundary layers. It is a convenient and 
reasonable approximation to assume that the 
concentration gradients across the lavers are linear 
and that the concentrations are uniform in the 
bulk of each phase. The small additional re- 
sistance due to eddy diffusion is allowed for by 
taking the effective thickness of the laminar 
layers to be somewhat greater than the true 
hydrodynamic values. 

In Parts I and II it was shown that, under the 
particular experimental conditions used, transfer 
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of mass to and from the interface took place by 
eddy and not by molecular diffusion, there being 
no laminar films since the turbulence. although 
diminished, extended to the interface. There 
was thus no qualitative difference between the 
transfer rates in the bulk of each phase and near 
the interface. However. an expression for the 
transfer rates exactly similar to that applicable 
in molecular diffusion controlled systems can be 
applied. When transfer is from solvent to water 
the two liquid phase transfer rates are given by : 


u 


With molecular diffusion controlled systems 
the two transfer coefficients are equal to the 
appropriate solute diffusivity divided by the 
thickness of the effective boundary layer. No 
simple relation holds when the transfer is entirely 
by eddy diffusion since the coefficients refer to 
transport across the whole liquid and not just to 
transport across a more or less well-defined laver. 

The transfer of uranyl nitrate across the actual 
solvent—water interface itself can be regarded as 
being a heterogeneous reaction between the 
uranyl nitrate solvent complex and water. As 
it is difficult to deal with such heterogenous 
reactions it is customary for practical purposes 
to consider such reactions as being equivalent 
to a homogeneous reaction eccurring wholly in 
one or other of the two phases. In most instances 
as the interfacial process consists of a reaction 
between a solute molecule and one or more 
molecules of the solvent and or water. which 
are both present in large excess, the reaction 
is of pseudo-first order. Since, however, uranyl 
nitrate dissociates into three ions in water, a 
third order reaction might possibly be found 
with this compound and such a reaction order 
was in fact reported by Murpocu and Prarr when 
evaluating the interfacial reaction in terms of an 
equivalent aqueous phase reaction. Since uranyl 
nitrate is unionized in the solvents used. being 
complexed with two molecules of the solvent, 
an equivalent solvent phase reaction would 
probably be of a pseudo-first order and thus 
more easily calculated. The two possible homo- 
geneous reactions are of course exactly analogous, 


the concentrations being related by means of the 
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distribution curve which is approximately cubic 
The equivalent homogenous solvent phrase 

reaction is expressed by an equation : 


Here 


c,;* the solvent phase concentration in « quilibrium 


is the reaction velocity constant and 


with the aqueous phase concentration a. oi 
can be seen that when k,, is very large ¢,, and 
¢,; are in equilibrium, 

Since in general the interfacial concentrations 
are not known, neither the individual film 
transfer coetlicients k, and k, nor the reaction 
velocity constant k, can be determined: it is 
necessary therefore to define an * overall ” transfer 
coelcient which can be experin ntally deter 
mined and used to measure the rate of transfer. 
This overall coeflicient is defined on the basis 
that the whole resistance to transfer is a diffusion 
resistance in one or other of the two phase s. Thus 
for transfer from solvent to water, if the resistanc« 
is assumed to be in the solvent phase : 


) K,,, (ce, (3) 


here A, 


cocthcient and e.* the solvent phase concentration 


is the solvent phase overall transfer 


in equilibrium with the uniform concentration. 
C.» In the aqueous phase. Analogous equations 
apply for transfer in the reverse direction and 
for the definition of an aqueous phase coctlicient 
Determination of the experimental overall mass 


transfer coefficient 


The overall solvent phase coeflicient can 
obtained as follows: at any time ¢ the rate of 


transfer is given by : 


henee from equation (3) : 


A dt (e, 


(4) 


OL . 
Q 
1958 
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A dt 
integrating over the time interval ¢ : 
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Since samples were taken from each phase during 
an experiment the variation of c, and c, with 
time was known. For each value of Cc. the cor- 
responding equilibrium solvent phase concentra- 
tion c.* was obtained from the distribution curve. 
Thus 1/¢, — ¢,* could be tabulated for various 
values of t enabling A,, to be obtained by 
graphical integration of equation (4) for each 
time interval. 

It can be shown from equations (1), (2) and (3) 
that the overall coeflicient is related to the 


individual coeflicients by the equation : 


l H 1 


where H is defined as H 


The relevant concentrations for transfer from 
solvent to water are represented graphically 
in Fig. 1 (a). Curve ABPJ is the equilibrium 
distribution curve between the two phases whilst 
point O represents the two bulk phase concentra- 
tions ¢, and ¢,, at any time ¢. Point J representing 
the two interfacial concentrations ¢,,; and ¢,, 
lies off the distribution curve because these 
concentrations are not in equilibrium, point P 
represents the solvent phase concentration, 
¢,,*, which is in equilibrium with the aqueous 
phase interfacial concentration. The “ driving 
force ” for the equivalent solvent phase reaction 
is IP MN (c.. c,,*). It can be seen 
that the solvent and aqueous phase diffusional 
driving forces are represented by OM = (c, — ¢,,) 
and NP = (¢,; — ¢,) respectively. From simple 
geometry NP = NB/H where H, defined as 
in equation (5), is the slope of chord BP and 
NB (e,* c,*). Hence the total driving 
force OM + MN + H & NP equals OB which 
is also defined as the “ overall’’ driving force 
(c, — c,*). 


Equation (5) can also be written in the form : 


(6) 
K can be calculated since both k, and k, can 
be computed from the correlation derived in 


Part I, and H can be obtained from the distribu- 


| 0, 
(a (a) 
me 
Cc; = 
A Cw Cwi 


| (b) 
p 
4 Cw, Cw 


Fic. 1. The relation between transfer resistance (a) trans- 
fer from solvent to water; (b+) transfer from water to 
solvent. 


tion curve as shown below. Thus k, can be 
determined when the experimental overall co- 
eflicient A,. is known. 


Determination of the effective distribution co- 
efficients H and H’. 

Referring to Fig. 1 (a), H has been shown to be 
equal to the slope of the chord BP. Points B 
and P lie on the intersections between the 
equilibrium curve and the two perpendiculars 
drawn from O and J to the abscissa. The position 
of O, which represents the bulk phase concentra- 
tion, is known but the exact position of J is not 
known. If there is no interfacial reaction the 
interfacial concentration would be given by L 
which can be calculated, since the slope of OL 
(= —k,/k,) can be computed, and L lies at the 
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Table 1. Transfer of uranyl nitrate between DBC and water 


Mass transfer coeff. 
10%) 


Reynolds Interfacial 
Time resistance 


number 
Calculated Obsd. Bulk phase Interfacial (mins) Ri 


(cm /sec Concentration (mg/ml) 


Temp. 


Expt., 


202 


Water 


Solvent 


2000 


2000 


1000 


2000 


wl 


(a) Transfer from solvent to water 


(b) Transfer from 
2-58 


384 2: 


4-64 


(d) Transfer 


250 314 1 


73 


water to solvent (interfacial turbulence observed) 


S87 


| 
| 


48 


1-50 
1-20 


2-08 
0-43 


O12 
0-65 
0-58 


1-18 
1-13 
1-56 
1-18 
0-90 
0-66 
1-56 
1-21 | 
1-06 
0-67 


| 2-67 


0-67 
0-33 


1-60 
1-54 
1-21 
0-92 
1-04 
0-62 


61-2 35-2 6o-7 
75°3 
80-0 
83-6 


2-00 


12-0 
140 


«101 
104 
108 
110 


iw 251 


278 


282 


276 
ore 


271 


260 14-0 54 


» 

oe 
235 
229 


~ 


267 


to te tye 


te 


59-2 


0-00 


250 o-0 238 
221 
206 
196 
187 
181 


370 
30-0 
28 


5-00 
8-00 
10-0 
13-0 


(c) Transfer from 1-5N, D.B.C. to 3N water 
2-62 | 1: 


55-9 
53-9 
53-0 


52: 


11-9 
29-5 
44-0 
53-6 
63-2 


69-0 


19-0 


from 3N water to 1-5N D.B.C, (interfacial turbulence obsd.) 
33 


136 


(sec 3) 


— 0-127 

— 0-108 
0-074 
0-34 
0-25 


0-86 


k, | | Ge | | Se | Go 
2-85 267 830 3-00 25 0-05 
21-1 50 0-09 
16-4 85 0-26 
12-8 127 0-65 
3-08 282 251 | S1 820 4:50 220 OBS 
273 20 2000 | 2000 3-33 | 215 27 0-02 
0-42 35 1-96 
JOL. 0-67 55 1-04 
Q | 1-00 a 7 | O-44 
0-55 96 | 1-36 
2000 1-00 1-29 1-26 3-17 93-2 41-0) 131 3-1 | 0-48 
1-17 38-0 20 0-06 
1-67 34-6 30 | 0-19 
1-00 32-4 0-21 
249 20 2300 = 200 58-0 | 9 7-70 
= 14 0-39 
a 18 0-64 
24 O-75 
15-5 33 1-79 
276 20 1000 | 4000 | 5-33 | 2-57 20-3 203 60-2 8 | 0-46 
28-6 | 0-35 
| 
| 35-0 28 0-25 
| 39-5 38 0-46 
42-5 48 | 0-73 
| | 15-4 63 1-12 
535 BO 4000 O00 (230 0 6-0 222 34 10 0-14 
10-0 20 0-33 
| | | 181 | 30 | 0-44 
| | | 10-1 40 | 1-00 
| 
553 20 | me (260 | 10 | —0-30 
| 5-25 | 20 | 0-82 
| 6-15 | 80 | 5-25 : 
775 | | 40) 
| 
511 20 2000 | 7 0 5 | 
| 13 | 
| 23 | 
33 | 
| | 43 
53 
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Table 2. Transfer of uranyl nitrate between hexrone and water 


Vass transfer co 
Interfacial 
Time resistance 


(mins) 


(em Concentration (mg/ml) 
Temp Reynolds 
Calculated Obsd Bulk phase Interfacial 
(sec om Wr’) 


Water Solvent 


(a) Transfer from hexone to water (interfacial turbulence observed) 
5 


41-7 


msfer from water to hexone 
2453 Om 241 


k. A _ Comp Cwi 
' 

O07 

Oso 2233 52.3 ooo 

Ooo 228 7-4 70 

ooo 20 wu 

20 20 OOS 

2-67 37-2. 38-3 18 

33-8 “) O05 

(b) T 

246 3-71 

250 306 1-82 

Mun 251 240 = 140 ©6277 Ta 3 

78 5-32 13 ool 

OSS 21 ool 

276 17-3 0-63 

576 see 1-8 6-73 2-73 2-73 3-38 187 2240 10 

270 18S 6-36 20 Ooo 

O29 95-3 0-08 2-27 

O22 92 O19 120 3-33 

O-28 950 O33 210 2-38 
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Table 3. Transfer of uranyl nitrate between TBP solutions and water 


Vass transfer coeff. Interfacial 


(cm. sec 10") Concentration (mg/ml) Time resistance 


Expt. 


Temp. Reynolds 
( 


©) Calculated Obsd. Bulk phase Interfacial (mins) Ri 


Water Solvent A Cop Cos st (sec 
(a) Transfer from steam stripped TV.B.P./O.K. to water 
3-24 2-3: 1) 68-2 
65-4 
63-4 
62-0 


(+) Transfer from distilled T.B.P./OK. to water 
6-06 246 25-7 574 28-1 5: 
30 20-0 


O17 


(c) Transfer from 'T.B.P./CCl, to water 
1 
TOS 


og 


10-6 


(d) Transfer from water to steam stripped T.B.P./O.K. 
B37 +70 ‘ 
1-00 
140 
1-20 


OSS 


(¢) Transfer from water to distilled T.B.P./O.K. interfacial turbulence observed 
26 34 3-80 183 5-22 170 137 640 
1-45 165 
163 


(f) Transfer from water to T.B.P./CCL, 
20 2-19 1-78 17: 11-5 171 
0-43 170 
O44 168 
O30 167 
OSS 165 
0-36 163 


(g) Transfer from water to T.B.P. /Benzene 
$-22 70 14 60 163 
0-28 193 72 
192 
0-22 191 


I 
Ooo 553 50 200 
632 25 Sooo 3740) 10 1-59 
20 2-04 
VOL. 30 1-59 
5 314 522 5-55 
1958 
627 25 MMM) 63 61-5 340 10 1-05 
20 1-54 
Oy 5-2 30 5°55 
0-50 67 680 1-54 
O17 | 67-1 5°55 
0-50 70 | 
0-33 80) 2-65 
613 25 woo 74 10 0-05 
20 0-27 
30 O45 
Ww) 
Os! 
238 Oo 0-04 
25 10 O-O4 
20 O40 
30 
O73 27-7 1-10 
20 1.76 
30 1-74 
is-0 2.78 
20-4 2-08 
22-4 80 2-22 
O45 25 77 10 O-74 
20 2-86 
3-70 
3-87 
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intersection of a line of this slope through O and 
the equilibrium curve. As the concentration 


gradients are assumed to be linear 7 must lie 
line OL. The tirst approximation is to 
assume that J lies at L so that the BP 
becomes BL with a slope of 1,. Using this value 


of i, A, is obtained and as the 


on the 
chord 


experimental 


value of A,. is known the first approximation to 
ke. can be obtained from equation (5). 
Also : 

¢.,*) 

which is the vertical distance LP’, so that by 


simple geometry the tirst position of r, Ps can 
be located on the 
P, a second value of H/, i. can be obtained from 
the slope of BP, and a second value of &,., &». 
This 
a third approximation for MH, Wey. 


distribution curs Knowing 


calculated. second approximation 


for LP and 
In practice it was found that this third approxi 


il, and 


mives a 


mation was identical with MH, and « 


Hl, were the same within the accuracy of the 
experiments. 
In the two systems uranyl nitrate — D.B.C, 


water and uranyl nitrate — hexone — water the 


equilibrium curve is tangential to the 


aqueous 
phase axis and H becomes very small for aqueous 
concentrations below 100mg ml of uranium, 
Thus < 1/k, 
accuracy 1/A 1 /k,, 


A similar series of approximations were mace 


and to a 


for transfer from water into solvent. The relevant 


concentrations are repre sented in Fie. 1 (b) and 
the appropriate distribution coetlicient, H’, is 


Initially HW’ 


a 


‘ 
was taken as the slope of BL and the same 
A third approximation 


the slope of 


method used as before. 


toH 


During an experiment O moves along a line O, 


was never necessary. 


which is equally inclined to both axes as the 
volumes of the two phases were equal. Thus a 
family of points O, I. ete. were obtained. Since, 
practice O] 
distribution curve 


however, in and practically 


coincided, and the was almost 
symmetrical between P and J the variations in 
some in- 


the position of I were negligible. In 


high degree of 


stances, e.g. at high aqueous concentrations in 


the uranyl nitrate - DBC - water system, the 
equilibrium curve was almost linear over the 
concentration range question and first 


approximation of 77 could be used with suflicient 
accuracy. 
The 


trons were 


values of the true interfacial concentra- 


obtained for each time interval and it 


was found that. although the bulk concentrations 


varied considerably, the interfacial values re 


san This variation 


maimed the to 10 per e nt. 


was considered to be too small to be significant. 


Since 1A 


resistance of 


, corre sponds to the total measured 


transfer Tk. ts the additional 


resistance present at the interface and it is con 
veniently expressed as R,, 
Values of A, AL. and obtained in 


Tables 1-3. 
set of eXyM rimental results are ob 
{ERE CE R 1366. Values 


of the interfacial resistance R have been plotted 


typical experiments are given in 
\ complet 
tainable as report 


as a function of time on Figs, 2—4. and interfacial 


5 and 6. 


concentrations on Figs. 


| | | 
Oo 0 2 30 40 50 60 70 60 
Tir 
Fic. 2. Variation of the interfacial resistance with time 


for the systems uranyl nitrate — D.B.C. — water and urany 
D.B.C. (15 N 3 N HNOg, 


from solvent to water; (+) transfer from water to solventr 


nitrate (a) transfel 


‘a 
VOL. 
195!) 
“Cy 9 
4 if C131 Cy 83) 
= | 4 
re 0 2 30 40 50 60 % 
- 
| ¢ | 
2: J 203} 
“ 
~ C136 4 34 xy G 
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Discussion 
The transfer across the interface was accom- 


panied by interfacial turbulence in some instances. 


In these particular experiments, except where the 


transfer was of uranyl nitrate from water to 
D.B.C., the observed rate was always greater than 
that calculated assuming no interfacial reaction. 
This accelerated transfer is undoubtedly due to 
the 
facial area, both of these being greater than the 
the 
the 
ellicients decreased with time, a small interfacial 


additional turbulence and increased inter- 


values used in 


Since, 


calculating individual co- 


etlicients. observed 


however, CO- 
barrier might have been present which became 
increasingly important as the interfacial turbu- 
died It appear 
reasonable to suppose that the transfer of uranyl! 


lence gradually away. would 


nitrate to D.B.C. is opposed by an interfacial 


resistance so large that it outweighs the effects 


of interfacial turbulence and increased area. 
The experimentally determined values of R, 


will of course be smaller than the true value. 


Fic, 3. 
For the 
HINO, 


water 


Variation of the interfacial resistance with time. 
systems 

D.B.C, 

interfacial 


uranyl nitrate — hexone 
(a) 


turbulence 


and 
to 
transfer 


water 
water. solvent 


(b) 


transfer from 
observed 4 


from water to solvent. 


| 
C..37 


--*—Pure TBP/ pure Ok 


R;, sec/emx10- 


Fic. 4. 
For 


Variation of the interfacial resistance with time. 
the T.BP 
transfer from (b) 


uranyl 
to 


water to solvent. 


svstems nitrate water 


(a) solvent water ; transfer from 


(a) 


BO 120 160 200 240 280 
Cy.9™mg/ mi 


Fic. 5. Variation of the interfacial resistance with inter- 
concentration ; 
water to T.B.P./OLK. 


centration ; 


facial transfer of uranyl nitrate from 
(a) effect of aqueous phase con- 


(b) effect of solvent phase concentrations 


c..48 
vA C, 62 
wie 
O et 50 oO ©) 
T me, rrur 4 
3} 
4 
, —T | 
“ 
Aor — — TAF K 
C..64 -~»—TBP/ benzene (b) 
Time,min 
NH hexone-water 2\ 
-——e—— UNH-hexone- water 50 *C 
D4 ev Time, min 
64 j 3-0} + + 10 
We 108 C..1 30 ; 
4. x +| 
2c 3 O 80 
Time “| 
UNH hexone- water 20°C Oo 40 
— INH hexone-water 5 | 
~ | [Time, min | 
—— — 148 | 
~ (b) 
0 10 20 30 40 50 60 Oo 80 
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Although in the nitric acid —D.B.C. — water 
system turbulence was only obtained during the 
transfer from D.B.C. to water, somewhat faster 
rates than calculated were obtained during 
transfer in the reverse direction. It is probable 
that these were due to experimental errors and 
are not significant. This is the only svstem 
where fast transfer rates were not accompanied 
by interfacial turbulence. 

The results generally are most conveniently 
interpreted by considering the variation of the 
interfacial resistance R, with time, concentration, 
ete. Variations of R, with time during some 
typical experiments are shown in Figs. 2-4. 
It can be seen that in many svstems the initial 
values of R. are small, indicating good agreement 
between the observed and calculated values. 
which may last as much as 30 minutes. However. 
R, then usually increased in each experiment as 
the observed transfer rate fell off with time. 
Attempts were made to correlate the values of 


with (a) the interfacial concentrations ¢, 


and (b) the bulk concentrations and 
(c) the driving foree | ¢, c,.* |, (d) the tem- 


y™ rature, (e) the Revnolds numbers of thy two 


phrase s. and (f) the time of the reaction. 


e, 
= 
x 6] ; 
| 
4 
| 
& ak 
2} 
+ 
20 30 40 50 60 70 «680 


Fic. 6. Variation of the interfacial resistance with inter 
facial concentration of the aqueous phase; transfer of 


uranyl nitrate from T.B.P./O.K. to water 


It might be thought that R. should be some 
function of ¢., or e, but, as shown in Figs. 5-7, 
there is no dependence whatsoever. The data 
for TBP OK are the most conclusive: thus in 
Fig. 5 at any time each point corresponds to a 
different experiment. In the other systems the 
experimental scatter obscures the picture some- 


what. Again, temperature did not appear to 


Lewis 
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have a very marked effect, (Figs. 2, 3), although 
there was a slight tendency for R; to decrease. 
No general relation between R. and the = six 
variables investigated was found. 

One possibility which was considered was that 


of surface contamination. Most surface active 


me, 
2 j 
t 
a ee | 
"Ff — — — | 
Cc /m 


Fie. 7. Variation of the interfacial resistance with inter 
facial concentration of the aqueous phase; transfer of 


uranyl! nitrate from D.B.C. to water. 


compounds would adsorb at the interface in a 
few seconds. but proteins, for example, might 
adsorb sufliciently slowly to give an apparently 
slowly increasing value of the interfacial resist- 
ance, Scrupulous care Was always taken to pre- 
vent contamination of the cell, ete... and both 
phases were run through an activated alumina 
column, and filtered, immediately prior to each 
experiment so as to remove any surface active 
agents present. With the exception of some of 
the T.B.P..O.K. solutions all solvents were 
fractionally distilled before use. 

From the experience gained in the earlier 
experiments in the transfer cell it was anticipated 
that unforeseen contamination would probably 
occur in a few isolated experiments giving very 
low transfer coeflicients. In a few experiments 
the transfer was indeed very slow, values of R, 
about 5— 20 times the average for that system 
being obtained. It was considered that the 
interface had become contaminated in these 
experiments so that the results were of no 
significance. In two experiments the interface 
was renewed by sucking away the liquids at each 
side after 10 minutes had elapsed and the rate 
had fallen off. However, the transfer rate 
remained unchanged and did not imerease as 


would have happened if the original surface had 


been contaminated. 
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It is not possible to compare these experiments 
directly with those of Murpocu and Prarr, and 
of Haun. The former authors used a wetted wall 
column where the duration of contact between 
the phases was very short, of the order of seconds, 
and the interface was renewed continuously. 
Consequently their results correspond to the 
initial period of this work when low values of 
R, were obtained. On the other hand, the data of 
Haun were obtained under completely quiescent 
conditions where transport to the interface was 
entirely by molecular diffusion. As this is a much 
slower process than turbulent transfer, chemical 
interfacial reactions are less evident and may 
not be detectable. 

The free energies and enthalpies of transfer of 
uranyl nitrate from water into T.B.P. and D.B.C, 
have recently been published, the values for 
T.B.P. being AG +58 and AH 3-62 and 
for D.B.C., AG 2.53 and AH 7-70, K. cals 
mole. No data are given for the transfer of nitric 
acid from aqueous to D.B.C. phases but very 
rough calculations based on the mass distribution 


coellicients give lower values. Although these 


AH values give no direct information as to the 
activation energies for transfer from water to these 
solvents they do indicate that the activation 
energies for transfer in the reverse direction must 
be greater than these values. It is perhaps 
significant that during the transfer of uranyl! 
nitrate from T.B.P. to water no instance was 
observed in which the initial value of R, was 
zero, in contrast to the transfer in the reverse 
direction when initial zero values of R, are 
obtained. In both Cases however, marked 
positive values of R; are obtained. Again, the 
absence of any appreciable interfacial barrier 
during the transfer of nitric acid between water 
and D.B.C, is in harmony with the fact that the 
AG value appears to be much less in this system. 

The question of mutual saturation of the two 
phases does not seem to be important. For 
example, the transfer of uranyl nitrate between 
water and D.B.C. results in a marked increase in 
mutual solubility unlike the corresponding 
transfer between water and T.B.P. In the former 
system the uranyl nitrate —- D.B.C. complex con- 


tains four molecules of water whereas the uranyl! 
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nitrate — T.B.P. complex contains no water and 
is unstable in water. Nevertheless, similar R, 


versus time, graphs are obtained in both cases, 


CONCLUSIONS 

1. Overall transfer coefficients were deter- 
mined for the transfer of uranyl! nitrate between 
water and the three solvents, dibutyl carbitol, 
methyl iso-butyl ketone, and tri butyl phos 
phate (which was dissolved in odourless kerosene, 
carbon tetrachloride, or benzene). The overall 
transfer coetlicients were also determined for the 
transfer of nitric acid between water and dibutyl 
earbitol, Experiments were carried out at 
various temperatures, 20°, 25° and 50°C. The 
experimental coeflicients were compared with 
those calculated using the correlation derived 
in Part I and assuming the absence of any slow 
step at the interface itself. 

2. Greater transfer coeflicients than predicted 
were obtained in a number of instances and. 
except for the transfer of nitric acid from water 
to D.B.C., all these were associated with the 
phenomenon of interfacial turbulence. The nitric 
acid coefficients were only slightly greater than 
predicted, 

3. In the majority of these experiments 
where interfacial turbulence was absent good 
agreement was obtained between the calculated 
transfer coetlicients and the initial values of the 
observed coetlicients. During the course of each 
experiment however, the values of the observed 
coeflicients diminished, corresponding to the 
build-up of an interfacial resistance, R.. 

+. The initial period during which there was 
good agreement between the calculated and 
observed coeflicients, i.e., when was negligible, 
varied considerably but was of the order of several 
minutes. These times are large compared with 
the residence times in most types of liquid 
liquid contactors so that the build-up of resistance 
afterwards is probably not of practical import- 
ance. Further, it appears probable that the 
increase in R, is associated with the age of the 
interface so that the effect will be absent in 
contactors, such as packed columns, where drops 
are continually coalescing and breaking up. 


5. All attempts to correlate the interfacial 


. 
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resistance R,; with the possible variables in the ’ = distribution coefficients defined by equation (5) 


systems failed. Nevertheless it is a real phenom- — transfer rate, g/em* see 


theoretical overall mass transfer coefficient 
enon and not due to surface adsorption of im- tant be @ — a 
defined by equation (6), em /sec 

purities. It is possible that some sort of time lag observed “ overall” mass transfer coefficient 
occurs in the build up or breakdown of the based on solvent phase, em /sec 

uranyl! nitrate-solvent complexes formed. individual mass transfer coeflicient, em 

velocity constant for first order interfacial 

ichnowledgements—Most of the experimental work 

was carried out by Messrs. DD. Beawes, 1D). Garay and 

J.D. C. My thanks are due to Dr. H.R. 


Peart for much helpful discussion and to Mr. AS. Warr 


reaction based on solvent phase, em /sec 
length of stirrer, em 

stirring rate, revolutions per sec 
interfacial resistance 1/k,,, see 
for his interest in this work. time. sec 
Volume of cell compartments, 
NOTATION Reynolds number = L* N/v. 
area of interface in transfer cell, em* kinematic viscosity, em* /sec 
concentration of uranium g (1000 cm (ie. mg/ml) Subscripts 

concentration in equilibrium with the other s = solvent phase 


phase s 


; = value in solvent phase at interface 
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Letters to the Editors 


(Received 18 November 1957) 


Note in connexion with the paper “ Boundary’ which we gave reference. We are glad therefore to note 
conditions of flow reactor *’.* the reference at this time. DamKOuter divided a plug 
IN connexion with our paper “* Boundary Conditions of flow reactor into fore section, reaction section and after 
Flow Reactor” it was drawn to our attention by Dr. section as we did. However, axial diffusion was permitted 
D. W. van Krevenen and Mr. H. A. G. Curr of the to vary between sections in our work, which was prompted 
by a classroom discussion of the logic for the choice of 


Staatsmijnen in Limburg, Netherlands that the paper 
by G. DamKounen, Z. Elektrochem. 1937 43, 1 dealt boundary conditions at the reaction section exit positiou. 


with the same general topic. We were not aware of this J. F. Wenner 
work, as apparently was also the case in the studies to R. H. Witue tm 


Department of Chemical Engineering, Princeton University. Princeton, New Jerser 
} - 


*Wenoner J. F. and Witnetm R. H., Chem. Engng. Sci., 1956, 6, 89. 
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Even the advent of atomic energy, in its peaceful applica 


will not soon dishodwe the il fire from 


hearth Ihe 


thoms at any rate 
the Koglishman’s smokeless 


fuel for 


provision « 


this fire, and for domestic purposes generally, 


main British interest in the researches described in 
This tith 


rs the 
“ Chemical Engineering in the Coal Industry 
wis ¢ hosen for the 


British 


report of the papers and discussions 
at an 
organised by the National Coal 
Research Estalbh 


in June 1056 


presented by French and German 
International Conference 
Board held at its Coal 


Stoke 


and shment at 


Orchard, near Cheltenham 
The production of pitch-bound briquettes of inthracite 


by the Phurnacite process is well known but anthracite 


umd «a large inerease in the supply 
solid 
processing the more 


higher heat ile 
partial dew 


is relatively 


smokeless fuel will require ecomonu 


of reactive 


methods of readily available non 


coking coals of ntent stayes 


mvolved are perature 


carbonisation, briquettin und in some cuses a mild 
= 


oxtdation to reduce the ikhing and swelling powers of 


the coal This last we mav be necessary because the 


caking power of many coals, though too bow for thre 


manufacture of coke by established methods. ti neverthe 


less suflicient to cause briquettes to agglomerate, and 


excessive swelling leads to their distortion 


hi « 


tion inte a satisfactory process presents muinmy pr 


f these stages has its difficulties and their mtegra 


This 
in 


The complexity of coal and the need for « 


process costs add to the difficulties The failure of many 


semicarbonization processes in the past ha nphasised 
the defects of the work 


effort 


empirical approach 


deseribed rests on a soumeer foundation 


has been devoted to the determination of necessary data 


amd modern techniques are being tested In particular 


attention is devoted to the use of fluidization both for the 
oxidation and stages 

rhe first paper, by A. I 
laborators 
G. 1 


per har 


Rover. ws devoted inly to 
the cond by 


i ile of 1 ton 


work on oxtdation 


Jenkins, to fhuidized oxidation on 


Production of briquettes is greater in Germany than in 
this « 


af ttes desorbed W 


country, and German experience in the irbontzation 


Reerink., In nteyrated 


wrocess, the briquetting stage may be mterposed between 
| 


two heat treatments, « oxidation and car nization 


ami DD ll 


to economies in heat and in pitch binder 


Grecory describes how hot briquetting leads 


A. Pevravy and P. Focu begin their paper by consider 


ing possible methods of heating coal fines im semi 


carbonization. External heating of fluidized beds 


lead to difficulties in scaling up, whereas internal heating 
except by recirculation of the gas evolved, dilutes the gas 
and reduces its valuc In spite of this, internal heating by 


air blast was chosen for an industrial plant, part of the 

heat being provided by partial combustion of the charge 
The last by G. H 

\. and by J. L. Sanarier 


produced in low-temperature carbonization 


Warson and 


deal with the 


two 


turs 
Summing up, the results of sone research are recorded 
areal the 


book In vddlition the re is 


in a convenient form discussion of each paper 


enhances the value of the 
spared the labour of translation which is usually necessary 
Netivity on the 


in following current technology abroad 


field discussed is considerable, and it is inevitable 
book of 


further che vek pinent of sore Processes and the 


this nature is superseded fairly rapidly by thy 


abandon 


ment of others 


M. J. G. Witson 


Ss. D. Hamann: Physico-Chemical Effects of Pres- 


sure. Butterworths Scientific Publications, London, 1057. 


pp is 


Tur study of the physical and chemical effects of high 


pressure has long been pursued in a number of universities, 


research institutes and, more recently in industrial 


research laboratories. The author, Dr. Hamann, is the 
leader of one such laboratory 


Screntittc 


which is part of the Australian 


Commonwealth and Industrial Research 


Organization's Division of Industrial Chemistry and is 


situated in the department of Chemical Engineering in 


the University of Sydney In this monograph he en- 
deavours to bridge the gap between physics and applied 
chemistry by surveying the effects of pressure in the field 
of physical chemistry. Most of the book is concerned with 


changes brought about by hundreds or thousands of 
atmospheres, but some data at lower pressures have been 
included in instances where the experimental pressure 
range was restricted or where the pressure effects were 
unusually large 

After 


pressure experimentation in a chapter which is influenced 


a general introduction, the author discusses high 


by much of his own practical experience. The designs can 
be used to develop apparatus working almost at the 
breaking point, to explore the extremes of measurements 
devoted to: the 


design of single and multiple cylinders, closures, packings, 


at very high pressures. Sections ar 


optical windows, electrical connections; also to transient 
pre SSsuUTes, shock waves and pressure measurement 


The compression of gases, liquids and solids in single 


and two component systems is surveyed in Chap. 3, 


810 


SS 


while Chap. 4 deals with phase transitions. Chap. 5 
discusses transport properties : viscosity, diffusion, heat 
conduction and sound conduction and Chap. 6 reviews 
dielectric and optical properties. These four chapters 
survey the changes in the physical properties of matter 
brought about by pressure. In the subsequent pages the 
author shows how these can be used to interpret the 
pressure effects observed in electrochemistry, chemical 
equilibria and kinetics of reaction. 

The considerable recent developments in these subjects, 
many originating from the authors own work. are here 
fully presented, together with a critical review of earlier 
relevant literature. It would perhaps have been useful 
without noticeably increasing the size of the book. to list 
more fully the experimental work and the magnitude of 
the observed pressure effects in each case. A chapter on 
miscellaneous effects, the most spectacular of which is the 
high pressure formation of diamonds (the history of which 
is not without a littl wry humour) concludes the text 
Appendices tabulate the sources of published literature 
on liquid compressibilities and the effect of pressures on 
the melting points of solids. A commendable feature is 
the extensive bibliography which is at the end of the book 
and so saves all cross-referencing. 

\s it is lucidly written and well produced, it is a pleasure 
to read this book. It is useful reading for those intending 
to take up research in this field, as a number of limitations 
to our present knowledge are pointed out, and prosst bole 
future lines of research are indicated, With the increasing 
commercial importance of high pressure processes, this 
volume presents the theoretical background for the 
chemical engineer wishing to « xploit this new technology 
For the experimental physicist and chemist studying 
high pressure phenomena the book will prove invaluable, 
and will take its place on his shelf alongside BripGMan’s 
classic The Physics of High Pressure 


W. Srrauss 


Messen und Regeln in der Chemischen Technik. 
lirsy. v. J. B. up. O. Winkier 
Mit 1331 Abb. Springer, Berlin u.s.w., 1057. NIX 1261 
S.. DM 136, 


der steigenden Bedeutung der Mess- und 
Regelungstechnik in der chemischen Lndustrie wurden 
auch an die hiermit betrauten Ingenieure und Physiker 


immer hOhere Anforderungen an Wissen und Kénnen 


yestellt. Da cine zusammenfassende Darstellung dieses 


Book reviews 


umfangreichen Gebietes bisher fehlte. haben sich die 
Herausgeber zu diesem Sammelwerk entschlossen. das 
zugleicn als Gemeinschaftswerk der NAMUR (Normen 
arbeitsgemeinschaft fiir Mess- und Regeltechnik in der 
chemischen Industrie) erschienen ist. 

Der Inhalt dieses iiber 1200 Seiten starken Bandes sei 
durch die Titel der Hauptkapitel angedeutet : Tempera- 


turmessung; Mengen-, Durchfluss- und Standmessung: 
Druck und Differenzdruckmessuny: Physikalische 
\nalvsenverfahren: Mlekt rochemiselhe Messmethoden: 


Speziclh Bet riebsmessverfahren: Revelungstechnik: 
Vorsehriften und Hilfsmittel fir Planung und Betrieb 
von Mess- und Revgelanlegan: zur Organisation der 
Bet riebskont roll 

Die 31 Mitarbeiter stammen ganz tiberw ievend aus der 
chemischen Industrie, wo sie eivene praktische Erfah- 
rungen tiber das Objekt ihrer Darstellung gvesammelt 
haben. So enthalten die einzelnen Abschnitte meist einen 
kurzen Abriss iiber die Theorie des Verfahrens. an den 
sich dann die Beschreibung und kritisehe Wiirdigung det 
vorhandenen Geriite anschliesst. Je nach Autor und 
Gegenstand finden sich hier naturgemiss Unterschiede. 
doch ist die Darstellung immer bewusst unkompliziert 
yehalten und auf das Wesentliche verichtet 

Das Buch wendet sich weder an den Theoretiker noch an 
den Hersteller, sondern an den Benutzer von Mess- und 
Reygelgeriiten. So herrscht auch der Gesic htspunkt det 
praktischen Bewihrung im Betriebe vor. wenn auch dix 
Grenz zur Laboratoriumstechnik yclegentlich zu 
liberschreiten war. Ein solehes Werk ist nur durch 
Zusammenarbeit zwischen Ingenieur und  Phvysiker 
moOglich. Unterschiede in der Handhabung der Massvstem« 
wirken dabei stOrend, aber es ist wohl noch zu friih. um 
auf diesem Gebiet Kinheitlichkeit zu erreichen. Hervor- 
zuheben ist noch, dass das Buch auch die kontinuier- 
lichen physikalischen Analvsenverfahren und die Par 
Messung ausfiihrlich abhandelt. In den oben venannten 
Hauptkapiteln sind einige Monographien enthalten, dic 
auch der Fernerstehende mit grossem Gewinn studic ren 
wird, wie denn tiberhaupt der klare und schlichte Stil sehr 
wohltuend auf den Leser einwirkt. Die dussere Ausstat- 
tung des Buches ist vorziivlich 

Den Herausgebern und ihren Mitarbeitern vebuhrt 
uneingeschriinkte Anerkennung, dass sie dieses miihe volk 
Unternehmen zum Erfolg brachten Der Zweck des 
Sammelbandes diirfte vollstiindix erreicht sein. dem 
Fachmann ein Nachschlagewerk und dem Lernenden ein: 
Kinfiihrung vegeben zu haben 
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Papers to be published in future issues 


S. Karz \ statistical analysis of certain mixing problems 


Heat transfer between a fluidized bed and a horizontal tulx 


Hi. A. Va : 


Adiabetic convection batch drving with recireulation of air 


DD. A. Mews 


Daksminamurry, G. Javanama Rao, M. C. Venkata Rao: Ternary vapour liquid equilibria system 


Benzene, cyclohexane methviisobutyvl ketone 


K. Kuamer, VAN and van : Studies on ion-exchange IV. Kinetic theory of ion-exchange 


J.C. H. Lanssen and D. W. van Kerevecen: Studies on ion-exchange V. Kinetic measurements 


KLAMER 


\. P. Jn, and C. A. Puank : Perforated-plate performance 


Rasa Rao and C, Venkata Rao: Flooding rates in packed liquid extraction towers 


Bexnxnerr: Diffusion in binary mixtures 


Kiamer and DD. W. van cen : Studies on ion-exchange Vi. Design and sealing up of ion-exchange columms 


D. B. Sraupine : A note on mean residence-times in steady flows of arbitrary complexity 


R. A. Sener : The thermal conductivity of catalyst particles 


Erratum 


The kinetic implications of an empirically fitted yield surface for the vapour phase oxidation of 


naphthalene to phthalic anhydride Chem, Engng. Sci, 1057 © 105 


The author regrets that errors were made in the integration in both Case 1 and Case 3 on p. 108, The equations should 


read as follows 


Case 1 


Disappearance of Naphthalene half order. Sulbsequent steps first order 


ky ky (hy + hy) ky thy Mot 
| 2h, (ky ky) 2h. k. 


Case 3 


Disappearance of naphthalene second order, Subsequent steps first order. (In this case integration leads to an intinits 


rapidly converging series. Only terms contributing more than 1 per cent to the total solution under normal operating 


conditions have been included) 
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